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Abstract 



We investigate the strong coupling limit of quantum chromodynamics on a lattice (lattice 
QCD) for systems with high baryon density. Our starting point is the Hamiltonian of 
lattice QCD for quarks with N c colors and Nf flavors at strong coupling. We regard 
this Hamiltonian as an effective Hamiltonian that describes QCD at large distances, and 
expand it in orders of the inverse coupling. In leading order one finds a Hamiltonian that 
looks like an antiferromagnet. It has interactions between nearest-neighbor and next- 
nearest-neighbor sites, and is invariant under U{Nf) x U (Nf). Physically this Hamiltonian 
describes meson dynamics with a fixed background of baryon density. Our goal is to 
extract the ground state and low lying excitations of this Hamiltonian for systems with 
non-zero baryon density. 

We first write the partition function of the Hamiltonian with boson coherent states and 
take the limit of large N c and Nf with the ratio r = N c /Nf kept fixed. In this limit 
we make a mean field ansatz to extract the ground state. For zero density we extend a 
condensed matter physics result to three dimensions, and show that above a critical value 
of r the ground state breaks global symmetry. Departing from zero density we show that 
for certain cases, where the baryon density is periodic, the critical value of r increases 
with increasing density. In cases where the baryon positions are random, we average over 
them (in a quenched calculation that is valid only for low density) and show that the 
mean field ansatz is misleading near the critical r. Nevertheless, following the condensed 
matter analog of this quenched calculation, we expect that in the presence of baryons, 
low lying excitations will have non-zero widths and reduced velocities. 

Writing the partition function with generalized spin coherent states leads to a nonlinear 
sigma model, which we investigate in the rest of the thesis. We first analyze the classi- 
cal counterpart of a toy sigma model with mean field theory, and find that the critical 
temperature, under which the symmetry of the model is spontaneously broken, decreases 
with increasing density. 

The quantum sigma model is treated in the limit of large N c with Nf fixed. We find that 
the ground state is classically degenerate. For zero density the degeneracy is discrete. 
In the case of non-zero uniform density the ground state is continuously and locally 
degenerate, like classical ground states of frustrated magnetic systems. 

When we consider quantum corrections of order l/N c we find the phenomenon of "order 
from disorder," also present in frustrated magnetic systems. For zero density we extend 
an old result to the next-nearest-neighbor theory, and show that fluctuations remove the 



discrete degeneracy to choose a ground state that spontaneously breaks U (Nf) x U (Nf) 
to U(Nf)y. For non-zero density we see that fluctuations remove the local degeneracy 
and choose a state that spontaneously breaks U(Nf) x U(Nf) as well as discrete lattice 
rotations in a way that depends on Nf and the baryon density. 

Finally we find a rich variety of low lying excitations that includes type I and type II 
Goldstone bosons. The latter emerge only for non-zero density and only when quantum 
fluctuations are taken into account. Their energies are of order l/N c , and are quadratic 
in momentum. Bosons of either type can develop anisotropic dispersion relations. 



Contents 



1 Introduction 1 

2 The effective Hamiltonian 9 

2.1 Hamiltonian approach to lattice QCD 9 

2.2 The strong coupling limit 15 

2.3 Strong coupling perturbation theory 17 

2.3.1 Second order: the antiferromagnet 18 

2.3.2 Third order: the baryon kinetic term 24 

2.4 Summary 26 

3 Mean field theory with Schwinger bosons 29 

3.1 Schwinger bosons 30 

3.2 Hubbard-Stratonovich transformation 34 

3.3 Mean field theory at zero density 36 

3.4 Mean field theory at non-zero density 41 

3.4.1 Periodic density with baryon concentration c = |, | 42 

3.4.2 Random baryon positions - the quenched approximation 47 



ii CONTENTS 



4 Nonlinear sigma model and the classical ground state 63 

4.1 Coherent states and sigma fields 63 

4.2 Partition function and action 67 

4.3 The classical counterpart of the sigma model 69 

4.3.1 Mean field theory 70 

4.3.2 The toy model of U(3) 72 

4.4 The semiclassical large N c limit 80 

4.4.1 Zero density 81 

4.4.2 Non-zero density 83 

5 Quantum fluctuations in the sigma model 89 

5.1 Removal of classical degeneracy at zero density 89 

5.2 Removal of classical degeneracy at non-zero density 92 

5.2.1 Zero modes at 0(1) 95 

5.2.2 Removal of zero modes at 0(1/N C ) - the self energy calculation . . 100 

6 Next-nearest-neighbor interactions in the sigma model 107 

6.1 Action and symmetries of next-nearest-neighbor interactions 108 

6.2 Ground state of the next-nearest-neighbor theory 109 

6.3 Excitations in the next-nearest-neighbor theory 116 

7 Summary and discussion 125 
A Anticommutation relations of baryon operators 131 
B Discretization of Euclidean time 135 
C Minimizing the next-nearest-neighbor term 139 



CONTENTS 



iii 



D The function $ for the U(3) toy model at non-zero density 141 
E The matrices C and D in the 0(1/ N c ) calculations 145 
F Feynman rules for the effective action of the sigma model 147 

G Integral equations for the self energy of the zero modes 151 

G.l Evaluation of diagrams 152 

G.2 Evaluation of frequency integrals 154 

G. 3 Solution of integral equations 158 

H Hamiltonian approach to order from disorder 163 

H. l The ground state at 0(1) 166 

H.2 Self energy and 0(1/N C ) effective Hamiltonian 167 

I Effective interaction between odd sites in the sigma model 175 

J Fluctuations on the odd sites at non-zero density 181 

K Coefficients in the dispersions of Goldstone bosons 183 

L Classification of fields according to unbroken subgroups 185 

Acknowledgments 197 



List of Figures 

1.1 The representation of U(AN f ) 5 

2.1 U(ANf) representation of a baryon operator with N c = 3 21 

2.2 SU(N C ) representation carried of a single site state 21 

2.3 S n representation of a single site state 22 

2.4 U(ANf) representation of a single site state 23 

3.1 The integration measure of 7 for d = 3 39 

3.2 f{r]) - right hand side of MF equation at zero baryon number 40 

3.3 Phase diagram of Schwinger bosons MF theory at T = and zero density. . 41 

3.4 A unit cell with c = | concentration of "magnetic" vacancies 43 

3.5 The integration measures of 71,3,5 for c = | 45 

3.6 Right hand side of the MF equation for c = | 46 

3.7 The integration measures of 71,3,5 for c = | 48 

3.8 Right hand side of the MF equation for c = \ 49 

3.9 The boson propagator 51 

3.10 The impurity vertex 52 

3.11 Scattering of bosons off three impurities 52 

3.12 The loop diagram that corresponds to the expression (3.93) 53 

v 



vi 



LIST OF FIGURES 



3.13 The loop diagrams that corresponds to p = 4, q = q e = 2, and s = 1, 2, 3. . 54 

3.14 Right hand side of the MF equation with c ^ in the quenched approxi- 
mation 58 

3.15 MF action for c = |; quenched approximation vs. exact diagonalization . . 59 

3.16 MF action for c = |; quenched approximation vs. exact diagonalization . . 60 

3.17 Right hand side of MF equation for c = |; quenched approximation vs. 
exact diagonalization 61 

3.18 Right hand side of MF equation for c = |; quenched approximation vs. 
exact diagonalization 62 

4.1 The most symmetric 5 = state in the U(3) toy model 72 

4.2 Adding a baryon in the U(3) toy model 73 

4.3 $ at zero density for N = 3 76 

4.4 $ non-zero density for iV = 3 79 

4.5 Classical phase diagram of the U(3) toy model 80 

4.6 A configuration with zero baryon number on all sites 81 

4.7 A configuration with alternating baryon number on adjacent sites 81 

4.8 Adding one baryon to the lattice 83 

4.9 Two sites that carry U(N) representations with mi and m 2 rows 84 

5.1 Coupling two 7-row reps of U(12) 93 

5.2 Weight diagram for iV = 3, m = 1, iV c = 2 94 

5.3 Schwinger- Dyson equation for the x fields self energy 101 

5.4 Solution of integral equation 103 

5.5 Dispersion relation of the x fields 104 

F.l Propagator of the 7r fields 147 

F.2 Propagator of the x fields 147 



LIST OF FIGURES vii 



F.3 Vertices corresponding to the quartic interaction 148 

F. 4 Vertices corresponding to the cubic interaction 149 

G. l Evaluation of Si 152 

G.2 Evaluation of S 2 153 

G.3 The first Brillouin zone (BZ) of the fee lattice 159 

G.4 Behavior of the approximated integral equation 160 



List of Tables 



4.1 Baryon number configurations in the U(3) toy model 77 

5.1 Field classification under the unbroken group of the nearest-neighbor theory. 100 

6.1 Realization of chiral symmetry for different densities 114 

6.2 Eigenvalues of n e o and m e o 121 

K.l The coefficients in the dispersion relations of Goldstone bosons 184 

L.l Classification of a fields on the even sites 186 

L.2 Classification of a fields on the odd sites 187 



ix 



Chapter 1 
Introduction 



The theoretical framework that describes the strong interactions of quarks and gluons is 
the theory of quantum chromodynamics (QCD). It is of fundamental interest to study 
what QCD predicts for nuclear matter at extreme conditions. This includes systems at 
very high temperatures (~ 10 12 K) accessible in ultra-relativistic collisions of heavy ions, 
and systems with high nuclear densities (~ 10 15 -10 16 g cm -3 ) that can be found in the 
cores of neutron stars. In this work we are concerned with the latter, nuclear systems of 
high density. 

The study of quantum chromodynamics at high density is almost as old as the theory 
itself, and goes hand in hand with the discovery of asymptotic freedom [1, 2]. The first 
to approach the problem were Collins and Perry [3]. They argued that at high baryon 
densities of the order of ten to a hundred times nuclear density, the hadrons overlap 
and lose their individuality. At these densities asymptotic freedom weakens the color 
interactions, and in the absence of long range forces the system can be considered as a 
weakly interacting "quark soup" , with quarks arranged in a Fermi sea. 

Next, Barrois [4, 5] suggested that the Fermi quark sea is unstable, and shows some 
phenomenon analogous to superconductivity. He suggested that the instability is caused 
by gluon exchange between quarks in the antisymmetric 3 channel, which is attractive. 
Nevertheless Barrois argued in [4] that the ground state does not break SU(3) gauge 
symmetry, and is rather characterized by condensates of six quark operators, which are 
color invariants. 

The next step is the work of Bailin and Love [6], that studied relativistic fermionic sys- 
tems in similar lines to the study of fermions in condensed matter systems that exhibit 
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non-perturbative phenomena. Using similar arguments to those in [3] they justify that 
quarks in neutron stars weakly interact and that they are the relevant degrees of freedom, 
rather than the neutrons themselves. They describe the superfluid, superconductor, and 
color superconductor ground states of this dense quark matter, generated by one-gluon 
exchange between quarks in the 3 channel. In contrast to Barrois they treat condensates 
of quark-quark cooper pairs, rather than higher condensates. These condensates are not 
gauge invariant and spontaneously break S77(3) co ior- This is the phenomenon of color 
superconductivity (CSC). 

Although these are fascinating theoretical predictions, one finds that the the gap and 
critical superconducting temperature T c are in the MeV range. This is disappointing, since 
it means that the color superconductive part of the QCD phase diagram is very narrow. 
For that reason the study of color superconductivity was at a halt for over a decade. The 
stimulus for the revival of the idea of CSC was the observation [7, 8] that the instanton- 
induced quark-quark interaction can be much stronger than that induced by simple one- 
gluon exchange, and can thus give a transition temperature on the order of 100 MeV. 
Subsequent work [9] showed that the perturbative color-magnetic interaction also gives 
rise to a strong pairing interaction. These and other dynamical considerations [10] underlie 
a picture of the ground state of high-density QCD in which the SU (3) gauge symmetry 
is spontaneously broken by a BCS-like condensate. The details of the breaking, which 
include both the Higgs (or Meissner) effect and the rearrangement of global symmetries 
and Goldstone bosons, depend on quark masses, chemical potentials, and temperature. 
Prominent in the list of possibilities are those of color-flavor locking in three-flavor QCD 
[11] and crystalline superconductivity — with broken translation invariance — when there 
are two flavors with different densities [12]. For a review see [13]. 

As noted, CSC at high density is so far a prediction of weak-coupling analysis. One expects 
the coupling to become weak only at high densities, and in fact it turns out that reliable 
calculations demand extremely high densities [14]. The use of weak-coupling methods to 
make predictions for moderate densities is thus not an application of QCD, but of a model 
based on it. It is imperative to confirm these predictions by non-perturbative methods. 
Standard lattice Monte Carlo methods, unfortunately, fall afoul of well-known technical 
problems when the chemical potential is made non-zero, although we do note remarkable 
progress made recently in the small-/i regime [15, 16, 17]. For a review of these methods 
see for example [18] . 

In this work we study high-density quark matter based on lattice QCD in the strong- 
coupling limit, which we regard as an effective theory that describes QCD at large dis- 
tances. We work in the Hamiltonian formalism, which is more amenable than the Eu- 
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clidean formalism to strong-coupling perturbation theory and to qualitative study of the 
ensuing effective theory [19, 20, 21, 22]. The fermion kinetic Hamiltonian is a perturba- 
tion that mixes the zero-flux states that are the ground-state sector of the electric term 
in the gauge Hamiltonian. In second order, it moves color-singlet fermion pairs around 
the lattice; the effective Hamiltonian for these pairs is a generalized antiferromagnet, with 
spin operators constructed of fermion bilinears. 

We depart from studies of the vacuum by allowing a background baryon density, which is 
perforce static in second order in perturbation theory. Our aim is to discover the ground 
state of the theory with this background. In third order (when N c = 3) the baryons 
become dynamical; we display the effective Hamiltonian but make no attempt to treat it. 

The symmetry group of the effective antiferromagnet is the same as the global symmetry 
group of the original gauge theory. This depends on the formulation chosen for the lattice 
fermions. Following [20], we begin with naive, nearest-neighbor fermions, which suffer 
from species doubling [23] and possess a global U (4iV/) symmetry group that contains the 
ordinary chiral symmetries [as well as the axial U (1)] as subgroups. We subsequently break 
the too-large symmetry group with next-nearest-neighbor (NNN) couplings along the axes 
in the fermion hopping Hamiltonian. A glance at the menu of fermion formulations reveals 
the reasons for our choice. Wilson fermions [24] have no chiral symmetry and make 
comparison of results to continuum CSC difficult if not impossible. Staggered fermions 
[25] likewise possess only a reduced axial symmetry while suffering a reduced doubling 
problem. The overlap action [26] is non-local in time and hence possesses no Hamiltonian; 
attempts [27] to construct an overlap Hamiltonian directly have not borne fruit. Finally, 
domain- wall fermions [28, 29] have been shown [22] to lose chiral symmetry and regain 
doubling when the coupling is strong. 

While the NNN theory still exhibits doubling in the free fermion spectrum, we are not 
interested in the perturbative fermion propagator but in the spectrum of the confining 
theory. We take it as a positive sign that the unbroken symmetry is now U{Nf) x U{Nf). 
1 This symmetry is what we want for the continuum theory, except for the axial U{\). 
The latter can still be broken by hand [32]. 

Our emphasis on the global symmetries is a consequence of the fact that the gauge field 
is not present in the ground-state sector and does not reappear in strong-coupling pertur- 
bation theory. In other words, confinement is a kinematic feature of the theory, leaving 
no possibility of seeing the Higgs-Meissner effect directly. This is but an instance of 

1 The breaking of the naive fermions' symmetry by longer-range terms is a feature [20] of SLAC fermions 
[30] and also occurs if naive fermions are placed on a bcc lattice [31]. 
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confinement-Higgs duality, typical of gauge theories with matter fields in the fundamen- 
tal representation [33]. Our aim is thus to identify the pattern of spontaneous breaking 
of global symmetries. For various values of N c and Nf, this can be compared to weak- 
coupling results [34]. 

Other groups have recently studied the strong-coupling limit of QCD in the Hamiltonian 
and Euclidean formalisms at non-zero chemical potential and baryon density [35, 36, 37, 
38, 39, 40, 41]. We differ from most approaches in eschewing mean field theory in favor of 
the exact transformation to a different representations, which is amenable to semiclassical 
treatment. Also we base our program on NNN fermions; we also work at fixed baryon 
density. 

The outline of this dissertation is as follows. In Chapter 2 we first give a brief introduction 
to Hamiltonian lattice QCD, and discuss our approach that regards the strong coupling 
limit of the Hamiltonian as an effective model for QCD at large distances. Next, in 
Section 2.3, we present the derivation of the effective Hamiltonian of lattice QCD in 
strong-coupling perturbation theory [20, 21]. The second-order Hamiltonian [0(l/g 2 )] 
is an antiferromagnet with U(ANf) spins; the global symmetry group is U(ANf) for the 
nearest-neighbor theory, broken to U(Nf) x U(Nf) by NNN terms. The baryon number 
at each site determines the representation of U(ANf) carried by the spin at that site. In 
second order, baryon number is static; it becomes mobile in the next order, where (for 
N c = 3) the new term in the effective Hamiltonian is a baryon hopping term. 

In the remainder of this thesis, we work only to 0(l/g 2 ), where the baryons are fixed in 
position. Motivated by the similarity of our Hamiltonian to the Heisenberg antiferromag- 
net, we apply condensed matter methods developed for that problem. Indeed, condensed 
matter physicists have generalized the SU(2), spin-1/2 Heisenberg model to SU(N) in 
many representations [42, 43, 44, 45, 46, 47, 48, 49, 50], which corresponds to adding fla- 
vor and color degrees of freedom to the electrons. 2 These are exactly the generalizations 
needed for our effective Hamiltonian. With iV c colors and N (single-component) flavors, 
a site of the lattice can be constrained to contain a color-singlet combination of mN c 
particles. The flavor indices of the spin then make up a representation of SU (N) whose 
Young diagram has N c columns and m rows (see Fig. 1.1). We set 

N = AN f (1.1) 

and the correspondence is complete (until we include NNN terms in the Hamiltonian). 

2 We refer the reader to the paper by Read and Sachdev [46] for a survey, including a phase diagram 
in the (N,N C ) plane. 
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Figure 1.1: The representation of U(N) carried by the spin in the effective antiferromagnet 
[See Chapter 2]. m is related to the baryon number B at the site according torn — B+2Nf, 
with \B\ < 2N f . 



In Chapter 3 we analyze the effective Hamiltonian by following Arovas and Auerbach [45] 
and write it with boson operators that obey appropriate constraints to restrict quantum 
states to the Hilbert space represented by Fig. 1.1. After introducing the boson operators 
and their corresponding path integral in Section 3.1, we apply the Hubbard- Stratonovich 
transformation in Section 3.2 and decouple the quartic boson interaction to a quadratic 
interaction. The result of integrating the bosons is a path integral tractable in the limit of 
large N c and Nf with a fixed ratio r = N c /Nf. We take this limit in Sections 3.3-3.4 and 
study a mean field ansatz for the ground state of a set of baryon number configurations. 
In Section 3.3 we extend the zero density results of [45] from dimension d = 2 to d = 3. 
As for d = 2, one finds a Bose-Einstein condensate (BEC) above some critical value of r. 
This BEC corresponds to a Neel ground state which represents a spontaneous breaking of 
chiral symmetry as expected. For QCD, with N c = 3, the system is safely in the ordered 
phase for any reasonable number of flavors. In Section 3.4 we depart zero density by 
inserting non-zero baryon number on some of the sites. More precisely we fill some of the 
sites with the maximum allowed baryon number per site. The U(N) representation on 
these sites then has m = N rows, and is a singlet. Analyzing again a mean field ansatz 
we discuss how these singlets affect the free energy and critical value of r. In this chapter 
we mainly restrict to the disordered side of the phase transition. 

In Chapter 4, and in the rest of the thesis, we treat the system deep in its ordered phase. 
Also the non-zero baryon number configuration we investigate are different than those of 
Chapter 3, and we mainly concentrate on uniform baryon configuration, with the same 
baryon number B on all sites. In Section 4.1-4.2 we follow Read and Sachdev [46] and 
employ spin coherent states [51] to write the partition function as a path integral for a 
nonlinear sigma model. Nf and m determine the target space of the sigma model to be the 
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symmetric space U(N)/[U(m) x U(N — m)\\ the number of colors N c becomes an overall 
coefficient of the action. As for the quantum Hamiltonian, the nearest-neighbor theory 
is symmetric under U(N) while the NNN terms break the symmetry to U (Nf) x U {Nf) 
(while leaving the manifold unchanged). 

In Section 4.3 we make a pause from the real sigma model, and treat its classical counter- 
part, with a variety of baryon number configurations. This analysis presumably describes 
the high temperature region of the system, where quantum fluctuations are not very im- 
portant. It also serves as a tool to see how the symmetry realization changes as we add 
baryons to the lattice. Although the classical sigma model can be treated numerically (its 
action is real), we postpone this to future research and use mean field theory instead for 
a toy model that has N = 3. In that case calculations can be done analytically, and we 
extract the classical phase diagram in temperature-density plane. 

Next we return to discuss the real quantum sigma model. The factor of N c that multiplies 
the action invites a large- N c analysis which we make in Section 4.4. We start by studying 
this limit for zero baryon number in Section 4.4.1. This leads us to an exercise proposed 
and solved by Smit [21], in generalizing the zero baryon number sector to allow baryon 
number ±B on alternating sites; this means specifying conjugate representations of U(N) 
on alternating sites, with respectively m and N — m rows. We find the ground state of 
the nearest-neighbor theory to be discretely degenerate as seen in [21]. We extend this 
result to the NNN theory and show how the degeneracy is removed by the NNN terms. 

We turn to non-zero baryon density in Section 4.4.2. As mentioned, we study homogeneous 
states, in which all sites carry the same representation of U(N), with m > N/2. We begin 
by studying the two-site problem, and we learn that its ground state is continuously 
degenerate. This degeneracy is additional to the global degeneracy of rotating the sigma 
fields on the two sites together. Replicating this ground state to the infinite lattice makes 
the lattice vacuum locally degenerate. This makes our system similar to some frustrated 
models of magnetic systems [52, 53, 54, 55, 56, 57]. In contrast to the zero density case 
this local degeneracy is not removed by the NNN terms. 

In order to remove this degeneracy we consider quantum fluctuations in Chapter 5. We 
first treat the case of zero baryon number in Section 5.1, where we see that 1/N C corrections 
to the free energy removes the discrete degeneracy of the nearest-neighbor theory 3 . We 
find that fluctuations choose the same ground state chosen by the NNN terms. 

We treat fluctuations for the non-zero density system in Section 5.2. In Section 5.2.1 

3 This result was obtained by Smit using a Holstein-PrimakofT transformation on the quantum Hamil- 
tonian. 
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we first show that the local degeneracy is expressed by having zero modes in the 0(1) 
dispersion relations of excitations around the classical ground state. These have zero 
energy for all momenta. In Section 5.2.2, we show that these zero modes get non-zero self 
energy of order 1/N C , thus removing the classical local degeneracy. This is the phenomenon 
of "order from disorder" also present in the frustrated spin systems mentioned above. We 
also find that the energies produced are quadratic in momentum, which makes the zero 
modes type II Goldstone bosons. 

The analysis in Chapters 3-5 is mainly concerned with the nearest-neighbor theory that is 
symmetric under the too large group of U(ANf). In Chapter 6 we take into account next- 
nearest-neighbor interactions, and treat them as a perturbation on the nearest-neighbor 
ground state. In Section 6.1 we introduce the action of the sigma model which is now 
invariant only under U(Nf) x U(Nf). In Section 6.2, we show that some of the global 
degeneracy of the nearest-neighbor ground state is removed, and discuss how the vacuum 
expectation values of the sigma fields break chiral symmetry and discrete lattice rotations. 
In Section 6.3 we calculate the effects of these next-nearest-neighbor interactions on the 
spectrum of the nearest-neighbor theory. 

In Chapter 7 we summarize the thesis and point to directions of future research. 

All the work presented is this dissertation, except for Section 3.4, and some of the technical 
appendices, has already been published. The analysis of Chapter 2 has appeared already 
in earlier work [19, 20, 21, 22, 30] and in [58, 59]. The results of Section 3.3 (without the 
analysis itself) and Chapter 4 were published in [58] and [59]. Sections 5.1 and 5.2 were 
also published in [59] , whereas the remainder of Chapter 5 was published in [60] . Finally 
the treatment of NNN interactions of Chapter 6 was published in [61] and [62]. 



Chapter 2 

The effective Hamiltonian 



In this chapter we derive the effective Hamiltonian we use to investigate QCD at high 
densities. We begin by giving a short review on the Hamiltonian approach to lattice QCD 
in Section 2.1, where we stress important points relevant to this work, and fix some of the 
notations. Aiming to work at strong coupling, and aware of the fact that the continuum 
limit is at weak coupling, we take the point of view that the investigated Hamiltonian is 
an effective one that describes QCD at large distances. We elaborate on this formal point 
in Section 2.2. We then proceed to derive the effective Hamiltonian by strong coupling 
expansion in Section 2.3. We conclude by emphasizing the similarity of the Hamiltonian 
to the t-J model Hamiltonian of condensed matter physics. 



2.1 Hamiltonian approach to lattice QCD 

In this section we introduce the Hamiltonian formalism of lattice QCD. This formalism 
is more amenable to strong-coupling perturbation theory and to qualitative study of the 
ensuing effective theory [19, 20, 21, 22, 30] than the Euclidean formalism. 

The Hamiltonian formalism was first introduced by Kogut and Susskind in 1975 [63], 
shortly after Wilson's Euclidean formulation [64]. This formalism defines the theory 
of strong interactions on a three dimensional space lattice with lattice spacing a and 
continuous time t. A lattice site is denoted by a three dimensional vector n = (n x , n y , n z ) 
taking integer values, and a lattice link is denoted by (nfi), where \i = ±x, ±y, ±z. 

The quantum fields that describe quarks are the fermion fields ip^/ a that live on the sites 
n of the lattice. They have color indices a = 1, . . . , N c , Dirac indices a = 1, . . . , 4, and 
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flavor indices / = 1, . . . , Nf. It is more convenient to group the Dirac and flavor indices 
to a single index a that takes values from 1 to 47V/. The Fermi fields obey the following 
anti-commutation relations 

(2.1) 



For the gauge fields, that describe the gluons, one chooses to work in the temporal gauge 
that fixes A = 0, and removes one degree of freedom (and its conjugate momentum). 
As a result, one is left only with the link operators (U n fj) ab , and {U^^j = (C/ Tl+ ^ ; _^) ab . 
These represent gauge fields on the link (nft), whose two edges are colored with colors 
a and b. There are 6 links emanating from each site n. On each link reside the gauge 
fields (U n fi) ab and their conjugate momenta E % n ^ % — 1, . . . , N% — 1. The following are the 
commutation relations of this set of operators, 



n V)ab 



\ l u. 



nfj. 



ab 



E i - E j - 



= if jk E k nfl , 



for the set of the links (nft), and 



ab 



E % - E 3 - 



(2.2) 
(2.3) 

(2.4) 
(2.5) 



for the set (n,—fi). Here A 1 are matrices that represent the generators of SU(N C ) in 
the fundamental representation, and f tjk are the structure constants of SU(N C ). For 
N c = 2, Eqs. (2.2)-(2.5) are exactly the commutations relations that appear in the classical 
problem of the rotating top, where E l n ^ and E l n _^ generate rotations of space fixed and 
body fixed coordinate systems. One uses this analogy to write 



E i - - = - (lJ Adi -) E j - 



(2.6) 



where U^' is the rotation U n p, in the adjoint representation of SU(N C ). For a more 
detailed discussion see [63]. 

We now discuss the lattice Hilbert space. A state |Q) is the direct product on all lattice 
sites, 

|fi) = |fi) G (8)|n>*. (2.7) 

Here the first factor is the projection of the state \Q) to the gauge field sector, while the 
second factor describes the fermionic sector. 



2.1. Hamiltonian approach to lattice QCD 
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Any state is the following direct product 

|n>» = II (!">*), 



(2.8) 



We now concentrate on the Hilbert space of each site. The "lowest" state is the no- 
quantum "drained" state |Dr)„, defined as 



dDr>„ = 0. 

Applications of various ip^ a create the corresponding quarks on that site, 

\aa) n ^ = ^ a \Di) n . 



(2.9) 



(2.10) 



For the free theory, a — 1, ... , 2Nf correspond to creation of positive energy excitations, 
i.e. quarks, whereas a = 2Nf, . . . , ANf corresponds to creation of negative energy exci- 
tations, i.e. annihilation of anti-quarks. To use the usual quark-anti-quark language we 
write for each site, color, and flavor, 



( 6 T \ 

{ d] ) 



(2.11) 



create a quark and & an anti-quark. From Eq. (2.9) we see that b and S annihilate the 
drained state, which means that this state is empty of quarks, and filled with anti-quarks. 
The operator of local baryon number is 



N c N f , 

E E E [W " d t af d a s f ] n = w 1>ti> n - 2N f 



N, 



c =i S =T,4. f=i 



(2.12) 



According to Eq. (2.12), the baryon number of the drained state is —2Nf, corresponding 
to filling the site with anti-quarks. The vacuum |0) is the state with no quarks and no 
anti-quarks. This state is the filled Dirac sea on a single site and obeys, 



6|0) = d|0) = 0. 



(2.13) 



The baryon number of this state is B = 0. Because of Pauli exclusion principle we cannot 
put too many fermions on a single site. The maximum number of local baryon number 
will be 2Nf, and is found only in the state | filled) 



6 f | filled) =d\ filled) = 0. 



(2.14) 
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Below we will see that gauge invariance puts more restrictions on the single site Hilbert 
space in order that it be color neutral. 

Moving to the gauge Hilbert space, we also write it as a direct product of the form 

l«> G = II, (MgW (2-15) 

Cg>Ti/i 



One denotes the state with no electric field E by |0) G , 

£*|0) G = 0, Vi. (2.16) 

Any application of the link operators {U^^j on |0) G , creates states which correspond 
to flux lines on the link {njx). The state |0) G is the only state with no flux at all. Using 
the fact that the electric field operators generate a SU (N c ) algebra, one can distinguish 
between the different quantum states created by the link operators as follows. Define the 
quadratic Casimir operator 

E 2 nfl = £ (2-17) 

i=i 

It is clear that the fluxless state |0) G is an eigenstate of this Casimir, with zero eigenvalue. 
Next the commutations of the Casimir with the link operators are verified from Eqs. (2.2)- 
(2.5) to be 

[Ei„(uU) ab }-C F (uU) ab , (2.18) 
where Cp is the Casimir operator in the fundamental representation, and is 

N 2 - 1 

Cf = ^~ (2.19) 



2N, 



c 



for SU(N C ). 

This means that the state 



Unit) |0)g * s a ^ so an eigenstate of (2.17), with eigenvalue 



equal to Cp. One can now classify the states in |fi) G according to their E 2 eigenvalue. 
The result is a Hilbert space with a ladder-like structure. The lowest state is |0) G , with 
zero flux, and is a singlet of SU (N c ) . Repeated applications of the gauge field operators 
(ylip) b cre& te states with higher and higher values of flux and the operators E^fi measures 
the flux on the link (nfi). Indeed we shall see shortly that it is proportional to the (kinetic) 
energy of the gauge fields. 



2.1. Hamiltonian approach to lattice QCD 
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To complete the picture we now discuss gauge invariance. First recall that the start- 
ing point of this formalism was to choose the timelike gauge. This leaves only time- 
independent gauge transformations as a symmetry. The fermion operators belong to the 
fundamental representation of the gauge symmetry and transform as 



r n -(v n r b ^ b n , 



(2.20) 



with V G SU (N c ) given in general by 



V = exp 



(2.21) 



Note that here, and in the following discussion, we suppress the Dirac-flavor indices, 
irrelevant to gauge invariance. Using the anticommutation relations (2.1), one can show 
that the quantum operator that realizes Eq. (2.20) in Hilbert space as 



(2.22) 



is 



Vf = exp 



N?-l 



n 



E £ O^o 



(2.23) 



The gauge fields transform according to 

(u n(l ) ab ^{v n u nf yl +p ) ab . 



(2.24) 



Using the commutation relations (2.2)-(2.5), one shows that the quantum operator Vg 
that generates these rotations is given by 



Vg = exp 



N?-l 



E E £ \E E l 



n i=l 



(2.25) 



Finally putting Eq. (2.23), and Eq. (2.25) together and using Eq. (2.6), we see that 
operator that induces gauge transformations is 



V = exp 



N 2 c -l 

EE 

n i=l 



£ {Ph + Pf 



(2.26) 
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where p Fn , and p Gn are the color charge densities of the fermions, and gauge fields. These 
two quantities are given by 

P Fn = ^Krtl (2-27) 
Ph, n = E (<A - K-U^Ua) " ( 2 - 28 ) 

Note that Eq. (2.28) can be written as the covariant divergence of the electric field, 

p Gn = D ■ K ee V L • K + £ (l - U^) Ei_^, (2.29) 

where the lattice divergence is given by 

V L ■ K = ^ «, - Ei_^) . (2.30) 
n 

Since the lattice Hamiltonian is gauge invariant, we know that the generators of the gauge 
transformations commute with the Hamiltonian 

H,ph n + P? Fn ]=0, V* Vn, (2.31) 

This means that we can choose to work with a basis that block diagonalizes Pc + p l p- This 
breaks the Hilbert space to separate sectors classified by their p l G + p l F representations on 
the lattice. Each set of operators p'xt,™ describes a different physical case, with a different 
external distribution of color charge (that can correspond, for example, to infinitely heavy 
quarks etc.). To describe the physics of quarks with finite masses, and zero external gauge 
fields, we work with the choice p* xt = 0. Working in this subspace means that all physical 
states must be color singlets, since all gauge transformations are trivial. 

Finally we describe the lattice Hamiltonian H of SU (N c ) gauge theory with Nf flavors of 
fermions. H is given by 

H = H E + H G + H F . (2.32) 
Here H E is the electric term, a sum over links (np) of the SU (N c ) Casimir operator 

He = \ 9 2 Y,EI»- (2-33) 
Next is the magnetic term, a sum over plaquettes, 

H o = ^- 2 E( N c-^U p ). (2.34) 



2.2. The strong coupling limit 
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Finally we have the fermion Hamiltonian, 

Hf = ~i E E (n tW/.) + ^ ( 2 - 35 ) 

In Eq. (2.35) we suppress Dirac and color indices, and denote explicitly the flavor index 
/. The function D(j) is a kernel that defines the lattice fermion derivative. It should be 
sufficiently local, namely at least falling exponentially with distance. For example it can 
yield a naive, nearest-neighbor action if D(j) = a next-nearest-neighbor kernel with 
D(j > 3) = 0; a long-range SLAC derivative [30] with D(j) = —(—iy/j will however not 
be good enough, since it falls only algebraically. 

The choice of the fermion kernel is important, since it fixes the global symmetries of the 
Dirac Hamiltonian (2.35), and in turn the global symmetries of the effective Hamiltonian 
that we derive in Section 2.3. Superficially, the Hamiltonian has only the chiral symmetry 
of U(Nf)n x U(Nf)i regardless of the kernel used. This is however not true. For a 
naive kernel, one can spin diagonalize the fermions, and show that the Dirac matrices 
disappear. This means that the Dirac a = 1, ... ,4 indices become just a different kind 
of flavor, and the symmetry is then U(4Nf). Applying the same procedure to the next- 
nearest-neighbor kernel, one sees that the Dirac structure does not disappear, and the 
symmetry is indeed only chiral symmetry U (Nf) x U (Nf). We postpone further discussion 
in symmetries to the level of the effective Hamiltonian calculated in Section 2.3. 



2.2 The strong coupling limit 

Formulating a quantum field theory on a lattice is in fact a way of regulating the theory 
by providing it with a cutoff in coordinate space. As in any other regularization scheme, 
the cutoff must be removed at the end in favor of observable quantities. On the lattice, 
the removal of the cutoff is basically the procedure of taking the continuum limit of the 
theory, i.e., taking the lattice spacing a to zero. In order to extract continuum physics, one 
calculates a physical observable as a function of the bare parameters of the theory (e.g. a 
and the bare coupling go), and demands that the a — > continuum limit of the observable 
will be its measured value. This procedure tells us how the bare coupling go behaves as 
a function of the cutoff a. In the case of asymptotically free theories (and in particular 
in the case of QCD), the bare coupling g (a) approaches zero as a — > 0. This means that 
ideally we should examine the lattice theory at small coupling and small lattice spacings. 
The regime where these both are small enough is determined by the scaling behavior of 
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physical observables one calculates. This approach has become standard, and is realized 
by doing Monte Carlo simulations of the lattice theory. These, performed close to the 
continuum limit, give us information on non-perturbative aspects of QCD. 

As mentioned in Chapter 1, these Monte Carlo approaches fail when the theory is defined 
at finite chemical potential. In order to examine the dense system, one is forced to resort 
to other methods of calculation. We turn to the strong coupling regime of QCD, which 
is far from the continuum limit. As a result one cannot trust numerical figures that come 
out of strong-coupling calculations (especially since we use only the lowest order and do 
not perform any extrapolation to small coupling). Despite this, examining the strong- 
coupling regime gives insight on qualitative features of the dense nuclear system. The 
physical picture we have in mind is that if one would perform a block-spin transformation 
of the lattice theory near its continuum limit, and integrate out high momentum degrees 
of freedom, then one will be left with a lattice theory for large distances (of the order of 
hadronic scale distances of 1 fm.) The resulting effective theory will be strongly-coupled, 
and have the symmetries of QCD. In particular it will have an £77(3) gauge group (or 
in general an SU(N C ), where N c is the number of colors). It will also have the global 
symmetry 

SU(N f ) L x SU(N f ) R x U(1) B . (2.36) 

Its dynamical physical features should include confinement and spontaneous chiral sym- 
metry breaking (at zero density). All these important features are present in lattice QCD 
with strong bare coupling, and a suitable definition of the fermion kernel, which we take 
to be our starting point. In this work we show that the Hamiltonian formulation of this 
theory is tractable even with non-zero baryon density. 

The only property that the theory explored here does not have is the anomaly of the 
axial U(l) symmetry. In QCD, this symmetry is classically conserved, but broken in 
the presence of gauge fields. This is expressed in the fact that the corresponding singlet 
axial current has a non-zero divergence. This fact is most easily seen in perturbation 
theory in continuum QCD [65], and is verified experimentally by the large mass of the 
rj' meson, which should have been as small as the pion masses in the absence of this 
anomaly. The same calculation can be done with weak-coupling perturbative methods of 
lattice QCD [66]. There, the lattice doubling of the fermions plays a crucial role and leads 
to conservation of the singlet axial current; the doubled fermions fall into two groups with 
opposite chiral charges, and the contribution of each group to the anomaly is canceled. In 
order to obtain the anomaly in a lattice theory, one must break chiral symmetry explicitly 
(for example by a Wilson term) to remove the doubling. The continuum limit is then taken 
carefully to restore all chiral symmetries except the axial U(l). 



2.3. Strong coupling perturbation theory 17 



We do not use Wilson fermions, and as mentioned also do not work in the continuum 
limit. Our fermion kernel leads to full species doubling in the weak-coupled continuum 
limit, but this fact by itself is not a real problem, since we work at strong coupling, where 
we do not have quark excitations. Our Hamiltonian will not lead to doubling of mesonic 
excitations, as long as it has the correct symmetry. Nevertheless, we can not avoid the 
absence of the anomaly; U(1)a is not anomalous, and is conserved at the quantum level 
as well. At this point we stress the fact that we regard our theory as an effective theory 
for QCD at large distances. In the process of the block-spin transformation that begins 
from a theory with no axial singlet symmetry, a U(1)a violating term will have to be 
generated. Therefore in order to have physical features related to this term, we will have 
to add it by hand, without knowledge of its exact strength. Note that in addition, we 
know neither the exact value of the coupling, nor the exact way the fermion kernel behave 
(we do assume that it falls fast enough). 

In practice we can add the U(1)a violating term in any stage of the calculation. At the 
level of the effective lattice Hamiltonian introduced below, it will be a 't Hooft vertex. For 
two flavors it can also be expressed in terms of the color singlet objects of Section 2.3.1, 
and for any number of flavors it is easiest to add this term at the level of the non-linear 
sigma model derived in Chapter 4 [32]. This term will have some simple effects on the 
low energy spectrum we derive, and we discuss them in Chapter 6. 



2.3 Strong coupling perturbation theory 

Our starting point is therefore the Hamiltonian (2.32) with large bare coupling. For g > 1 
the ground state of H is determined by Hp alone to be any state with zero gauge flux, 
whatever its fermion content, 

(2.37) 

G 

These states have energy eo = and are degenerate with respect to all the fermionic 
degrees of freedom. We consider perturbation theory in V = H G + Hp. Both Hq and Hp 
are sums of operators that are strictly bounded, independent of g except for the explicit 
coefficient in Hq. We can dismiss first-order perturbations by noting that H G and Hp 
are multilinear in link operators U and W , which are raising/lowering operators for the 
electric field; thus there are no non-zero matrix elements within the zero-field sector. 



|0> 



n = °) 



nil 
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We proceed to higher orders, and seek an effective Hamiltonian that acts in the zero-field 
sector [67]. Define Pq to be the projector onto the subspace of all the E = states. Then 
perturbation theory in V gives an effective Hamiltonian, 

#eff = PoVQDVPo + PqVQDVQDVPo + ■■■. (2.38) 

Here Q = 1 — Pq projects onto the subspace orthogonal to the E = states; the operator 
D = (e — He) 1 supplies energy denominators, so that 

QD= Yl |A>— ^— <A|. (2.39) 

E^O states C ° Ca 

The intermediate states |A) contain flux excitations. In second and third order the patterns 
of flux can only be strings of length j in the fundamental representation of the color group. 
Thus the energy denominators are 

eo-ex = -\g 2 C F \j\. (2.40) 

The perturbations Hq and Hp are explicitly of 0(1/ g 2 ) and 0(1), respectively; each 
energy denominator gives a factor of l/g 2 . Thus to 0(1/ g 4 ) we can forget about Hq- Our 
result to this order is 

H eS = P H F DH F P + P H F DH F DH F P . (2.41) 

Since H F has no non-zero matrix elements within the E = sector, we have dispensed 
with Q in Eq. (2.41). The first term in Eq. (2.41) arises for any value of N c and is 0(1/ g 2 ); 
the case N c = 2 must be treated carefully, but all cases N c > 2 are generic. The second 
term is special to N c = 3 and is 0(l/g 4 ). 

2.3.1 Second order: the antiferromagnet 

We calculate explicitly the first term in H c g. Each term in H F creates a string of flux of 
length j, which must be destroyed by the conjugate term. Thus 

H$ = 2J2[-K(j)} £ (^a^V .J (0 | (jj u) (JJ ^) | ) (^a^*) , (2.42) 

j>0 n^i 
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where we define 

KU) = 0T\ > (2 ' 43) 

The matrix element of the gauge fields yields ^-<5 a cAc, independent of j. 

As they appear in Eq. (2.42), each ip* is next to a ip on a different site. This invites a 
Fierz transformation on the product of fermion fields, which we write generally as 

(Vjo^,) (4^i) = <W&i - 7 E s a (Vl^) (l^Vi) • (2-44) 

Here k, I are combined site, flavor, and color indices, and we have assumed that k and 
/ are always different while j and k might be equal [as in Eq. (2.42)]. The matrices T A 
are the 16 Dirac matrices, normalized to (T A ) 2 = 1, and we have defined 

4 = ^Trr^r A a M = ±l. (2.45) 

This sign factor is ±1 according to whether T A commutes or anticommutes with a M ; it 

(2) 

will be a constant companion in our calculations. As they appear in H^ s , the indices 
i, I are the same site and color but different flavors, and likewise j, k. Leaving the flavor 
indices explicit, we obtain 

H< $ = 7^ E * OX (^ f/ rV) (^*rV) + -- d E*0')E(^V) . (2.46) 

Each fermion bilinear in parentheses is a color singlet located at a given site. The second 
term contains the baryon density 1 B' n = N~ 1 ipl l ip n , and the sum J2 n B' n is the total 
baryon number B'. A final important note is that one can easily verify that the effective 
Hamiltonian H cS commutes with the local baryon number B n (and B' n ) . This is a result 
of the fact that each of the local operators in Eq. (2.46) does not change the fermion 
occupation of a site. 

We now combine the Dirac indices with the flavor indices and write 

(^t'Fy 9 ) (if>ter A if> f ) = 8 (Vm^)^ (VaT^)^ . (2.47) 

1 This baryon number is positive semidefinite, and is zero for the drained state. The conventional 
baryon number B n is zero in the half-filled state, and thus B' n = B n + 2Nf. 
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We have defined new matrices M v as direct products of the 4x4 Dirac matrices and the 
U (Nf) flavor generators, 

M v = T A <g> A a , (2.48) 
and we have normalized them conventionally according to 

TrMW' = i<5 w . (2.49) 



The M v generate a U(N) algebra, with N = 4JV). 
An alternating flip 



ij n (2.50) 



- n 

(spin diagonalization [20]) removes the matrices from the odd-j terms in Hp, and 
hence removes the sign factors s^ A from the odd-j terms in H^. We have finally 

H « = W C ^ KU) (<)- y J ^ MV ^^ M ^) n+jPl - (dNc^m^J B'. (2.51) 

The odd-j terms are of the form M ■ M which can be written in any basis for the U (4Aj) 
algebra. The even-j terms, however, contain which is defined only in the original basis 
(2.48). 

We now turn to describe the zero-field sector Hilbert space in which the Hamiltonian (2.51) 
works. In this sector in which we work, Gauss' Law constrains the fermion state at each 
site to be a color singlet. The drained state |Dr), with ip%/ a \Dr) = for all (a fa), is the 
unique state with B' — 0. The other color singlet states may be generated by repeated 
application of the baryon creation operator, 

blp... = e ab ...^ n aa ^ ■ ■ ■ , (2.52) 

with N c operators ip. (Here and henceforth, the indices a, (3, . . . combine the flavor and 
Dirac indices.) 

At each site n the operators 

Ql = ^ n M^ n = tfrKfiti? (2-53) 

generate a U(N) algebra, with = 4Aj. The drained state is obviously a singlet under 
this algebra. The creation operator 6^... is symmetric to any exchange of its A c Dirac- 
flavor indices. As a result it is in the symmetric representation of U (N) with one row and 
A c columns (see Fig. 2.1). 



2.3. Strong coupling perturbation theory 
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Figure 2.1: Young diagram of the representation of U(ANf) carried by the baryon operator 

Repeated application of v a p... to the drained state gives the state 

6 t 6 t ---|dr). (2.54) 

The form (2.54) places the single site state within an irreducible representation 
of SU(ANj). Let us elaborate on this point. The single site state (2.54) is generated by 
an operator that has two sets of n indices; SU (N c ) color indices, and U(ANj) Dirac-flavor 
indices. The number n is the number of fermions occupying the site. Gauge invariance 
restricts the color indices to be in an SU(N C ) singlet, which means that 

n = m-N c , (2.55) 

where m is an integer. In fact, m is the number of applications in Eq. (2.54). More 
precisely, gauge invariance restricts the color indices to belong to the singlet representation 
s given in Fig 2.2. 
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Figure 2.2: The representation of SU(N C ) carried by the color indices of the singlet site 
state, that has mN c fermions. The baryon number is B = m — 2Nf. 



Next, we recall that states in an irreducible representation of any special unitary group 
SU(k), are also an irreducible representation of the permutation group. In our case, 
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this means that the n color indices belong also to the s irreducible representation of S n . 
Our aim is to determine the irreducible representation r of the Dirac-flavor indices, with 
regard to the SU(4Nf) group. Since this is also the S n irreducible representation of the 
n Dirac-flavor indices, we find that the state (2.54) can belong to any possible part of the 
reducible representation r <g> s of S n . Finally, we consider the fermionic nature of the state, 
and realize that the state belongs only to the completely antisymmetric representation A 
of S n included in f <g) s and given in Fig. (2.3). 



mN 

c 



I I 

Figure 2.3: The representation of S n carried by the complete set of indices of the singlet 
site state. 



For given s and r, A occurs in s (g> r if and only if r is the S n conjugate representation of s 
(which means flipping the Young tableau's rows and columns). Moreover the representa- 
tion A occurs only once in the decomposition of s <8> conjugate (s). The conclusion is that 
the S n (and therefore the SU(4Nf)) representation of the Dirac-flavor indices is given in 
Fig. 2.4. 

We conclude by emphasizing the importance of the result of this section. Since the second- 
order effective Hamiltonian H^ s preserves B n , the baryon number on each site, then any 

(2,) 

distribution of B' n defines a sector within which H c q is to be diagonalized. In other words, 
baryons constitute a fixed background in which to study "mesonic" dynamics. The baryon 
number at each site fixes the representation of U(N) at that site, which is the space of 
states in which the charges Q n n act. 
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Figure 2.4: The representation of U(ANf) carried by the spin in the effective antiferro- 
magnet. m is related to the baryon number at the site according to m = B + 2Nf, with 
\B\ < 2N f . 

We now describe the global symmetries of the Hamiltonian. The j — 1 terms in Eq. (2.51) 
are of the form QnQn+p,' anc ^ ^hey commute with the generators 

Q v = J2Ql (2-56) 

n 

of a global U(N) symmetry group. This symmetry is in fact familiar from the lattice 
Hamiltonian of naive, nearest-neighbor fermions: Spin diagonalization of Nf naive Dirac 
fermions transforms the Hamiltonian into that of ANf staggered fermions. In the weak 
coupling limit, there are in fact 8Nf fermion flavors — the doubling problem. This doubling 
is partially reflected in the accidental U (4iVj) symmetry, which is intact in the g — > oo limit 
and is respected by the effective Hamiltonian. Retaining terms in the fermion Hamiltonian 
(and thus in H^g ) that involve odd separations j does not break this symmetry. 

The Nielsen-Ninomiya theorem [23] guarantees that any fermion Hamiltonian of finite 
range will possess the full doubling problem. This is a statement, however, about weak 
coupling only, since the dispersion relation of free fermions is irrelevant if the coupling 
is strong and the fermions are confined. It is interesting that the accidental U(4Nf) 
symmetry nonetheless survives into strong coupling as a vestige of doubling. 

The terms in Eq. (2.51) with even j, on the other hand, break the U(N) symmetry, as do 
even-j terms in the original fermion Hamiltonian. It is easy to see via spin diagonalization, 
which leaves the even-j terms unchanged, that the only generators left unbroken are the 
Q v corresponding to 

M v = 1 <g> A a and 7s <8> A a , (2.57) 

which form the U(Nf)i x U(Nf)n chiral algebra. This of course makes no difference 
to the Nielsen-Ninomiya theorem, which will enforce 8-fold doubling in the perturbative 
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propagator even without the U(ANf) symmetry. If we are interested in the realization of 
the global symmetries of the continuum theory, though, we can study this lattice theory 
which has the same symmetry. The simplest theory one may study is thus one containing 
nearest-neighbor and next-nearest-neighbor terms. We shall proceed to discard terms 
with longer range; we shall begin with the nearest-neighbor theory, with its accidental 
doubling symmetry, and later break this symmetry to U(Nf) L x U(Nf) R with the next- 
nearest-neighbor terms. 

Two essential differences will always remain between this lattice theory and the continuum 
theory. One is the presence of the axial U(l) symmetry on the lattice. This symmetry is 
exact, broken by no anomaly, and may make the drawing of conclusions for the continuum 
theory less than straightforward unless it is broken by hand. The other difference is the 
fact that the effective Hamiltonian for baryons (see below) is also a short-ranged hopping 
Hamiltonian. If the baryons were almost free, we would say that they are surely doubled 
like the original quarks. The fact that the simplicity of the hopping terms is only apparent, 
and that the baryons are still coupled strongly to mesonic excitations, offers the possibility 
that doubling may not return. 

2.3.2 Third order: the baryon kinetic term 

Finally we show here the form of the third-order contribution to H e g, which only exists 
in the case of N c = 3. This term is calculated via 

H$ = P V F DV F DV F P . (2.58) 

For a single link, we have 

(0\U ab U a , b ,UjO) = l -e aa ,,e whl (2.59) 
and the same can be proved for a chain of links, 

iO\(Uu)jYlu) c , b ,(Ilu)jO)= l -e M e wi ,. (2.60) 
Thus [f,g,h,... are here flavor indices], 

h$ = -i e m £ (^v^-a) K^^a) (#«^-a) + h - c - 

(2.61) 
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The kernel is 

K(j) = [J,) (2.62) 

s(y 2 c F \j\) 

Again, spin diagonalization simplifies the odd-j terms, but not the even-j terms. The 
result is 

H$ = H(% + //(3) n , (2.63) 

with 

^% = -i E E + ^ ( 2 - 64 ) 

j>0 Tl^, 

j odd 

where rj^n) is the usual staggered- fermion sign factor, and 

#c ( l = -i E E # K ® «m ® «,J/r CaO) + (2.65) 

j>0 nfj, 

jeven 

where (^(n) = (— l)^ 1 ^™". The baryon operators are 

b 1 = e afec ^VV 7 , (2.66) 

where we have written / to represent the compound index {afij}, taking values in the 
symmetric three-index representation of U(N) in Fig. 2.1. More precisely, 



a, /? and 7 are indices of the fundamental representation of U(ANf), and the numbers are 
color indices. Note that these baryonic operators do not obey the canonical anticommu- 
tation relations of fundamental fermions [68, 69, 70, 71, 72, 73], but rather obey 

{bi, = A IP 6 nn , ^ 5 IP 5 nn , (2.67) 

The operator A n > has a complex index structure and is given in Appendix A. 

The odd-j part of H^, like that of , is symmetric under the U (N) doubling symmetry. 
The even-j part breaks U(N) to U(N f ) L x U(N f ) R . is a baryon hopping term. As 
mentioned in the introduction, however, its simplicity is deceptive. 
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The separation of if eff into a canonical kinetic energy and an interaction term is a chal- 
lenge for the future. 

The baryon operators responsible for the hopping are composite operators that do not 
obey canonical anti-commutation relations. If this were not the case, then the effective 
Hamiltonian in third order would strongly resemble that of the t—J model of condensed 
matter physics [50], 

Ht-j = -t^clcjs + (Si ■ Sj - ^f-) + J'. (2.68) 

W> (ij) v 4 J 

s 

Here Cj S is an annihilation operator for an electron at site j with spin s, and the number 
operators rij = cfcj and spin operators Sj = |c|crcj are constructed from it. The added 
term J 1 is a more complicated hopping and interaction term. The t-J model describes a 
doped antiferromagnet; it arises as the strong-binding limit of Hubbard model, a popular 
model for itinerant magnetism and possibly for high-T c superconductivity. The model is 
not particularly tractable and, absent new theoretical developments, does not offer much 
hope for progress in our finite-density problem. It is nonetheless worth pondering the fact 
that a model connected, however tentatively, with superconductivity appears in a study 
of high-density nuclear matter. 

A first step will be to replace these baryon operators with simpler operators that mimic 
only some of their features. Reintroducing all the properties of the real baryons can be 
done using methods of Hilbert space mapping (see for example [69]). These methods have 
been introduced in the context of atomic and nuclear physics to treat quantum fields that 
are composite operators and therefore do not obey canonical commutation relations. One 
adds to the physical, original Hilbert space an auxiliary Hilbert space that is occupied by 
canonical operators, that replace the composite (non-canonical) operators. The price for 
this replacement is mainly very complicated interactions between the composite operators 
(baryons in our case), that has up to six-dimension operators. One can simplify these 
interactions by solving a non-Hermitian effective Hamiltonian, whose spectrum can pro- 
vide us with insight on the real spectrum only if its projection on the physical spectrum 
is clustered at low energies. 

2.4 Summary 

To conclude we gather the results of the previous sections and present the effective Hamil- 
tonian derived from the strong coupling limit of lattice QCD. We restrict to the second 
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order contribution to the effective Hamiltonian. At this order the local baryon number 
operator commutes with the Hamiltonian. As a result, the baryon number on each of the 
sites is fixed, and the Hamiltonian acting on the degenerate subspace of fluxless states 
takes a block diagonal form, in which each block corresponds to a certain distribution of 
baryon number on the lattice. This also means that the last term in Eq. (2.46) is just a 
constant and can be dropped. One must now diagonalize within each of the blocks, 
i.e., for a certain distribution of baryon number, and therefore a certain distribution of 
representations of the Q n operators. 

The effective Hamiltonian has the structure of an antiferromagnet with interactions that 
connect sites over j lattice links, with j = 1, 2, 3, . . .. It has the following form (we drop 
a constant proportional to B'). 

H= E JjQlQl +jfl - (2-69) 

3=1,2,3,... 

Here Q\ are U(N) charges at site n, with r] — 1, ... , N 2 . This Hamiltonian moves mesonic 
excitations around the lattice, leaving the baryon density fixed. The states at n comprise 
a U(N) representation whose Young tableau has N c columns. The number of rows m 
depends on the baryon number B at n according to 

m = B n + 2N f . (2.70) 

The sign factors are given by 

= 2Tr M v a fl M v a fl . (2.71) 

Here M v = T A <g> A\ where T A are the 16 Dirac matrices, and A* are the Nj flavor 
generators. This makes M 71 ; r] — 1, ... , 16iV|, the generators of U(N) in the fundamental 
representation. We fix their normalization by, TrM^M 71 ' = ^S m >. The matrices are 
the 4x4 Dirac matrices times the unit matrix in flavor space. 

In terms of the fermion kernel D(j) the effective couplings Jj are 

j W)P > ( 2 72) 

Recalling that we regard the strong coupling Hamiltonian as an effective Hamiltonian that 
describes QCD at large distances, then a physical fermion kernel must fall with j at least 
exponentially. In view of the symmetry structure of the even and odd terms in H e s it 
is sufficient to take into account only nearest-neighbor and next-nearest-neighbor terms, 
and fix D(j > 2) = 0. The real values of the exchange coupling should be determined by 
the block-spin transformation discussed in Section 2.2. 



Chapter 3 



Mean field theory with Schwinger 
bosons 



In this chapter we follow Arovas and Auerbach [45] and express the nearest-neighbor 
Hamiltonian derived in Chapter 2 in terms of Schwinger bosons 1 . We write the partition 
function as a path integral whose action is quartic in the bosons. Using a Hubbard- 
Stratonovich (HS) transformation, the interaction is decoupled, and the action becomes 
a quadratic form. After integrating the bosons we get an effective action for the HS 
fields. This action is proportional to the number N = ANf of boson flavors, and can be 
approached in the large N limit. (To be precise the number N c of colors is also large 
in this limit, and N c /N is kept fixed.) We restrict to uniform and time-independent 
configurations of the HS fields, and replace the action by the mean field (MF) action Smf- 
The latter depends on the MF values of the HS fields, and its minimization yields the 
ground state, which we analyze. 

We first extend the analysis of [45] to three dimensions and show that at zero density and 
for large enough values of N c /N the ground state spontaneously breaks the global U(N) 
symmetry. The ordered phase is the Neel phase of the antiferromagnet and in terms of 
chiral symmetry respects only vector rotations. 2 

Next we examine how the MF approach changes when we insert a non-zero concentration 
c of baryons into the lattice. The baryons cannot move, and are represented by sites 

1 This paper is concerned with SU(N) antiferromagnets, in order to understand the SU{2) antiferro- 
magnet. For more work in that context see [43, 44, 46, 47, 48, 49, 50]. 

2 This is a subtle point, since the initial global symmetry U(N) is much larger than chiral symmetry. 
We discuss this point in chapters 4 and 6. 
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filled completely with baryon number. On such a site there is only one state, which is 
an SU(N) singlet. This means that the generalized "spin" operators defined in Eq. (3.3) 
obey Q v = 0, and that there are missing links in the nearest-neighbor antiferromagnetic 
Hamiltonian (2.69). Therefore the system becomes similar to an antiferromagnet with 
magnetic vacancies. We calculate the MF action to this system by two methods. For 
concentrations c = | and | we exactly minimize the MF action and see that the Neel 
phase shrinks as c increases. For other values of c 1 we disorder the baryon background 
and use a quenched approximation [74]. 



3.1 Schwinger bosons 

Our starting point is the effective nearest-neighbor Hamiltonian derived in Chapter 2, 

with Ji/N = J\ from Chapter 2. For brevity we drop the tilde sign in the rest of the 
chapter. The operators Q 11 generate the U (N) algebra 

[Qni Qm\ — iP^Qn^nm, (3.2) 

and are defined as 

Ql = E tiFKptft- (3-3) 

o/3,a 

Here a and f3 are Dirac-flavor indices that range from 1 to N, while a is a color index that 
ranges from 1 to N c . The matrices M v are the generators of U(N) in the fundamental 
representation. The Hilbert space on site n, occupied by B baryons, is an irreducible 
representation of U (N) that corresponds to a rectangular Young tableau of the form of 
Fig. 2.4. It has iV c columns, and its number of rows m obeys m = B + N/2. Following 
[45, 46] we first restrict to configurations with zero average baryon number. In particular 
we examine configurations with baryon number B on the even sites and — B on the odd 
sites. This means we put conjugate U(N) representations on adjacent sites, represented 
by Young tableaux with m rows on the even sites, and with (JV — m) on the odd sites. 
The case of exactly zero baryon number on each site is the case of m — N/2. 

We now show how to construct the Hilbert spaces using bosons. On an even site n, we 
define a boson operator b aPtn . a is again a Dirac-flavor index that ranges from 1 to N, 
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but the index p is now an auxiliary "color" index that ranges from 1 to m. The bosons 
belong to the fundamental representations of both SU(N) and SU(m), and obey 

(3.4) 

The Hilbert space on this site is constructed by mN c applications of the operators b^ n on 
the no-quantum state |0) 

^•••;pg-)~Okn-|o>- (3.5) 

We denote the SU(N) representation of the Dirac-flavor indices (a, (3, . . .) as r, and the 
SU{m) representation of the color indices (p, q, . . .) as r'. This means that the state (3.5) 
belongs to a representation r" e r x r' of the permutation group S m N c - In particular, 
since the state \a(3 ■ ■ ■ ;pq---) is created with bosons, r" is the completely symmetric 
representation of S m N c . This means that r' must obey r = r'. As a result we find that 
in order to constraint the SU(N) indices of Eq. (3.5) to belong to Fig. 2.4, it suffices 
to restrict its SU(m) indices to a representation with the same tableau. This is easy. 
Because Fig. 2.4 has m rows, it is a singlet of SU(m), so any state \Q) must obey 

N 



On the odd sites we choose the bosons to be in the conjugate-to-fundamental represen- 
tations of SU(N) and SU(m), and consider them as hole operators. The commutation 
relations (3.4) are the same, and so is the restriction (3.6), which is now interpreted as 
counting the number of holes created by applying the hole operators to the filled singlet 
state with N rows in its Young tableaux. 

In term of the bosons, the operators are given by 

n ~ I bl ■ (-API)* -b n ne odd. [6J) 

Here we have suppressed all indices. The difference between the even and odd defini- 
tions of Q\ above is inserted in order that on both sublattices the operators Q\ belong 
to the adjoint representation of SU(N); the fact that b transforms in the fundamental 
representation on the even sites and in the conjugate-to-fundamental representation on 
the odd sites is compensated by the fact that while M v generates the fundamental, —M v * 
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generates its conjugate. Using these definitions, the Hamiltonian becomes the following 
normal ordered quartic form, 

H = -^£^[.^.4,*], (3.8) 

N 

(•^■nfijpq = ^ap,n,^ag,Ti+/i- (3-9) 

«=1 

Here we have used the hermiticity of M v and the relation 

EK^aV^VV' (3-io) 

We now introduce bosonic coherent states (see for example [50]). This enables us to write 
the partition function as a path integral. 

For each bosonic operator b (we suppress all site, Dirac-flavor and color indices), one 
defines the bosonic coherent state \z) as 

|^) = e 2bt |0). (3.11) 

Here |0) is the state annihilated by b, and z is a complex number. This state obeys 

(3.12) 
(3.13) 

^ 2 \z)(z\, (3.14) 
-\ z \\z\d\z). (3.15) 

Here O is an operator that depends on b and . We use Eq. (3.15) to write the partition 
function 

Z = Tr'e- f,H , (3.16) 

where (5 = 1/T, as 

Z = J' d 2 z (zo\e- pH \zo) e" |2 ° |2 . (3.17) 

Here the trace in Eq. (3.16) and integral of Eq. (3.17) are constrained to contain bosonic 
states that obey Eq. (3.6) only. 



b\z) 


= z\z), 


( Z l\ Z 2) 




1 




Tr (6) 


= ! <e 
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Next, we regard the inverse temperature (3 as Euclidean time and slice it N T times, 

(zo\e- fiH \zo) = (z \e- eH e- eH e- eH ---e- £H \z ). (3.18) 

Here e = (3/N T and is small enough to write 

e - eH = (l-eH). (3.19) 

Inserting the unity (3.14) in between each of the 1 — eH factors in Eq. (3.18) we have 

n t -i . 

Z= [] / d 2 z T e- z * ZT {z T \l -eH\z T+1 ), (3.20) 

r=0 J 

where z^ T = z . Each of the unities inserted is denoted by a different Euclidean time 
index r = 0, . . . , N T — 1, hence the notations. Next we evaluate the matrix elements of 
Eq. (3.19) to 0(e) using Eqs. (3.12)-(3.13). The result is 

(Z T \1 - eH\z T+1 ) = e z>T + i-eH(z* T ,z T+1 ) _ ( 3 21 j 

Here H(z*, z T+ i) is given by replacing all with z* and all b with z in the quantum 
Hamiltonian. This procedure gives in the following path integral. 



Z 



f DzDz*e- s l[5lj2z^ n z aq>n -N c 5 pq ), (3.22) 

Jz —z Nt npq y a j 

N T -1 

S = e -z*d T z T + H(z*,z T ), (3.23) 

T=l 



h = —E^KW^W]. ( 3 - 24 ) 

N 

(•^■n^( T )) pq — ^ ^ap,Ti^Q(j,Ti+/i- (3.25) 
a=l 

Here one assumes smooth Euclidean time configurations of the fields, and we have sup- 
pressed site, Dirac-flavor and color indices in the first term of Eq. (3.23). A note is needed 
regarding the definition of the temporal derivative of Eq. (3.23). To be unambiguous we 
write 

ed T f(r) = f T+1 -f T . (3.26) 

In appendix B we show that ignoring this slicing procedure, and taking the continuum 
limit of the derivative, leads to slightly different results for energies. The difference is 
usually ignored, but in our context it remains important. 

Finally let us note the 5 function in Eq. (3.22). This constrains the boson fields and keeps 
the Hilbert space explored by the path integration within the quantum constraint (3.6). 
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3.2 Hubbard-Stratonovich transformation 



We proceed to implement a Hubbard-Stratonovich (HS) transformation to decouple the 
quartic interaction Eq. (3.23) into quadratic forms. We use the identity 



// p(3 ON \ 

DQDQ ] exp f -Tr ^ dr (Q n , - ( Ji/27V)^)t _ (Q^ _ (J^N)^)) = 1, 

(3.27) 

where we suppress all auxiliary color indices, and symbolically write the temporal sums 
as an integration over the variable r. The HS fields are complex matrices in the color 
indices p and q, and live on links of the lattice. Eq. (3.27) leads to 



expQ^Tr J <iT^(T)^(r)) = J DQDQ^ exp 



Tr 



Q 



(3.28) 



ap,n Z aq,n N c 5 pq 



(3.29) 



We also write the constraint in Eq. (3.22) as 

L pg \ a 

where A^(r) is a real valued matrix field. Using Eq. (3.28) and Eq. (3.29) the partition 
function becomes 

Z = J DzDz*DQDQ^DXe~ s , (3.30) 
r ( 2N 

S = drTr ]T— Ql^Q nil + iY,N c \ n -Y, z l\ d r + ^]z n 

(, n^i "1 n n 



[^Ti/j^Ti/i + ^linQn/j, 



(3.31) 



The action is now quadratic in the bosons, and they can be integrated. This gives the 
following effective action for the Q and A. 

S = N [ drTr yQ^Q^ + i£ + logG" 1 } . (3.32) 

Here G~ l is the propagator of a single boson in the presence of Q and A, while the 
general propagator is proportional to unity in the Dirac-flavor indices. Therefore the 
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boson integration gives an overall factor of N, and logG" 1 is of 0(1). Also note that we 
have introduced k = N c /N. 3 This will be useful in the following sections when we take 
the large N limit and keep k fixed. 

Let us make a final remark on the HS transformation. This exact transformation is based 
on Eq. (3.27), which is always true for Ji > 0. As a result we see that we can perform 
a HS transformation only if the interaction term in Eq. (3.23) can be written in the 
form (3.8). The latter was received with the aid of Eq. (3.7), which is true only for the 
special baryon number configuration of conjugate representations on adjacent sites (that 
has zero total baryon number). There are two additional simple cases in which one can 
write the magnetic interaction in the form of Eq. (3.8) [45]. The first is the ferromagnet 
with identical representations on all sites, and the second is the antiferromagnet, written 
with fermions and zero baryon number on all the sites. The HS transformation for the 
fermions was examined in the large N limit by many authors (for example [46]), and 
seems to lead to a disordered phase. This limit is of no interest to us, since we expect 
that at zero density, our Hamiltonian describes an ordered phase with spontaneously 
broken chiral symmetry. 

Antiferromagnetic interactions for other baryon distributions can also be written in the 
form of Eq. (3.8). The key point is that for each link one must have a representation of the 
generators analogous to Eq. (3.7). This however makes calculations more difficult. For 
example let us take uniform baryon density, with the same baryon number B on all sites. 
This distribution is realized by having the same U(N) representations on all sites, and 
therefore the same number of rows to in the Young tableaux of all sites. We now build the 
Hilbert space on even and odd sites differently, making the boson operators on the even 
and odd sites belong to the fundamental and conjugate-to-fundemental representation of 
SU (to) and SU (N — m) respectively. This choice of bosons indeed leads to a relation 
similar to Eq. (3.7). This means that while the bosons on the even sites have an auxiliary 
color index p that takes values in [1,to], the auxiliary color index q on the odd sites will 
take values in [1, N — to]. The bosons on the two sublattices are of a different kind and 
calculations become more difficult. For example, the HS field will now be rectangular 
matrices with to rows and N — to columns. 



3 k is equal to r/4, where r is the ratio N c /Nf mentioned in Chapter 1, and in the abstract and 
summary of the thesis. 
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3.3 Mean field theory at zero density 

Taking the large N limit with k fixed allows a stationary approximation for the path 
integral. We therefore proceed to work with the following ansatz, 

QW = Q 5P9 ' ( 3 - 33 ) 

A«(r) = i\S", (3.34) 

where we take Q and A to be real. Using this ansatz we find that the MF action Smf 
defined by 

S MF = {NmflNs)' 1 log Z MF , (3.35) 

becomes 

Smf(Q, \) = ^-k\ + -±-Tr logG-^Q, A). (3.36) 

To evaluate the last term, we perform the integration over the bosons with the constraint 
and HS fields replaced by their MF values. (Here we have a single flavor and "color" 
field.) 



det G(Q, A) = J DzDz ] exp J J dr ^ 



d 

Z n (&T — A) Z n + Q ^ ( z nZn+p, + %n Z n+n 



(3.37) 

We redefine the bosons on the even and odd sites as z^, where N is an index denoting 
sites on the even sublattice, which is an fee. We now have two degrees of freedom on each 
site of the fee. The first is the original boson z e on N, and we use the convention that 
the second is the boson residing on the odd site N + z/2, one link above the even site 
occupied by z e N . A We now apply the Fourier transform, 



*M = ^g*V.e^— >x|l_, s , /2 HI , (3.38) 

where u n = 2im/T are Matsubara frequencies, and k belongs to the first Brillouin zone of 
an fee (i.e. to the unit cell of a bec lattice.) In momentum space the inverse propagator 
is 



We use the convention that the lattice spacing of the fee lattice is 1. 
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with the 2x2 internal space defined by 

/ 

fk,n = 



(3.40) 



V Z -°k,-n ) 



and 7fc = | X^=i cos(A;^/2). The spectrum of this mean field theory is given by the poles 
of the propagator, 

co 2 n + A 2 - ±d 2 Q 2 ll = 0. (3.41) 

Near = 0, the energy described in Eq. (3.41) is E\ = A 2 + c 2 k 2 , with the squared 
mass A 2 = A 2 — 4d 2 Q 2 , and the velocity c 2 = dQ 2 . This spectrum becomes massless for 
A 2 = 4d 2 Q 2 . Evaluating the determinant of G we get 

2dQ 2 / 1\ 2 . //W 



'MF 



J, - A ( K + D + ^? teg2s " lh (^)' < 3 ' 42 > 



Here cjfc = ^A 2 — 4rf 2 <5 2 7 2 3 , and A is assumed to be positive. Note that the transition from 
Eq. (3.36) and Eq. (3.39) to Eq. (3.42) included performing a Matsubara sum. We stress 
that a naive summation will not lead to the +| in the second term of Eq. (3.42). This 
factor is received only by performing the determinant of the inverse propagator exactly, 
with the correct discretization of Euclidean time given in Eq. (3.21). We give a detailed 
explanation in Appendix B, where we derive the following rule of thumb for evaluations of 
a two-sublattice determinants of the form (3.37): Evaluate the determinant "naively" by 
replacing the partial derivative d T with its simple continuum form —iu> n , and multiply the 
result by a factor of e /3A / 2 for each degree of freedom. For example, in the case considered, 
the number of degrees of freedom is NmN s , hence the factor | in Eq. (3.42). 

We proceed to write the mean field equations for the action (3. 42), 5 

dSyiF 9Smf 



d\ ' dQ 



0, (3.43) 



that become 



Q - 2 -«^ W + i), (3.44) 



2dJ 1 N s ^ u k \ y J 2 

1 2 _ A / . . 1\ . 
" + 2 - N.Z- (»<<*> + 5 ) ■ <*«> 

5 Note that in principle one should also check whether the point (Qmf,Amf) is a minimum. If it is 
not, however, one can take some of the variables (in this case Q or A) to be imaginary, and make all 
curvatures positive. This is actually the reason for the ansatz (3.34) 
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Here n(oo) is the Bose function (e^ — 1) 1 . From this point on we restrict to T = 0. The 
first solution (solution A) to Eqs. (3.44)-(3.45) is 



Qa = 0, 
2 

\ A = -coth -1 (2re+ 1). 



(3.46) 
(3.47) 



This gives Q A = X A = at T = 0, and a value of S M f = -p' 1 log 
The second solution has Q ^ and Eq. (3.45) becomes 

" nV) + N s /2A' 



K + l 



— 5 — 7 



(3.48) 
(3.49) 



with i] = X/(2dQ). Here the momentum integration was replaced by an integration over 
the function 7^, provided with a integration measure D(j) given in Fig. 3.1, normalized 
to/d 7 D(7) = l. 

Note that we have separated the fc = (7 = 1) mode from the others, to allow a 
Bose-Einstein condensate (BEC). The second term in Eq. (3.48) will be non-zero at the 
thermodynamic limit only if n(A)/\/l — i] 2 ~ 0(N S ), which in turn means that 77 — > 1. 
Note that since 7 2 < 1, the range of r\ is [0, 1]. The function f{rj) is given in Fig. 3.2. It 
is monotonically increasing, and at 77 = 1 takes the following values: 



/(I) 



00 d — 1, 
0.19 d=2, 
0.0787 d = 3. 



(3.50) 



We divide the solution of Eq. (3.48) into two cases. 
Disordered phase (solution B) 



Here one assumes that A 7^ at the thermodynamic limit and drops the second term in 
the right hand side of Eq. (3.48). The solution then reduces to 



K = f{r} B ), 
Xs/dJi = J d-yDia) 
Qb = riEXs/^d. 



r 



\J\ — 77 2 7 2 ' 



(3.51) 
(3.52) 
(3.53) 
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Figure 3.1: The integration measure of 7 for d = 3. 



This solution is guaranteed as long as k < /(l). Recalling Eq. (3.50), this solution exists 
always for d — 1, but only for k < k c , with k c ~ 0.19 for d = 2, and k c ~ 0.0787 for 
d = 3. The MF action of this solution is -4dQ 2 B /Ji < 0. 

Ordered phase (solution C) 



For /t > /t c , Eq. (3.48) cannot be solved without the extra condensation term. A consistent 
solution will be 

Ac = A B (l) + dJi(/c-K c ), (3.54) 
Qc = A c /2ci, (3.55) 

= \h - AT , J ( 3 - 56 ) 

V N s \ c (k - K c ) 

Here A^(l) means the value of Eq. (3.52) at rj = 1. The action of this solution is 
—AdQg/Ji, so solutions B and C have lower action than solution A, and are continu- 
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ous at k — k c . Nevertheless, the derivative of the mass A at not continuous. 

Recalling that k = N c /4Nf, one finds that for d — 3, and at T = there is a Neel phase 
for ^ > 0.31. This means that for N c = 3, and Nf < 9.6 the U(4Nf) is spontaneously 
broken (for calculations of the staggered magnetization that show that the BEC indeed 
corresponds to a Neel phase [45, 50]). 



We summarize the zero density problem in the schematic phase diagram in plane N c —Nf 
in Fig. 3.3, borrowed from [45]. The large- Nf phase is disordered, and will not concern 
us. The location of the phase boundary is a line of constant slope, which turns out to 
be the line N c /Nf = 0.31; this means that QCD is safely in the ordered phase for any 
reasonable number of flavors. 
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Large N c sigma model 
(Chapters 4-6) 



Ordered phase 




Disordered phase 



Large N f 



Figure 3.3: Phase diagram for the nearest-neighbor antiferromagnet at T = 0, and zero 
density in N c -Nf plane, for MF theory with Schwinger bosons. 



3.4 Mean field theory at non-zero density 



In this section we investigate how a non-zero density of baryons influences the MF equa- 
tions and changes the boundary of the Neel phase in the N c -Nf plane. The MF approach 
discussed above can be used to investigate the following configuration of non-zero baryon 
number. On M of the N s lattice sites we replace the U (N) representations with the singlet 
of SU(N), that has B = 2Nf. This means that the Q 1 ^ operators in these sites are zero 
(except Q° which corresponds to baryon number). The density p of these configurations 
is p = mc, where 

M 

(3.57) 



c = 
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is the concentration of the singlet representations. In the presence of the "impurities", 
the effective Hamiltonian (3.1) suffers from a set of missing links and becomes 

H =^k £ QnQn+pi- (3.58) 
n,n+(ig{r} 

{r} denotes the locations of the M impurities. We again make the MF ansatz (3.33)- 
(3.34) and turn to calculate the MF action and MF equations. First note that cN s sites 
do not participate in the Hamiltonian. Also if there are no neighboring impurities, then 
there are 2cdN s missing links. This means that the first two terms in Eq. (3.36) are now 
multiplied by (1 — c) and (1 — 2c). 

Q 2 -> (l-c)Q 2 , 

k\ -> (1-2c)kA. (3.59) 

Also, the evaluation of the functional determinant is different. We have analyzed this 
problem using two techniques. The first is an analysis of certain configurations of impuri- 
ties, in which the system does not lose translation invariance. We solve only the simplest 
cases of c = |, \ in section 3.4.1. In section 3.4.2 we make a quenched average of the 
impurities' positions which restores translation invariance. In this method we find that 
the corrections to the MF action for low values of c can be written as a power series in c, 
and we the present results of the first order. 

3.4.1 Periodic density with baryon concentration c = |, | 

In general a random distribution of baryon impurities destroys translation invariance 
found at zero baryon number, and makes the methods used in the previous section in- 
applicable. In particular, a Fourier transform does not diagonalize the propagator of 
the bosons, making the evaluation of the determinant very difficult. There are however 
densities that can be realized with a translation invariant configuration of impurities. In 
three dimensions the simplest configurations have c = |, and c = \. In this section we 
calculate the MF action of these configurations, and proceed to solve the corresponding 
MF equations. 

Concentration c—\ 



We divide the lattice into unit cells of the form shown in Fig. 3.4. Each includes seven 
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3 1 

Figure 3.4: A unit cell for impurity concentration of c = |. 



ordinary sites and one impurity. We place the latter on the upper-left-back corner of the 
cell, and denote the bosons on the other sites by z njP . Here n is the location of the cell in 
a simple cubic lattice of spacing 2. The index p = 1, . . . , 7, is an internal index denoting 
the sites in each cell, corresponding to the numbering in Fig. 3.4. Using these notations, 
the MF action is 



■S'mf 
detG 

S'iQ) 



^--A/t+^-Tr logCT 1 , 

2Ji 8 N s (3 B 

[ DzDz* exp - [drJ2 z * np (~d T + A) z np - S\Q) 

— Q dT*^2 [ z n,2{Zn,l + z n-x,l) + z n,z{ z n,l + z n-y,l) 
J n 

+ z n,4:( z n,l + z n-z,l) + z n,b{ z n,2 + z n-z,2) 
+ z nfi( z n,3 + z n-z,Z) + z nj( z n,2 + z n-y,2) 
+ z nj( z n,3 + z n-x,z) + z n,5( z n,4 + z n-x,i) 
+ z nfi{ z nA + + CC. 



Next we introduce the Fourier transform 



-n,p 




^ ' z kn,p& 
k,n 



i(n-k—w n T) 



(3.60) 
(3.61) 



(3.62) 



(3.63) 



where k belongs to the first Brillouin zone of the simple cubic lattice. In momentum 
space S'(Q) becomes 

S'(Q) = -\ (QQ z l T fc z^ k + c.c.) , (3.64) 
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where z and T are 



Zk = 



V Z k,7 J 



( o 

Y* 
Z* 



V 



X Y Z \ 

Z Y 

Z X 

X Y 



and Xi = |(1 + e" ifc 



Z* X* 
Z* Y* 

o r* x* o ooo 

. We find that the eigenvalues 7 P of T k obey 
= 1 [ 

i 

ai = 3Epf + 4E^P?-2Re((Xn 2 + ^(^ 2 + ^ 



(3.65) 



7 |7 6 + fl 27 4 + Oi7 2 + a 



i>j 



E PfP? - E A* - + 2 Re [(Xr) 2 ^ + - p*) 



+(x 2 -Y*)( P y x -py y ) + (x*+Y>)p* 



(3.66) 
(3.67) 

(3.68) 



(3.69) 



with pf = \Xi\ 2 . We have verified that the solutions to Eq. (3.66) are real and obey 
7o = 0, 7 2 = —73, 74 = —75, and 7 6 = —77, and are all even in the momentum k. Since 
the first term in the action of Eq. (3.61) is unity in the 7-dimensional internal space, we 
write 



detG 



= f DzDz* exp - E / dT ( z l P z -k,P ) ( 



<9 r + A 



-d T + \ 



%kp 

7* 

z -kp 



(3.70) 



Here the factor 1/2 of Eq. (3.64) is replaced by summing over only half of the Brillouin 
zone (the remainder of the Brillouin zone gives the same contribution due to the coupling 
of Zk and z-k). As Eq. (3.70) stands, all the fields are independent and can be integrated 
over. This gives 



'MF 



9Cf 
2Ji 



7 ^ ( 1 



-E 1 



d^Dpipf) log 2 sinh 



(3.71) 
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Figure 3.5: The integration measures for c—\. 



where 0^(7) = \/X 2 — 36Q 2 7 2 , and -D p (7) is the integration measure of 7fc iP given in Fig. 3.5, 
normalized to / d^Dpij) = 1. The factor | accounts for the N s /S sites found in the simple 
cubic lattice. Also note that the third term, — |A|, is added according to the rule of thumb 
prescription of Appendix B. 

We proceed to concentrate on the MF equation 



MF 



0, 



(3.72) 



which, excluding BEC, becomes 



1 {K + i ) = l („(A) + I) + li;J dlDfb)7T ± w („(.) + I) . (3^73) 
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At T = this simplifies to 

In Fig. 3.6 we show the right hand side of the MF equations (3.51) and (3.74) for compar- 
ison. First note that for c = |, rj takes values up to 1.14. This is because the maximum 
value 7 can have is 0.877, and not 1 as for c = 0. Next we find that the right hand side 
has a maximum of 0.0881 at i] = 1.14. This is the point where BEC occurs. The physical 
result is therefore that the Neel phase at c = | shrinks, and symmetry restoration occurs 
above n c = 0.0881. For QCD with N c = 3 this means that the symmetry is restored 
already at Nf ~ 8.5 instead at Nf ~ 9.6 at c = 0. 



\/l — ?7 2 7 2 



(3.74) 




Figure 3.6: Right hand side of the MF equations Eq. (3.51) and Eq. (3.74), with c = 
and c = |. The solution t](k) is given by intersection of the curves with horizontal lines. 
Note that the curves have an end at rj — 1 and rj = 1.14, marked by the full circles there. 
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Concentration c—\ 



In this case we put an impurity on site no. 1 in Fig. 3.4 as well. This means that c = |, 
and Eq. (3.60) is changed to 

Smf = ^ - \\k + j^Tt \ogG-\ (3.75) 

In evaluating the last term, we again use Eq. (3.61), but with the boson field z n> i omitted. 
This means that the matrix T k in Eq. (3.65) has lost its first row and column and is now 
6-dimensional. We find that its eigenvalues obey 

7 6 + a 2 7 4 + ai7 2 + a = 0, (3.76) 

with 

a 2 = -2j2p 2 i, ai=(j2p^j , «o = -^P 2 x pIp 2 z . (3.77) 

Again we find that the solutions 7 P are real and obey 71 = —72, 73 = —74, and 75 = —76. 
At T = 0, the MF equation becomes 

p=l,3,5 

with the measures given in Fig. 3.7. 

In Fig. 3.8 we show the right hand side of the MF equations (3.51) and (3.78) for com- 
parison. In this case, because the maximum value of 7 is 0.66, r\ takes values up to 1.5. 
We find that the right hand side at the BEC point rj = 1.5 is 0.117. The physical result 
is therefore that the Neel phase at c = | shrinks even more than in c = |, and symmetry 
restoration occurs above n c = 0.117, or for QCD with N c = 3 already above Nf ~ 6.4. 

3.4.2 Random baryon positions — the quenched approximation 

In this section we solve the MF equation of a quenched random distribution of baryon im- 
purities. We adopt the formalism of [74] developed to treat the affect of magnetic impuri- 
ties on an SU (2) magnet. Here the analog to the impurities are sites with SU(ANf) singlet 
representations. We denote the number of of singlet representations on the even and odd 
sites as M e and M Q , and present calculations to first order in the concentration M e /N s , and 



- 1 



(3.78) 
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Figure 3.7: The integration measures D p {p(). The oscillations in D 3 are an artifact of the 
resolution of the momentum mesh. 
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0.04 



0.02 



Figure 3.8: The right hand side of the MF equations Eq. (3.51) and Eq. (3.78), with c = 
and c = \. Note that the curves have an end at r\ — 1 and r\ = 1.5, marked by the full 
circles there. 



M a /N s . Again the tricky part is to calculate the functional determinant (3.37), because 
now the boson propagator is no longer diagonal in momentum space. 

We begin by modifying Eq. (3.37) as follows, 



det G(Q, A) 



J (DzDz*) ne{r} e- s '^ = J DzDz^exp{-S -S lmp }, 



where S is the 'pure' action of Eq. (3.37) 
S = I drJ2 

J n 

The auxiliary action S' is 



Z n (~ ®t + A) Z n — Q (z n Z n+ p, + Z^Zn + {i) 
fi>0 



S' = Jdr £ z * n (-d T + X + X') Zn , 



(3.79) 



(3.80) 



(3.81) 



ne{r} 
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where A' is real and obeys A + A' > for stability. Also we choose the action S' to 
be defined in the continuum temporal limit, i.e. with the "naive" definition of d T (in 
Matsubara space it is just —iuj n ). This means that the expression in the square brackets 
in the left hand side of Eq. (3.79) is given by 



logZ' = -ciV s log 2 sinh 



/3(A + A') 



(3.82) 



The A' term was inserted in order that Si mp will have the following form, 



'Imp 



ne{r} 



±d 



A Z* n Z n + Q ^2 {znZn+ix + 

|U=±1 



" nrt -m 



z* n V nrn z m . (3.83) 



We now write the MF action for the Hamiltonian (3.58) as 



S MF = (l-2c)^- C A(« + i)-clog(l- e -^')) 



Ji 

Tr logG _1 + AS, 



AS 



N s (3 



Tr log (1 + GV) . 



(3.84) 
(3.85) 
(3.86) 



Here G _1 (Q, A) is the same inverse propagator that appears in the zero density problem 
in section 3.3, and we proceed to evaluate AS. Using the Fourier transform (3.38) we 
write V in momentum space as 



"Simp - EE fk,n V k,k'fk',n, 



(3.87) 



u n k,k' 



Vfc,fc' = 2dQ — 



E e 

m<a{r e } 



b 2 



:(k'-k)m u Ik' , v-^ 



J(k'-k)m ( 7fc 

{ l'k V 



(3.88) 



Here fk, n is given in Eq. (3.40), m are the positions of the impurity sites in terms of 
the original simple cubic lattice, {r e } and {r Q } denote the impurities on even and odd 
sublattices, and b 2 = X'/2dQ. Fortunately, the vertex V^w is reducible to the form 6 



V 1 \ ^ -ye ire „i(k'—k)m . ^ \ ^ -7-0 rro „i(k'—k)m 
k,k' - Jj- 2^ L kJk' e + Jj~ 2 fc<^fc' e > 



(3.89) 



s m,e{r} e 



s me{r} 



3 This reducibility is achieved only for 6^0, and is in fact the reason we introduced A'. 
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T " = ^ ( llh -llh ) ' Ji = ^ ( o ) ' (3 ' 90) 
2! = v^Q ( " 7 fe/6 7 f ) ■ •» = ( "iJil I ) ■ < 3 ^) 

With the above definitions the evaluation of AS* becomes tractable, and takes the form 
of a power series in the impurity concentration. We proceed to write 

N s pAS = ~ — - — Tr (GV f, (3.92) 
P =i P 

where we implicitly assume that the sum (3.92) is convergent. The trace is over the 
momentum and internal even-odd indices. Using Eqs. (3.90)-(3.91) and Eq. (3.39) we 
find that Eq. (3.92) contains contributions of the form 

Tr Al 1 Al 1 ---Al/ [(kl - k2)mi+ -- Mkp - kl)mp \ (3.93) 



ki ,. .. ,kp 
m 1 ,...,m p e{r} 

Here A% can be any one of the four matrices 

A ee = ^G k ll A°° = ^ 3° k G h Tl (3.94) 

A eo = jj-JRGJZ, A° e = ±-JZGJZ. (3.95) 

A diagrammatic method is useful to treat the evaluation of Eq. (3.93), and we define the 
following Feynman rules. 

1. The propagator G given by Eq. (3.39) is represented by Fig. 3.9. It is diagonal in 
momentum and Matsubara frequency. a,b = e,o are the even-odd indices. 

(G k M) ab = a » b 

k 

Figure 3.9: The boson propagator. 



2. The interaction are given in Eq. (3.89) is a sum of M e + M Q vertices, one for each 
of the impurities. The vertex that corresponds to the impurity on site r is given in 
Fig. 3.10. 
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Figure 3.10: The vertex of the impurity that resides on the site r. Here X, J are defined 
by either Eq. (3.90) or Eq. (3.91). 

fci 



Tr J2 *2 

ki,k 2 ,k 3 

Figure 3.11: The loop diagram that corresponds to the expression in Eq. (3.93) with 
p = 3. 



3. The expression (3.93) is represented by the diagram in Fig. 3.11, which describes 
the scattering of the bosons off three impurities. Note that this diagram stands for 
all possible ways to scatter three times off M e + M Q available impurities. 

4. We adopt the notations in [74], and write the diagram in Fig. 3.11 as the sum in 
Fig. 3.12. A diagram with q external points (each denoted by an x) corresponds to 
a process in which the bosons scatter off q different impurities. The lines connecting 
the external points with the vertices denote which of the impurities was hit in a 
given vertex. For example in the leftmost diagram, the boson collides three times 
off the same impurity, while in the rightmost diagram, it collides three times off 
three different impurities. Finally note that each of the diagrams in Fig. 3.12 stands 
for all possibilities to choose the impurities that the bosons scatter off. 

The next step is to use the quenched approximation, and average Eq. (3.93) over the 
position of each impurity. This assumes that the impurities cannot move on the time 
scale that the bosons move, and is exact in the order we work to within the strong 
coupling expansion. We define the averaging procedure over the locations of the set {r} 
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X X X X X 




Figure 3.12: The loop diagram that corresponds to the expression (3.93) with p = 3. The 
leftmost diagram describes a boson scatters three times off the same impurity, while the 
rightmost diagram describes a boson that scatters three times off three different impurities. 



as 

/ O \ M e +M N 3 /2 

(0) = ( w ) E 0(({r e };{r }). (3.96) 

K,o}=i 



This relies on the assumption that the concentration of impurities is sufficiently small 
that there are no links between any two impurities. This assumption is invalid at high c. 

Applying the averaging procedure to each of the diagrams in Fig. 3.12 results in the 
following power series for expression (3.93), 

Min[p,M] q 

N s E E W e B< Me) . (3.97) 

q=l q e =0 

Here q e and q — q e are, respectively, the number of external points that correspond to an 
impurity residing on an even or odd site, and c e o are defined as 

c a = M a /N s , a = e,o. (3.98) 

We proceed to prove this point. First we denote a diagram with p vertices and a total 
of q external points, out of which q e represent impurities on even sites, by 2^ ge - s is 
an additional index that denotes the different diagrams that have the same values of p, 
q, and q e . For example for p — 4, and q = q e = 2, s takes the values 1,2, and 3, and 
corresponds to the three diagrams in Fig. 3.13. 

Looking at the structure of the exponentials in Eq. (3.93) we see that when we evaluate 
T^pqq e , where m.j take q different values, the exponential turn into 

e i{k' 1 ri+k' 2 r 2 +...+k' q _ 1 r q -- L ) _ (3.99) 
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X 
I 



X 



i \ 





Figure 3.13: The loop diagrams that corresponds to p = 4, q = q e = 2, and s = 1, 2, 3. 



Here {fc^} are linear combinations of the momenta that flow in the diagram. The averaging 
procedure is in fact applied only to the exponentials (3.99), and results in q — 1 momenta 
conservation rules, 

( e i(fcir 1 +*4r 2 +...+fc' e _ 1 r,_ 1 )) = £ fe ,£ fc , . . ,g (3.100) 

Substituting Eq. (3.100) into Eq. (3.93), and taking into account of the number of pos- 
sibilities to choose which impurities are represented by the external points in we 
have 

WD ={ M q e ){ q M °q e ) T * £ < ' ' - " (3-101) 

The next step is to assume that the terms in Eq. (3.92) with p > M ej0 are negligible. 
(M ej0 are of the order of N s ) Recalling that A% oc iV" 1 , we find the following continuum 
limit of Eq. (3.101) for p, q e , q — q e <C M e o . 

/ \ M qe M q ~ qe { r \p-i+ 1 , s 



Qe ] - (q~q e 

{,c e y*{ Co y° 



= 4ifFi>i(/*) F «* ({fe}) ' (3 ' 102) 

where the (p — q + 1) momentum integrals are over the Brillouin zone and are finite. We 
can now identify that 

(r ) qe (r ) q ~ qe / r \p-i+ 1 , s 

Si = gj(g - ge)! ? (Z ^ F - Lm) > (3 - 103) 
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for p < M Ct0 . Q.E.D 



MF action to 0(c) 



We now evaluate AS* to first order in the concentration c. We therefore take into account 
diagrams with q = 1 only. In that case q e = 0,1, and s = 1 (there is only one possible 
diagram for each value of q e ). In this case the diagrams that contribute to the p-th term 
in Eq. (3.92) correspond to p scatterings off the same impurity, which resides on either 
sublattice. The structure of Eq. (3.92) is now simple. For q e — 1 all A 1 = A ee , and for 
q e = all A i = A°°. The result is 



N S (3AS 



P 



5(1,1) , g(l,0) 



= ^EE c « Tr E 

p<M e , D 



(_1)P+1 



u n a=e,o 



p 



jv s EE^io g 



i + 



/( 



V 27r , 



(3.104) 



We have verified that the contribution of the Tr is independent of a, so AS becomes 



(3AS(\, rj;b) = c E log \ det 



1 + J d 7 D( 7 )- 



7/ 



+ e 2 - rj 2 -f 2 



b 2 {l - ie) + 2r^ 7 2 + 7 2 (1 + ie n )/b 2 - 7 2 (r? + (1 + *e n )/6 2 ) 
1 2 ( v +(l + te n )/b 2 ) _(l + ien ) 7 2/ & 2 



(3.105) 



Here the overall concentration c = c e + c D = M/N s , e n = u> n /\, and the momentum 
integration was again replaced by an integration over the variable 7 with the measure 
shown in Fig. 3.1. We have calculated AS(X, rj; b) and confirmed that for \r)\ < 1 it is real. 
Substituting Eq. (3.105) into Eq. (3.86) we have 



Smf — 



2rf(1 - 2c) -g 2 -(i- c) A. 



Ji 



+ 



( d 7 .D( 7 ) log 2 sinh 
p J 



Mi) 



5 A(l-c) 



c / n . , /3A(l + r/6 2 ) 0/1 A 
+ - ( - log 2 smh ^ 1 2 ; + AS(X, m b)J. 



(3.106) 
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The first row of Eq. (3.106) is the result of (3.59). The second row is Tr logG -1 for 
the pure system. (Note the term proportional to (1 — c) is exactly the rule of thumb 
prescription from Appendix B) The contributions to the third row are the continuum 
evaluation of AS and — log Z', evaluated in the continuum limit of Euclidean time. We 
have verified that the contributions in the third row do not depend on the value of the 
unphysical constant b. 



MF equations to 0(c) 



We concentrate here on the T = case. There we have 



jV s 7Vm/T MF J, 

= 2 / dlD ^ V 1 -V 2 ! 2 + c f Hm 



2d(l - 2c) 



Q 2 -\ «(l- c ) + - )+XJ r (v), 



—AS(X,i]; b) 




(3.107) 



(3.108) 



Where we have evaluated the (3\ — > oo limit of the second term in F(rj) as an integral 
over e. In order to investigate the boundary of the Neel phase with non-zero value of c, 
we concentrate on the equation 

9SmF = 0, (3.109) 



and assume no condensates. This becomes 



1 



K 



[f(v) +c£loop] 



Here f(rj) is the right hand side of the MF equation for c 
Fig. 3.2. Q\ oov , includes all other loop contributions. 



(3.110) 

0, Eq. (3.48), as plotted in 



We plot the right hand side of Eq. (3.110) in Fig. 3.14 and see that as the concentration c 
increases, the maximum value of k decreases. This means that the Neel phase is expanding, 
which appears to contradict the results of Section 3.4.1. 

In order to resolve this problem we compare the results of the quenched calculation for 
c = | and c = | with the exact diagonalization in Section 3.4.1. We present the results 
of the two calculations in Figs. 3.15-3.18. The results are very similar. We see that the 



errors of the quenched results are larger for c = \ than for c 



This is presumably 



because of the restriction to first order in c. Also we note that the quenched perturbative 
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approximation is inapplicable for rj > 1, since the unperturbed ground state does not 
admit such values of 77. Nevertheless, the result of Section 3.4.1 lets rj be larger than 
1, and in fact it is in this extra interval that the right hand side of the MF equation 
continues to grow, and the Neel phase shrinks. For both c = | and c = |, there is a value 
of ^max > 1 at which BEC occurs. This regime of MF parameters is unfortunately not 
accessible in the quenched approach. 

To conclude we see that the MF quenched approximation in first order in c does not 
reproduce the MF theory near the phase transition. The deviation between the quenched 
approximation and the exact diagonalization grows as one approaches the BEC. This 
approach is not applicable to extract information about the phase transition, simply 
because it is a perturbative approach whose starting point is the c = system. One can 
use this approach to extract quantitative properties of the disordered phase, and to see 
how these change when baryons are included. For example, the propagators of bosons 
will then obtain a self-energy contribution also calculable via the diagrammatic tools 
described above. It is expected that this will show that the baryon impurities slow down 
the bosons, and even make them decay. This is in analogy to what happens in condensed 
matter magnetic systems with lattice vacancies [74, 75, 76]. 

Finally we note that although this MF analysis was devoted to the regime k < k c and 
explored the disordered phase (while emphasizing the location of the phase transition), it 
is also possible to investigate the ordered phase. This can be done in the large N c limit 
by using a Holstein-Primakoff representation of the Q n operators (which is presented for 
m — 1 in Appendix H). The quenched approach can be also applied in this limit, and in 
analogy to [74] we expect that the Goldstone bosons in this phase (the pions and kaons 
in our system) will have a smaller velocity, and will obtain finite decay widths. 
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Figure 3.16: The MF effective action for c = | in the quenched approximation, and the 
exact diagonalization. 




Figure 3.17: The right hand side of the MF equation for c = g in the quenched approxi- 
mation, and the exact diagonalization. 




Figure 3.18: The right hand side of the MF equation for c = \ in the quenched approxi- 
mation, and the exact diagonalization. 



Chapter 4 



Nonlinear sigma model and the 
classical ground state 



In this chapter we use a coherent state basis to write the partition function of the Hamil- 
tonian (2.69) as a path integral. The result is a nonlinear sigma model. We first analyze 
the properties of the classical counterpart of a toy sigma model that has N = 3. We 
use mean field theory to find its ground state as a function of temperature for a set of 
baryon number configurations. Next we approach the full quantum sigma model. The 
action of the model has an overall factor of N c , and is tractable in the large N c limit. We 
take this limit and examine the classical ground state for a set of uniform baryon number 
configurations. 



4.1 Coherent states and sigma fields 

We employ a generalization of spin coherent states [51] to derive a path integral for the 
spin model of i? eff . We recall that a given site carries generators Q n n of U(N) in a 
representation with N c columns and m rows, with B n = m — 2Nf. 

First we choose a basis for the Lie algebra of U(N). This consists of the generators Sp, 
with a, (3 — 1, . . . , N, whose matrix elements in the fundamental representation are 

(S^S = SajSps- (4.1) 
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The corresponding charges are 

a Z 

= Et^-k", (4.2) 

a Z 

where we have subtracted a constant for convenience. The Cartan subalgebra consists of 
the operators 

H a = Q a a . (4.3) 
We build the coherent states from the state of highest weight |^ ) 

m 

l*o>= II «i< 2 ---<^ c )|Dr), (4.4) 

a=l 

where the numbers are the color indices. This state is an eigenstate of the Cartan gener- 
cttors 

H °\* o) ~ { -(N c /2)\*o) fora = m + l,...,N. (4 " 5) 
The expectation values of the generators in this state obey 

(*o\Qp\*o) = \N c A a p, (4.6) 

with 

A = ( V -iL ) ■ (4J ) 

In order to obtain the other states in the rep, we perform the following U (N) rotations 
on |^ ), 

(m N \ 
E E W-"?Ql))\9 ). (4.8) 
A=l /u=m+l / 

The only generators that appear in Eq. (4.8) are those that lower an H a that starts 
from N c /2 in Eq. (4.5) while raising another H a that starts from —N c /2. Any other 
generator would annihilate |^o) and thus give no effect in the exponential. 

The coherent states are normalized, 

(a\a) = 1, (4.9) 
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and over-complete, 

fdfji(a) \a)(a\ = 1. (4.10) 

In Eq. (4.10) the integral is over the coset space U(N)/[U(m) x U(N — m)} (see below). 
Matrix elements of the generators are given by 

H W = \N c a aP) (4.11) 
where the matrix a a p is given by a unitary rotation from A, 

a = U(a)AU( a y. (4.12) 
The matrix U(a) is built out of the m x (N — m) matrix a*, 



U = e A = exp 




(4.13) 



a is both Hermitian and unitary. 



The manifold of matrices a is the coset space U(N)/[U(m) x U(N — m)], a sub-manifold 
of U(N). This is because for any matrix U(a), one can generate an orbit U(a)V by 
multiplying with a matrix 

where X e C/ (m) and Y & U (N — m) . All matrices in the orbit will give the same matrix 
er when inserted into Eq. (4.12), and thus in integrating over the configuration space of a 
one must choose only a single representative of each orbit. This set of representatives, the 
coset space U(N)/[U(m) x U(N — m)], is the quotient space of the non-invariant subgroup 
U(m) x U(N — m). Basically the origin of this manifold is the fact that the quantum 
state |\&q) is invariant (up to a phase) to U(N) rotations that belong to the subgroup 
U (m) x U(N — m). This can be understood by the form of the Young tableau in Fig. 2.4. 
This representation of SU (N) is a singlet of a certain SU (m) x SU (N — m) subgroup, 
more precisely, of the rotations that do not involve indices from the [m + 1, N] range and 
that leave the [l,m] indices of the state invariant. 

The measure over the coset space must be invariant under unitary rotations, 

\a)^R(V)\a), (4.15) 
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where R(V) represents the rotation V in Hilbert space. Equation (4.11) shows that this 
is a rotation 

a^VaV j (4.16) 

and by Eq. (4.12), this means that a measure in U must be invariant under U — > VU. This 
fixes the measure uniquely to be the Haar measure in U(N), and thus one can integrate 
over the coset space by integrating with respect to U over U(N) and using Eq. (4.12). 

A representation whose Young diagram has N — m rows is the conjugate to the represen- 
tation with m rows. Its coherent state space can be constructed to look the same, with 
only a sign difference. To do this we start with the lowest-weight state, which satisfies 
[cf. Eq. (4.5)] 

' ; I (N C /2)\V ) for a = m + l,...,N. 1 } 

This introduces a minus sign into Eq. (4.6). The subsequent steps are identical, with only 
the replacement of Eq. (4.11) by 



lfl\Q$\a) = --N c a aP . (4.18) 
Here, too, a is given in terms of A and U by Eq. (4.12). 

Using the Taylor expansion of the exponentials in Eq. (4.13) and noting that the expo- 
nential argument A anticommutes with A we have 1 



a(a) 



( , v sin (2v / aJa) \ 

cos(2vaaM —a - 



sin ( 2v a)a 



V 



,t 



aa 



t 



cos 



V at a 
(2v / a^ 



(4.19) 



Here a) a is a square matrix of dimension N — m and aa) is a square matrix of dimension 
m. Using the notation 

sin (V at a) 



at a 



we have 



u{4>) 

a{4>) 



^t yr^ 



1 This generalizes a parametrization found in [77]. 



-20^1 - 0t0 

-2^1 - 0t00t _! + 2 0t0 



(4.20) 

(4.21) 
(4.22) 
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The representation given in Eqs. (4.21)-(4.22) is useful to treat fluctuations of order 1/N C , 
as will be seen in next chapter. Also for m = N — 1 it coincides with the form given in 
[49] , which in turn can be used to infer a generalized Holstein-Primakoff realization of the 
U(N) generators (the <fi operators then become the Holstein-Primakoff bosons). We note 
in passing that another useful representation given in [78] is related to Eqs. (4.21)-(4.22) 

4> by identifying — a tan VaJa/ \J a) a — <f>/ J 1 — 



4.2 Partition function and action 

The partition function Z = Tre _/3// with (3 = 1/T, can be written as a path integral by 
inserting the completeness relation (4.10) at every slice of imaginary time. This gives 

Z = J Da exp -S, (4.23) 

where the the action is 



f 13 

S= dr 
Jo 



l-(a(r)\a(r + dr)) + } 
dr 



(4.24) 



The Hamiltonian H(a) is a transcription of the quantum Hamiltonian to the classical a 
matrices. Starting with the quantum operator we have 

Ql = ^ n M^ n = M2^ n S^ n 

= M2pQ a p (n) + ±N c TrMV. (4.25) 

Expressed in these variables, the quantum Hamiltonian is 2 

H$ = E JiQlQl** « +1 • (4-26) 

where Jj = (2/N c )K(j). We transcribe this according to Eq. (4.11) to obtain the classical 
Hamiltonian, 

V / nt± 

3^0 

2 The B' term from Eq. (2.51) indeed disappears. We have dropped an additive constant that is 
independent of B 1 . 
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where 



al = Tr M^Or,, 



(4.28) 



Recall that each a n is an N x N matrix ranging over the coset space appropriate to site 
n. We find it useful to write the interaction as 



H(o) 



1 fN r \ 2 



^ Tr cr n cr n+ jfj, + ^ JjTr <j n a^cr n+ jfj. 



a. 



(4.29) 



j'Godd j Seven 

where the three Dirac matrices have replaced the sign factors s^. 

The time- derivative term in S is a Berry phase [46]. It can be expressed as follows 

Sb = / dr duTr \(j(t,u) <9 u <t(t, iz) d T a(r,u)] . 

A Jo Jo 

Here a(r,u) = U(t,u)AW(t,u) with 

U(t,u) =e uA ^ T \ 

and obeys 

a(r,0) = a(r',0), Vr,r', 
<t(t,1) = <t(t), Vr, 
a(0,u) = cr(/3,u). Viz, 



(4.30) 



(4.31) 



(4.32) 
(4.33) 
(4.34) 



This means that the t-u plane is a disk with unit radius, r plays the role of an angular 
variable, while u is the radial coordinate. Using Stokes theorem Eq. (4.30) is brought to 
the form 

Sb(U) = jy^Tr (AU\t) d T U(r)) . (4.35) 

It is interesting that Eq. (4.35) gives the same result even if U is a general U(N) matrix, 
not restricted to the form of Eq. (4.13). To see that we write U'{t,u) = U(t,u)V(t,u). 
Here U' E U(N), U is given by Eq. (4.31), and V E U{m) x U(N - m), 



V 



exp (uA£\t)) 

exp (uA 2) (t)) 



(4.36) 



U' and therefore V are continuous functions on the disk. Substituting U' into Eq. (4.35) 
one finds that 

AS B = S B (U') - S B (U) = yJ^ ( Tr " Tr KO) ' ( 43? ) 
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This means that ASb is a an integer multiple of inN c . A non-zero value of ASb is a result 
of winding the U(l) factors of V on the disk boundary u — 1. Such winding is however 
forbidden, since it means that V must have a discontinuity inside the disk. In fact the 
continuity of V means that det(V) must also be continuous, so \/u G [0, 1] one has 

= J drd T det V(r, u) = J dr d T (log det exp (uA^) + log det exp (uA^)) . (4.38) 

Substituting Eq. (4.38) with u — 1 into Eq. (4.37) we find that ASb is an integer multiple 
of i2irN c . Since N c is an integer, exp ASb = 1, and expS B (U) = exp S B (U'). 3 



Finally the lattice action is 



- Ill f 

2 Jo 



dr 



^krJJld^ + ^-H^r)) 



(4.39) 



A„ will vary from site to site if m does. If one takes the route of Eq. (4.18) for a 
representation with N — m rows, then the kinetic term for that site acquires a minus sign 
(see below). 

The number of colors has largely dropped out of the problem, since a is just an N x N 
unitary matrix field. The explicit factor of N c /2 that multiplies the action invites a 
semiclassical approximation in the large- iV c limit which we do in Section 4.4. This of course 
neglects the ^-dependence of the interaction term, but we recall that for N c 1, the 
couplings of the quantum Hamiltonian in Chapter 2 are Jj ~ 1/ (g 2 N^). This means that 
the couplings of the classical Hamiltonian in Eq. (4.29), multiplied by 2/N c in Eq. (4.39) 
scale like l/g 2 N c , which we then fix when taking iV c — > oo. 

We conclude this section by noting that the form of the Berry phase of Eq. (4.39) is 
exactly the reason why the representation of a and U in term of the variable <j) given in 
Eq. (4.21) becomes useful. In terms of this variable the kinetic term becomes 

N r r i 

Sb = \ dr ]T Tr . d T <f> n - 8^1 ■ 4> n ] . (4.40) 

This simple quadratic form is not obtained by expanding U in terms of a or Z. 



4.3 The classical counterpart of the sigma model 

In this section we make a pause from the discussion on the real quantum sigma model 
and investigate its classical counterpart by ignoring the Berry phase Sb- The classical 

3 Nevertheless we will concentrate on U that has the form given in Eq. (4.13). 
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sigma model provides us insight about the way the system behaves at high temperatures, 
where Sb should not be very important. Also, thermal fluctuations may partly replace 
the quantum fluctuations of the real system which, as discussed in Chapter 5, are very 
important. 

The classical model is much more accessible than the real quantum one because its action 
is real, and relying on convexity arguments, we can look for a mean field theory that 
mimics its dynamics. This classical action is also amenable to Monte Carlo simulations, 
but we defer any progress in that direction to future research. 

The partition function of this classical model is given by 

Classical = J Da s «*~*- , (4.41) 

in d = 3, with ^classical = H(a) defined by Eq. (4.29). In particular we concentrate on the 
nearest-neighbor (NN) theory only, and have 

^classical = IT \ ~~q~ ) J2 Tl ^n^n+ft- (4.42) 

Using mean field (MF) theory, we analyze a toy sigma model that has N = 3. We restrict 
to a set of configurations with different values of baryon number density n#, and derive the 
classical MF phase diagram in the temperature-density plane. We will show that the MF 
ground states at zero temperature break the U (3) symmetry, and that as we increase the 
temperature, thermal fluctuations melt the condensate and restore the symmetry above 
a certain temperature T c . We will show that T c decreases with increasing 

4.3.1 Mean field theory 

Our starting point to mean field theory is the partition function Eq. (4.41) that we rewrite 

Z = J Da e~ s °^ x e _ ' Sclassical ( CT - ) ~ s ' o(cr - ) l = Zq • ^ e _ ' s ' classical ^ _ ' So ^' ) . (4.43) 

Here is a trial MF action, and Z is its partition function. Also the brackets (, ) denote 
the MF average of any functional A that depends on the sigma fields. 

( A ({a n }) ) ee Z^ J Da e~ s ^A ({a n }) . (4.44) 
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Using standard convexity arguments one has 

^ g— [^classical — So] ) < g~ ( [^classical - So] } (4 45) 

This means that the free energy F — — log Z has a upper bound, and obeys the Bogoliubov 
inequality, 

F < F + ( S classical -S ) = $. (4.46) 

This bound is useful since we can now guess a MF action So, minimize it with respect to 
its parameters and obtain an evaluation for the real free energy of the system. The value 
of the upper bound $ determines how close the MF guess is to capture the dynamics of 
the system. 

In our case we choose the simplest possible MF action that corresponds to decoupling all 
the nearest-neighbor interactions in ^classical- We first choose to proceed in the basis of 
U (N) generators M v that obey 

Tr M*M* = \ V, £ (M% (M*% = ±5 a 6 jk . (4.47) 

In this basis the sigma field is given by its iV 2 components a v , 

a = 2j2^ V M riT , (4.48) 



and we write the MF action as 



2 

parametrized by the set h a . The function $ is then 



, N \ 2 

Smf = PJi(^) (4.49) 



-£log{/ da n e - K9 ^\ + K £ ■ fin - K £ fin ■ L, (4.50) 

where the first contribution is F , and the next two terms correspond to ( S^ia^icai ) and 
(Sq) q . K = (3J\ (pf} an( i A* is the magnetization 

f da p- KS ^nff 

^ ]_a^e a„ 

^ J da n e~ Kffnhn V ; 
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Note that the integration measure da n depends on the U(N) representation chosen on 
site n. We minimize <3> with respect to {h n } to obtain the MF equations, 

fl n = - - ^, 4.52 

Jda n e- Kabfn ^^)^ 

which are extrema (in jl space) of $. This allows us to write $ as a function of jx only, 
^-^logff^e^^U^ fl n -fl m . (4.53) 

n ^ J > (nm) 



The MF equations Eq. (4.52) and the function <3> in Eq. (4.53) constitute our classical 
MF theory. For each value of K ~ 1/T, $ has a set of extrema points that obey the 
MF equations. A subset of these are minima that correspond to classical ground states 
of the MF theory. The magnetizations, m n = 2/2 n • M T , of each minimum determine the 
symmetry properties of the corresponding ground state. 



4.3.2 The toy model of £7(3) 

The Nf = 1 theory is a {7(4) antiferromagnet and contains many degrees of freedom in 
which to do mean field theory. We treat the simplest non-trivial model which is a toy 
model with {7(3) symmetry. It does not correspond to an actual value of Nf, but allows 
for analytical evaluation of the MF quantities. The most symmetric B = state in this 
model cannot have B = on every site, but must alternate between the representations 
corresponding to B = ±1/2 (see Fig. 4.1). Increasing B is done by replacing some of 



B=-l/2 B=l/2 B=-l/2 B=l/2 B=-l/2 

Figure 4.1: The most symmetric B = state in the Z7 (3) toy model 



the B = —1/2 representations with B = +1/2 representations as in Fig. 4.3.2. Another 
option is to replace either B = —1/2 representations or B = +1/2 representations by 
i7(3) singlets with B = +3/2 and three rows. 
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B=-l/2 B=l/2 B=l/2 B=l/2 B=-l/2 

Figure 4.2: Adding a baryon in the U(3) toy model. 



As mentioned, several MF theory quantities such as $ can be calculated analytically for 
this toy model. The reason is that for N = 3, and m — 2, Eq. (4.13) takes the simple 
form 

( 

exp 
V -at 




We use the definitions 



ai 



pcos(#/2)e* 



a 2 = psin(fl/2)e ia2 , 
with 0<p<7r, 0<#<7T and < ai j2 < 2ir, and get 

sin p sin# cos0 
— sin 2 psin# sin0 



-2 



sin 2 p cos# 
sin2pcos6»/2cos(|(^ + 0)) 
-sin2p cos6»/2 sin(|(^ + 0)) 
sin2psin#/2 cos(|(-0 — 0)) 
-sin2psin#/2sin(!0 - 0)) 

-V3(C08 2 P - |) 



(4.54) 



(4.55) 



(4.56) 



here = «i — a 2 and ip — a± + a 2 , with < ip < An and < < 2n. 

At zero density we assume that the magnetizations on even and odd sites are pi and p 2 
respectively. This results in the following MF equations 



pi = 



/*2 



Jd<7e- 6K °-^a 

J dae +6Kff pl a 
J dae +6Kff h 



(4.57) 
(4.58) 
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The measure of the manifold U(3)/U(2) x U(l) is the measure of the manifold of 3 x 3 
hermitian matrices da modified to account for the dependence (4.56) of a on the manifold 
coordinates. The induced metric of the space is 



induced 

9af3 



d6 a 86 ^ 



(4.59) 



where g^ v = 5^, and the coordinates 9 a are the coset coordinates that appear in ct m . The 
measure of integration is thus 



da = d^det (induced)- 
For U(3)/U(2) x 17(1), we find that 

da = d(sm 4 p) d(cos9) — 

y H> y ' 2vr 4tt 



(4.60) 



(4.61) 



In order to solve the MF equations we choose a basis that diagonalizes m\ = 2 J2 V (J>iM vT . 
Since the measure in the angles \I/ and is flat, m 2 = 2J2 V f I 2^ r ' T * s diagonal as well, 
and only and /x| are non-zero. Since /x° is fixed [see Eq. (4.56)], one is left with 

the following MF equations for fli = (fxf, /if) with % — 1, 2, 



a*1 



/ 3 (6X^,6X^), 
-/ 3 (-6A>?,-6A>!), 



-f s (-6Kfi,-6Ktf). 



(4.62) 
(4.63) 



Here 



h{x,y) 



e x+ e x ~ 

X-\- X— 



- 1 



- 1 



X- 

e x 



x u 



- 1 



X. 



Z r- 

fs(x,y) = ^ + v3 



x+ 



X+ 



x + 



X- 



i e 
- + — 



1) 



X 



+ 



X_ 



i r 



x 



- 1 e x + - 1 



x 



(4.64) 



(4.65) 



with x± = \^3y ± x. 



4.3. The classical counterpart of the sigma model 



75 



We use Eq. (4.62) to evaluate ®(fa) as a function of fa, 



$(fa)/N s = -3K fa- fa- 2KV3(fa - fa) - (p(6Kjt 



<p(xi) = 



log 



-X-A 



- I - log 1 



IJ-I 



e x +' 2 - 1 



(4.66) 



(4.67) 



Here N s is the number of sites and x±j = Vs^f ± /if. It is easy to check that the MF 
equations have a symmetric solution with fa = (and hence fa = 0) for all values of K. 
The magnetizations of this solution are 



mi = — m 2 



1, 



(4.68) 



and are invariant under U(3). As we increase K, three additional nontrivial minima 
appear. In the fa-fa plane, the solutions are found on the following three lines. 



to 



0, 

-Vsfa, 



(4.69) 
(4.70) 
(4.71) 



They have the same $. The magnetizations of each solutions are invariant under a certain 
U(2) x U(l) subgroup of U(3). For example the solution with fa = has 



mi ~ — 777,2 ~ 



/ 1 \ 

1 
V -2 J 



(4.72) 



For K > K c these solutions have lower free energy than the symmetric solution (4.68) 
and £7(3) is spontaneously broken to U(2) x U(l). In Fig. 4.3 we present $ on the line 
fa — for K/K c = 0.98, 1, 1.02. We find that there exists a first order phase transition 
at K r = 0.9. 
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Next we performed the MF analysis for a set of configurations that have non-zero baryon 
number. We construct these by putting net baryon number on the sites of 2 x 2 x 2 unit 
cells. This set is listed in Table 4.1, where we give a pictorial description of the unit cells 
of each configuration. 

Table 4.1: The set of configurations with different baryon 
number density for the U(3) toy model. Here p max = 3/2. 
A ± represents B = ±1/2, and a • represents B = +3/2 
with a — 1. 



Configuration 



P/Pn 



j_ 

12 



_5_ 

12 



Unit cell 



+ i+ 



+ + 



K r . 



0.9 



0.94 



1.03 



1.14 
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Configuration 



Table 4.1: (continued) 
p/Pmax Unit cell 



7_ 

12 



1.72 



1.5 



The numbering of the sites in the unit cells denotes the magnetization distribution of the 
ansatz used for each configuration. For example, in configuration no. 1 the ansatz consists 
of having two magnetizations denoted by m 1 and m 2 on the even and odd sites denoted 
by 1 and 2. We present the MF equations and corresponding $'s for all configurations in 
Appendix D. 

Except for configuration no. 2, the MF equations are consistent with the assumption that 
all magnetizations commute and are diagonalized simultaneously. Therefore as for zero 
baryon number, the MF equations and <3>'s can be written for /i 3 , and fi 8 only. Also, in all 
these configurations $ becomes a function of two variables, allowing a simple graphical 
minimization of $. For configuration no. 2 $ is a function of four variables, and we make 
an ansatz that only /z 8 is non-zero in order to treat this case graphically as well. 

Except for the symmetric solution that exists for all values of K, nontrivial minima of $ 
appear on the lines (4.69)-(4.71). In configurations 3, 5, 6 one has three equivalent minima 
that have the same symmetry and $. As for zero baryon number U(3) is spontaneously 
broken to U(2) x U(l) above a certain value of K, and the phase transition is first order. 
We give the values of K c in Table 4.1. In configurations 4, 7, 8 we find six minima. 
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Four of these correspond to breaking U(3) into U(2) x U(l), and the remaining two 
correspond to breaking U(3) to U(l) x U{1) x U(l) 4 . All six solutions have the same $, 
and become preferable over the symmetric solution above the same value of K c . Here the 
phase transition is second order, and we illustrate it in Fig. 4.4 where we plot $ on the 
line (4.32) for K/K c = 0.8, 1, 1.16 (for configuration 4). 

As mentioned, for configuration no. 2 we take an ansatz that /x 3 = 0, and find similar 
behavior to that of the zero density system. In view of what happens in the other cases, 
we believe that there are at least two more equivalent minima with the same free en- 
ergy and the same symmetry properties. There is also a possibility that there are other 
nonequivalent minima, with different symmetries as in configurations 4, 7, 8. 




Figure 4.4: <3> for non-zero density near the second order phase transition. 



We give the results of these MF calculation in the phase diagram of Fig. 4.5. Increasing 

4 This is the symmetry of the ground state of the quantum sigma model at these configurations, as 
predicted in the large N c limit. See Chapter 5. 
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the baryon number density decreases the transition temperature. Note however that we 
do not find a restoration of the symmetry at zero temperature for any density since this 
is a MF theory which always breaks the symmetry at T = 0, even at one dimension. 



1.25 



0.75 



0.5 



0.25 



O 1 st order PT 
+ 2 nd order PT 




1 1 1 1 1 1 1 

2 4,6 8 10 12 



Figure 4.5: Classical phase diagram of the U(3) antiferromagnet at finite temperature and 
density. We find that when the phase transition is second order, U(3) is spontaneously 
broken either to U{2) x U(l) or U{\) x (7(1) x U(l). In all other cases U(3) is broken to 
U(2) x 1/(1). 



4.4 The semiclassical large N c limit 

We now treat the full quantum sigma model in the large N c limit by performing a saddle 
point evaluation of the ground state of the system. We will see that the ground state 
is degenerate. For B = this degeneracy is discrete, while for B > the ground state 
has a huge continuous degeneracy, exponential in the lattice volume. In general these 
degeneracies are removed only by l/N c corrections that we calculate in the next chapter. 
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4.4.1 Zero density 

The simplest way to set the baryon density to zero is just to choose B n = on each site, 
meaning m = N/2 = 2Nf (see Fig. 4.6). It turns out to be just as easy to consider a 



Figure 4.6: B = configuration, for the case of Nf = 1. 

slightly generalized case [21], in which B n is chosen to alternate, B n = ±b, on even and 
odd sublattices (see Fig .4.7). This means to choose a representation with m = N/2 + b 



Figure 4.7: B = state, for the case of Nf — 1, and b = ±1. 

rows on even sites and N — m = N/2 — b rows on odd sites, which gives a pair of conjugate 
representations of U(N). In view of Eq. (4.18), we can substitute a — > —a on the odd 
sites and thus have identical manifolds on all sites. The Hamiltonian is then 

= ( E ^(-W^A • ( 4 - 73 ) 

m is a new parameter in the theory, and we can ask what value of m gives the lowest 
energy for the ground state. 

In the large- iV c limit, we seek the classical saddle point 5 of S. We assume the saddle is 
at a configuration cr„(r) that is independent of time, and so we drop the time derivative. 
We begin with the nearest-neighbor Hamiltonian, 6 

=_- -y (y) 2 ETr^ +A - (4-74) 

5 Note that the kinetic term is pure imaginary. 

6 This classical Hamiltonian is ferromagnetic; the antiferromagnetic nature of the quantum Hamiltonian 
(4.26) is preserved by the alternating signs in the time-derivative term. 



82 



Chapter 4. Nonlinear sigma model and the classical ground state 



Again, the matrices a n are Hermitian and unitary, and the expansion coefficients a\ are 
real and satisfy XXa'') 2 = N/2. The minimum of H is clearly at a constant field, a n = ao, 
which can be diagonalized to o n = A by a U(N) rotation. This is a Neel state in the 
original variables. The U (N) symmetry is broken to U (m) x U (N — m) and there are 
2m(N — m) Goldstone bosons. 

The classical energy density (per link) is e = —JiN^N/8, independent of m. Thus 
at leading order in 1/N C , the optimal value of m is undetermined, and any alternating 
background of baryon number is equally good. 

Since it turns out that the 1/N C corrections select m = N/2 [see Section 5.1 let us 
consider the effect of the next-nearest-neighbor (NNN) term in H for this case only. The 
perturbation is 

#' = ^(y) 2 |>^a4 ( 4 - 75 ) 

It breaks the U(N) symmetry to SU(N f ) L x SU(N f ) R x U(l) A x U(l) B . Assuming that 
J2 *C Ji, w e again seek the minimum energy configuration in the form of a constant field; 
we minimize 

e 2 = 5>VX ( 4 -76) 

m 

among the U (iV)-equivalent a = a states that minimize the nearest-neighbor action. It 
is not hard to show (see Appendix C) that €2 is minimized for a — 70 <S> 1- This is a 
condensate that is symmetric under the vector generators M v = 1 ® A a but not under the 
axial generators 75 <g> A a , and thus it breaks the chiral symmetry to the vector subgroup, 
SU(N f ) v x U(1) B . 

It is interesting to note that the discrete degeneracy in the number m can be removed at 
0(1/N®) as well by the NNN terms themselves. To see this we calculate the NNN energy 
Eq. (4.76) for the configurations that minimize the NN enegry Eq. (4.74). It is easiest to 
do the calculation by writing the NNN energy in terms of the Dirac matrices as 

e 2 = \ Tr [oa-iiOOL^ . (4.77) 

We show in Chapter 6 [Eq. (6.13)] that Eq. (4.77) has a lower bound 

3 

e 2 > -(4m -37V), (4-78) 
which means that the energy from NN and NNN interactions is bounded by 



^ x 3iV c 2 
E m) > — s- 



-J ± N +^(4m-3AT) 



(4.79) 
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Since m > N/2 then E(m) > E(N/2). In Appendix C we show that the ground state for 
m = N/2 saturates E(N/2) = —3N/2. This means that the ground state of m — N/2, is 
also the overall ground state (out of values of m). As a result we find that the discrete 
degeneracy of the classical NN energy in the number m is removed by NNN as well as by 
1/N C corrections. 



4.4.2 Non-zero density 

The zero-density theories considered in the preceding section were defined by selecting 
representations with m and N — m rows on alternating sites. For any m, this led to a a 
model with identical degrees of freedom on all sites — after redefinition of the spins on the 
odd sublattice — and ferromagnetic couplings. 7 We eventually settled on m = N/2 as the 
background that gives the ground state of lowest energy. 

Introducing non-zero baryon density means changing m on some sites of the lattice. 
Choosing a different representation adds (or subtracts) baryons on a site-by-site basis. 
For instance, we can add a single baryon by adding a row to the Young tableau as in 
Fig. 4.8. Since in general there will be representations on different sites that are not 



B=0 B=0 B=l B=() B=0 

Figure 4.8: Adding one baryon to a site in the Nf — 1 theory [a U(4) antiferromagnet] 

mutually conjugate, different sites will carry a variables that do not live in the same 
submanifold of U(N). We here limit ourselves to the simpler case of uniform m, where 
adjacent sites carry identical spins — but the coupling is antiferromagnetic. 

In order to learn how to work with such a theory, we begin by studying the two-site 
problem. The results of this study will lead directly to an ansatz for the ground state of 
a lattice with a fixed density of baryons. 

7 The (—1)" factor in the kinetic energy retained information about the antiferromagnetic nature of 
the quantum problem; it did not affect the classical analysis. 
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t 



m 



m 



Figure 4.9: Two sites that carry U(N) representations with mi and m 2 rows. 



Two site problem 

Consider, therefore, two sites with quantum spins Q 1 and Q 2 that carry representations 
of U(N) with mi and m 2 rows, and N c columns (see Fig. 4.9). 

The quantum Hamiltonian is 

H = J 1 J2Q r lQl (4-80) 

v 

an antiferromagnetic coupling. The corresponding classical a model has the interaction 
Hamiltonian 

J / AT \ 2 

H ^ = i{Y) Tr ^, (4.81) 

where 

a t = UiAiUl (4.82) 
The two A matrices reflect the different values of m 8 according to 




(4.83) 



The N c — > oo limit is the classical limit, in which we seek values of <7i ; 2 that minimize 
H(a). A global U(N) rotation, viz., 

a % -> VaiV\ (4.84) 

can be used to diagonalize Oi so that a 1 = A 1 . Now we have to minimize Tr A ± a 2 - The 
case of conjugate representations, m 2 — N — m±, is easy: a 2 is a unitary rotation of A 2 , 
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which (in this case) can be rotated into — Ai. This is the unique antiferromagnetic ground 
state. o\ and o 2 can be copied to the odd and even sublattices of an infinite lattice to 
give the classical Neel state considered in the preceding section. 

The case rrt\ = m 2 = m is more complex. We consider m > N/2 for definiteness; the 
other case is similarly handled. We write explicitly 



02 



\ 



sin 



sin 



—a 



aa> 



(2vW) \ 

V at a 

— cos (2\/ at a 



(4.85) 



Since a 2 is a rotation of A, 

2m- N = Tkt 2 = Tt cos (2v / aat) - Tr cos (2v / a^) , 
and hence the energy is 



E 



2 



Tr cos (2v / ota) +m- N/2 



(4.86) 

(4.87) 
(4.88) 



E is minimized when all the eigenvalues of a) a are equal to 7r 2 /4. This means that the 
N — m column vectors a, form an orthogonal set in m dimensions, with 



a] ■ a,- = I r I 5. 



(4.89) 



Since m > N — m by assumption, such a set of vectors can always be found. 

Since a)a = (it 2 /4)ljv- m , we have sin (2V a+a) = and so the off-diagonal blocks of 
Eq. (4.85) vanish. The lower-right block of 02 is the unit matrix l^-m- We know that 
(C2) 2 — 1 since A 2 = 1 and thus the upper-left block must have eigenvalues ±1. Equating 



traces of a 2 and A, we find that the upper-left m x m block must take the form 



CT (m) _ Jj(rn) ]^(m)jj(m)\ ^ 



where 



1-2m-N 






-ljy-m 



(4.90) 
(4.91) 
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and U (m) E U(m). represents the coset U(m)/[U(2m - N) x U(N - to)]. 

We conclude that the classical ground state of this B ^ two-site problem is degenerate, 
even beyond breaking the overall U (N) symmetry. The solutions can be written as 




(4.92) 



(to which a global U(N) rotation can be applied). A particular instance of is A^ m \ 
given by Eq. (4.91). The symmetry of the vacuum is the set of rotations that leaves both 
<j\ and <t 2 invariant, namely, U{2m — N) x U(N — to) x U(N — to). 



Infinite lattice 

At N c = oo we seek the saddle point of the action, which we assume to be a time- 
independent configuration. The classical Hamiltonian of the a model is 

#=y(^) 2 ETr^ +A . (4.93) 



n,fi 



Seeking an antiferromagnetic ground state, we set a n = A on the sublattice of even sites. 
The odd sites are then governed by 

#odd = Jid f c £ Tr Aa n . (4.94) 

V 1 ' n odd 

This is just the two-site problem studied above, replicated over the lattice. As we saw 
above, the ground state energy is degenerate with respect to the configuration at each 
odd site, 

/ rr( m ) D \ 

( o iL)- (495) 

and we find that at N c = oo there is a large number of ground state configurations 
with the same classical energy, that differ in the choice each of the odd sites. 

These ground states are continuously connected by the coordinates of that furnish 
U (m)/[U (2m — N) x U (N — to)]. They have different realizations of the global symmetry 
of U(N); A uniform choice for the odd sites, = A*™) for instance, breaks the U(N) 
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symmetry to U(2m — N) x U(N — to) x U(N — to); a non-uniform choice can break 
the symmetry all the way to U(N — m). Furthermore, one can contemplate making a 
non-uniform ansatz for the even spins as well. The entropy of this classical ground state 
is evidently proportional to the volume. 

The situation is reminiscent of that in the antiferromagnetic Potts model [79] . The crucial 
difference is, however, that here the degeneracy is continuous, and the sigma fields are 
allowed to oscillate around the classical ground state, even at zero temperature. In fact 
this makes the problem more similar to problems encountered in condensed matter physics 
in the context of frustrated magnetic systems. These systems are known to experience the 
phenomenon of "order from disorder", [52, 53, 54, 55, 56, 57] where fluctuations remove 
a classical degeneracy and choose a ground state. We discuss how order from disorder 
happens in our problem in Chapter 5. There we show that when N c = oo, and one 
keeps a n on the even sites fixed to be A, it costs no energy to change a n on an odd 
site by changing its direction within U(m)/[U(2m — N) x U(N — to)]. Fluctuations in 
these directions are "zero modes" and have zero energy for any momentum. In principle 
they can lead to instabilities of the classical ground state, but as we will show, order 
from disorder gives the zero modes non-zero energy of 0(1/N C ), removes the classical 
degeneracy, and stabilizes a canted ground state in which all the even sites point to A, 
and all the odd sites align uniformly to have the same a^ m \ 

A final note here is that different from the zero baryon case, the NNN interactions do not 
remove the degeneracy of the NN theory. In particular the zero modes are not removed 
when one takes a non-zero value for J 2 . This can be seen in Section 6.3 that gives the 
dispersion relations of excitations in the sigma model with NN and NNN interaction. 
Taking N c = 00 and keeping J 2 7^ still results in zero modes. 



Chapter 5 



Quantum fluctuations in the 
nonlinear sigma model 



In this chapter we investigate the effects that quantum fluctuations of 0(1/N C ) have on 
the classical ground state of the NLSM introduced in Chapter 4. We restrict ourselves to 
nearest-neighbor interactions only and show that the degeneracies present in the classical 
N c = oo limit are removed at this order. We first show that for zero baryon density 
fluctuations remove the discrete degeneracy in the parameter m [See Section 4.4.1] and 
choose the configuration that has zero baryon number on each site. We proceed to perform 
a similar calculation for non-zero density and see that in that case fluctuations remove the 
infinite continuous degeneracy present in the classical ground state. This degeneracy is 
represented by the fact that in the classical level, the sigma fields on the odd sites can move 
freely independent of each other in certain flat directions. Fluctuations in these directions 
have zero energy at 0(1) for all momentum and are zero modes. We find that quantum 
fluctuations give these modes non-zero energy of 0(1/N C ), and stabilized the classical 
ground state. Furthermore we find that these energies are quadratic in momentum, which 
makes them type II Goldstone bosons. 



5.1 Removal of classical degeneracy at zero density 



At zero density we consider fluctuations around the a n = A minimum of S. 



2 Jo 



- E("l)"Tr AUld T U n - y E Tr <^»+£ 



(5.1) 
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with J[ = 2Ji/N c ~ l/g 2 N c . (From here on we drop the prime) Recalling Eq. (4.13), we 
write 



U r . 



(5.2) 



where A n is anti-Hermitian and anticommutes with A. It is more convenient to work with 
the Hermitian matrix 



L n — 2A n A, 



in terms of which we expand 



a n = A + L n - -L 2 A. 
If we further expand L n in the basis of generators of U(N), 



(5.3) 



(5.4) 



(5.5) 



we find that the l v corresponding to generators of U (m) x U(N — m) vanish; this is the 
subgroup under which the vacuum is symmetric. The field L n thus contains 2m(N — m) 
real degrees of freedom, corresponding to the Goldstone bosons. 



We expand U n and a n in powers of L n \ using Eq. (5.5), we obtain to second order 

/ dr^ ——(;"" - - 



r 



' 1 \n j 



The coefficient matrix is 



C w = Tr (A[M V , M v ']), 



(5.6) 



(5.7) 



and the classical energy is = — ^ (^) NN s d{3. C is antisymmetric and purely imagi- 
nary; we show in Appendix E that C has eigenvalues ±1, each with degeneracy m(N — m). 
We change basis so as to diagonalize C, and write the index r\ as the compound (a, ±) 
with a — 1, . . . ,m(N — m) and the ± corresponding to the eigenvalue of C. Since the 
original l v are real, we have 

(r+)* = r-. (5.8) 



Thus we eliminate l a and write 



(-lrnmc^c + ^EtoA-c 2 



(5.9) 
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The alternating sign in Eq. (5.9) is what makes the theory antiferromagnetic. It forces 
us to differentiate between even and odd sites, and we transform to momentum space as 
follows (dropping the + superscript): 



n V n. v 1 ^v* 



n even, 
ln - i] n odd. 



(5.10) 



The even sites comprise an fee lattice with lattice constant 2, and the momenta k take 
values in its Brillouin zone. We obtain 



Here 



M(k) 



' -d T + J x d - Jid7(fe) N 
v -Jid-f(k) +d T + Jid J 



(5.11) 



(5.12) 



and 7(fe) = ^E^cos/^. 



The Gaussian path integral over the / field is very similar to the path integrations that 
appear in Chapter 3 and gives the free energy, 



F = F + m(N - m)— - ^ 



(5.13) 



where u(k) = dJiN cy fl - 7 2 (fe) and F = -\ J 1 dNN s (N c /2) 2 . For the ground state 
energy, we take j3 — > oo and obtain (restoring all constants) 



E = -J ± N s Nd 



1 + 



1 m(N — m) 



N 



(5.14) 



This is exactly the result of Smit [21]. The 0(1/N C ) corrections lift the degeneracy 
of the ground states with different values of m. The integral in Eq. (5.14) is positive 
and its coefficient contains the number of Goldstone bosons. Thus the state of lowest 
energy is that with m = N/2 = 2Nf, and the symmetry breaking scheme is U(4Nf) — > 
U(2Nf) x U(2Nf). Further breaking by the next-nearest-neighbor terms was discussed in 
Chapter 4 where we extended this result to the next-nearest-neighbor theory and found 
that its U(Nf) x U(Nf) chiral symmetry is broken to the vector U(Nf) flavor subgroup. 
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5.2 Removal of classical degeneracy at non-zero den- 
sity 



Recall the classical solution of the nonzero-density two-site problem, given in Chapter 4. 
We found that the sigma fields are 



and a m is given in Eq. (4.90). The symmetry of the vacuum is the set of rotations that 
leaves both o~\ and o~2 invariant, namely, U (2m — N) x U (N — to) x U (N — to). There are 
two kinds of degeneracy in this two-site problem: that which results from breaking the 
global U (N) to U (to) x U(N — to), and that which comes of breaking the U (to) subgroup 
to U (2m — N) xU(N — to). The latter degeneracy is connected with freedom in choosing 
the orientation of a 2 relative to a\. It is instructive to see how quantum mechanical 
fluctuations lift the classical degeneracies. 

The quantum two-site problem is easy to solve. We rewrite the Hamiltonian (4.80) as 



Q\ and Q\ are constants, the quadratic Casimir operator in the TO-row, iV c -column rep- 
resentation of U(N). The first term in Eq. (5.16) is minimized by coupling Qi and Q 2 to 
the representation that minimizes the Casimir, which is the representation with 2m — N 
rows and N c columns (see Fig. 5.1). The ground state has discrete degeneracy equal to 
the dimension of this representation. 

The exact quantum solution naturally shows no sign of spontaneous symmetry breaking 
and hence it is not of much relevance to the infinite volume problem. More interesting is 
the problem where the state of Q\ is fixed and Q2 is allowed to vary. This breaks by hand 
the global U(N) while allowing quantum fluctuations to lift any remaining degeneracy 
in the relative orientation of the two spins, so it can be regarded as quantization in the 
presence of spontaneous symmetry breaking. In effect, this is mean field theory. 

We replace the Hamiltonian (4.80) by 




(5.15) 



h = y [(Qi + Q2) 2 -Ql-Ql 



(5.16) 



H 



MF 



(5.17) 
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X 




Figure 5.1: Coupling two spins in a 7-row representation of U(12) to the representation 
with minimal Casimir operator 



To minimize the energy we maximize the mean field by choosing the state of Q\ to be the 
highest-weight state. This state diagonalizes the generators H a of the Cartan subalgebra 
while other generators of U (N) have expectation value zero. Thus 

N 

H MF — Ji ^2(H la )H 2a . (5.18) 

Q=l 

The operators H a all commute, and their eigenvalues make up the weight diagram of 
the representation. 1 H MF is a dot product of the weight vectors of the two spins. The 
energy is minimized by choosing for Q 2 a state that lies opposite the highest-weight state 
in the weight diagram. As shown in the example of Fig. 5.2, this still leaves a degeneracy, 
albeit a discrete one. We stress that this degeneracy comes from freedom in the relative 
orientation of Q\ and Q 2 - 

In the N c — > oo limit, the discrete degeneracy becomes infinite and presumably it is well 
described by the continuous degeneracy of the classical problem. As mentioned this is 

1 More precisely, the weight diagram shows eigenvalues of the N — 1 traceless diagonal generators of 
SU(N). These can be obtained by isolating the U(l) member of the set H a and taking linear combinations 
of the rest. 
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Figure 5.2: Weight diagram for N — 3, m — 1, N c = 2 [the sextet of 577(3)]. The highest- 
weight state lies at the top of the triangle; all states at the base minimize H MF . All 
values of A" c give the same triangular shape, with multiplicity of one along the boundary. 
There are iV c + 1 states at the base. For N = 4 the triangle becomes a tetrahedron, and 
there are (A" c + l)(A r c + 2)/2 states at its base. 



similar to problems of frustrated magnetic systems, where a local group of spins will have 
a degenerate ground state, with a continuous degeneracy. 

This degeneracy, however, is an artifact of the mean- field approach that, like the classical 
ansatz, assumes a fixed state for the spins on the even sites. In fact the sigma fields will 
fluctuate as soon as they are allowed to do so, i.e., when N c < oo. A given odd spin 
will not be surrounded by a uniform fixed background; the neighboring even spins will 
be influenced by all their odd neighbors, and will induce an interaction among the odd 
spins. This should reduce the entropy of the ground state to zero. 

The systematic way to see this effect is to carry out an expansion in 1/N C and verify the 
stability of the classical ansatz. We do that for the canted ground state, which is the 
solution to the two-site problem, replicated uniformly over the lattice: 



(5.19) 



where = U^A^U^ m ^ is now a global degree of freedom. 

The meaning of the classical local degeneracy is that at 0(1) there is a set of directions 
in which one can move from one classical vacuum to another without cost of energy. 
Fluctuations in these directions are zero modes, that have zero energy for all momentum, 
and we discuss them in Section 5.2.1. The sign for stability must therefore be the removal 
of these zero modes at 0(1/ N c ), which we demonstrate in Section 5.2.2 along the lines of 
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the phenomenon of "order from disorder" [52, 53, 54, 55, 56, 57] in frustrated magnetic 
systems. 

The essence of order from disorder is that fluctuations of quantum, thermal or even 
quenched nature can lift a classical degeneracy. A famous example is the Kagome antifer- 
romagnet [55] whose classical ground state energy is invariant under correlated rotations 
of local groups of spins, leading to a degeneracy exponential in the volume. This system 
also has modes with zero energy for all momentum. These zero modes obtain non-zero 
energy due to quantum fluctuations. Most discussions of order from disorder work in a 
Hamiltonian framework (but see [56, 57]). We employ instead the Euclidean path in- 
tegral. In Appendix H we also present a Hamiltonian calculation using a generalized 
Holstein-Primakoff transformation. The latter can only be constructed easily for a subset 
of the cases we consider [21, 49]. The Euclidean calculations are less restrictive in their 
applicability. 

Finally, as an exercise, we present a derivation of the effective ferromagnetic interaction 
among the a n on the odd sites in Appendix I, which in fact motivated us to investigate 
the canted ground state in particular. 



We begin to investigate the stability of the canted ansatz, where the vacuum expectation 
values (vevs) of the a fields on the odd sites point in the same direction (a ferromagnetic 
alignment). Without loss of generality the ground state can be chosen to have 



5.2.1 Zero modes at O(l) 



(Oeven) OC A, 




(5.20) 



with 



V = 



( 1 \ 
1 



(5.21) 



\o -i o/ 



These vevs break U(N) to U(2m - N) x U(N - m) x U(N - m), with 
tlq = 2 (3m — N)(N — m) broken generators. Here we write an effective action for the 
fluctuations around the up to 0(1/N C ). Working at 0(1) first, we identify the zero modes 
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that correspond to the classical degeneracy. We then show how quantum fluctuations give 
them non-zero energy. 

In this and next sections we choose to work with the parametrization of the coset space 
U(N)/[U(m) xU(N-m)] [49] ofEqs. (4.21)-(4.22), in terms of the m x (N-m) complex 
matrix 0. As mentioned, the main advantage of this parameterization is that the kinetic 
part of the action is quadratic in [see (4.40)]. 

The fields 0, suitably shifted, can be identified with Goldstone bosons around the ground 
state (5.19). For the even sites, = indeed gives a = A. For the odd sites, the vacuum 
is at 

*=d° )• 

\ i-N-m J 

where the upper part of is a (2m — N) x (N — m) matrix. We therefore shift on the 
odd sites according to 



U' = VU, 
cr'(0) = Va(<j>)V*. 



(5.22) 



This gives cr'(0) = A2. We drop the primes henceforth. For later convenience we write 

as 

>-(*.)■ 



(5.23) 



Here 7r is an (A — m)— dimensional square matrix while X has (2m — N) rows and (N — m) 
columns. Both are complex. Thus 



and 



O"odd 



V 



1 - 2xX f 

-27T X t 
-2yjl - 0t0 X t 

1 - 2XX 1 

-2y/l ~ 0t0 X t 
27T X t 



-2 X 7rt 

1 - 27^ 
-2Jl- 0t07T t 



27T./1 



-2 X y/l - 0+0 ^ 
-2nyjl - 0t0 
-1 + 20+0 J 

2 X vrt N 

1 + 20t0 2y/l - 0t07T t 
1 - 27T7T t 



(5.24) 



-2Xa/1 



(5.25) 



Each submatrix in Eqs. (5.24)-(5.25) has dimensions as indicated in (5.19). 
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After scaling — > 4>/\fW c the action takes the form 



nfj, 



Jl 



[XriXfi+A 71 "'"- Xn+flXn^n+p, + ^-C 



(5.26) 



n/j. 



. ^i r t t _ t t _ t t 



N, 



nfj. 



^n-^nXn+fiXn+p, ^ n+ fj^ n + fiXnXn 



+ 



n+fj^n+fj^n + ^ n+ jj^ ■ 



1 f t 

^n+fjXn+fjXn+fi^n + 7T Tl+A 7r Ti xlaXr l + h.C?j 



+0(l/iV c 3/2 ). 



The AF structure of the action demands the introduction of an fee lattice. We write the 
fields 4> n = 4>n w ith A = (even, odd) and A?" belonging to an fee lattice. We Fourier 
transform according to 



AT 



T 



— /-V 



LA i(k-N-wr) 



X 



1 A = even 

e -ik z /2 A = odd ' 



(5.27) 



Here k = (fc,o;). In momentum space the action is given by S = S2 + 5*3 + S4 with 
(discarding Umklapp terms) 



S2 = TrE(< f ^ 



— £<j + 2dJi —2Jid^k 



Til. 



-2 Jid^k iuj + 2 Jid / I 7r 



of 



5 3 
S 4 



(^) T )(T ° )((!))• 

2 WiV?5 Tr ^ (x?Xk<-k7fe - X? Xl-p^plk-p + /i.e.) , 

V s c < kp 



(5.28) 
(5.29) 



-'" /Ui/ 'r r ~~Tr ^2 (xlXp XqXq+k-p ^k^p^q^q-k+p 71 "fc ^ 7r p 7r g 71 'q-k+p 

-XkXl^q+k-p ~ KrfXqXq-k+p) lk-p (5.30) 
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( of o i of o\ ( C O i of ct \\ 

(Tfc + Xfc Xp) [^k-p-q^q + ^-'^-(fc-p-,) J J 7fc-p-<r 

7 fc is given by 

1 d 

7fc = iE cos y 2 ' ( 5 - 31 ) 
At large iV c , S 3 and S4 are small perturbations. The bare propagators can be read from 

(r m\ _ 1 ( iu + 2Jid 2J x d^ k \ 

[ k)AB ~ u* + 4J?<PEl\ W lk -iu + 1hd) A ^ [b - 6Z) 

(The propagators are diagonal in the internal group indices.) Here E k = ^1 — 7^. The 
poles of the propagators give the dispersion relations of the various bosons at 0(1). The 
result is 

uj 2 + {2J l dE k f = for vr, (5.34) 

and 

cu = for x- (5.35) 

There are 4(N — m) 2 real fields per momentum of the 7r kind on the even and odd sites, 
and the propagator (5.32) has 4(iV — m) 2 poles that obey (5.34). We pair these into 
2(N — m) 2 physical massless excitations with 

iu = ±2J!dE k , (5.36) 

which are 2(N — m) 2 spin waves. It is easy to verify that the Euclidean dispersion relation 
is uj 2 ~ — \k\ 2 at low momentum. We note that at B — they are the only excitations. 

Let us pause to see how the above discussion would change if the interaction between the 
7r fields were ferromagnetic. First the propagator (5.32) would be replaced by 

G^ = iw-2Jid(l- 7fc ), (5.37) 

where k is now a momentum in the first BZ of a simple cubic lattice (the fee structure is 
redundant for a ferromagnet). Here one has 2(iV — m) 2 fields per momentum, all of them 
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are massless and obey iuj ~ J\\k\ 2 . Pairing negative and positive frequencies together, one 
has (N — m) 2 massless excitations only. To see the connection to the antiferromagnetic 
case write the ferromagnetic action and propagator in terms of fee degrees of freedom as 

Sfm ~ E(4 e <°)Gvl(^\), (5-38) 



! ( iuj - 2J x d 2J 1 d'y k 

™ 1 2J x d lk iu-2J x d 



(5.39) 



Here the zeros obey (icu — 2J\d) 2 = (2Jid^k) 2 , or 

iw = 2Jid(l±7fc). (5.40) 

So one has again 4(iV — m) 2 fields per momentum, but half of them have a mass equal 
to 4J 1 d. Next we pair positive and negative frequencies together and find that half of 
the remaining 2(N — m) 2 massless real fields describe (N — m) 2 massless particles. The 
difference between the two cases is that while the number of broken real fields is the 
same ng = 2(N — m) 2 , the number of physical excitations is different. It is uq for the 
antiferromagnet, and nc/2 for the ferromagnet. 

The x fields are 4(iV — m)(2m — N) zero modes ("soft" modes). Their energy is zero 
for all momentum, a sign of the local degeneracy of the classical ground state discussed 
above. 

We conclude this subsection by classifying the different fields according to their U (2m — 
N) x U(N — m) x U(N — m) representations. We denote a representation by 

(ri,r 2 ,r 3 ) (91 ' 92 ' 93) . (5.41) 

Here r\ denotes the representation of SU {2m — N), and r 2i3 denote the representations of 
the two SU (N — m) subgroups, qi are the charges of the excitations under the remaining 
U(l) factors of the unbroken subgroup. These are generated by the following diagonal 
matrices. 

/ l 2m . N \ / \ 

1 N , , and l N _ m . (5.42) 

V / V -ljv_ m J 

We give the different representations in Table 5.1. The crucial point is that the AF spin 
waves and the zero modes reside in completely different representations. This means that 
the separation between these excitations in the 0(1) spectrum will survive at higher orders 
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Field 


Representation 


Dimension 


7T 


(l,N -m, N -m / 


(0,U,+2) 


(N - mf 


Xcvcn 


[2m - N, 1, N - mj 


(0,+l,+l) 


(2m - N)(N - m) 


Xodd 


(2m -N,N-m,lj 


(0,+1,-lJ 


(2m - N)(N - m) 



Table 5.1: U(2m - N) x U(N - m) x U(N - m) representations of the excitations. The 
conjugate fields belong to the conjugate representations. n(n) stands for the fundamental 
(its conjugate) representation of SU(n) and 1 stands for a singlet. 



and that mixing can not occur (unless there is further spontaneous symmetry breakdown 
in higher orders). 

To be precise we give the transformation laws of each of the fields with respect to a 
U (2m — N) x U(N - m) x U(N - m) element given by 

/ ui \ 
U=l u 2 . (5.43) 
V u 3 J 

Here -u 2 ,3 € U(N — m), and U\ G U(2m — N). The fields 7r e and transform as 

7T £ -> M 2 7T e 4, 7T t0 -> M 2 7T t0 4, (5.44) 

while the x fields transform as 

X e -> WirtL x° -»■ ( 5 - 45 ) 

From Eq. (5.26) we see that the action is invariant under these transformations. 



5.2.2 Removal of zero modes at 0(1/N C ) — the self energy calcu- 
lation 

In this section we calculate the self energy of the x fields to first order in 1/N C . We shall 
see that the poles in the soft modes propagators move away from zero energy in this 
order. The Feynman rules for this problem are presented in Appendix F. We treat S3 
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Figure 5.3: Schwinger-Dyson equation for the self energy of the x fields in one loop. 
Solid lines correspond to the 7r propagators, and double dashed lines represent the full x 
propagator. Single dashed lines are amputated external legs. 



and S4 as perturbations to S 2 and present the contributions to the self energy in one loop 
in Fig. 5.3. 

Since the bare x propagator is divergent at uj = for all k, we insert its self energy 
self-consistently, making the replacement 

(Gir 1 = (G&)- 1 - E*. (5.46) 

G is the bare propagator given in Eq. (5.33) and E is a matrix coupling odd and even 
degrees of freedom according 2 to 

E =(| e;)' <" 7 > 

Thus the full x propagator is 

nx = \ [ iuj ~^i,k E 2 ,fe \ , ,„s 

k (-iu - ^(iu - EJ >fc ) - |S 2 ,,f v S* ;fe -iu - E lifc ) ■ { ■ ) 

Here we assume that both matrices and S 2i a; are diagonal in group space and that 

ReS 1>fc <0 (5.49) 

for the stability of the path integral. 

At this point we note that the ansatz (5.47) is invariant under (5.45) only for S 2 = 0. A 
non-zero S 2 leads to the following non-zero vev 

(X e (k) m x°(-k) pj ) ~ A 2 (k)S aP S ij , (5.50) 

which is of 0(1/N C ). This vev is invariant under Eq. (5.45) only for u\ = uj and w 2 = u^. 
This in turn means that when the Schwinger-Dyson equation has a solution with S 2 7^ 0, 

2 We assume that the x propagator has a similar structure as Eq. (5.32). 



102 



Chapter 5. Quantum fluctuations in the sigma model 



then the stationary subgroup U{2m — N) x U(N — m) x U(N — m) of the vevs in (5.19) 
is broken by the 0(1/N C ) correlation function (5.50) to 

0(2m-N) x U(N-m). (5.51) 

A different ansatz can lead to a different symmetry. For example, a self-energy structure 
that couples to x° leaves the U{2m — N) subgroup intact. In principle one should 
compare the free energy of these two ground states. We do not investigate in this direction 
any further since the only difference will be in the symmetry structure of condensates of 
0(1/N C ). The symmetry structure at 0(1) will be the same. 

In fact because the additional condensates here are of 0(l/N c ), the physics related to 
it (Goldstone modes corresponding to the generators it breaks) is weakly coupled to the 
0(1) physics related to the vevs in (5.19), and we ignore it in this work. 

In Appendix G we derive self-consistent equations for E 1; ^ and S 2j fc. Defining tanh#fc = 
—Re S 2i (o,fe)/Re £i,(o,fc), we write these equations in the form of a single integral equation, 

/ (p) I 2 (q,k)smh9 q 
(N - m) tanh 9 k = BZ V J ' = ™ . (5.52) 

r)(k) and Ii,i{q, k) are defined in appendix G. The poles of the propagators may then be 
obtained from Eq. (5.48) [see Eq. (G.15)]. The dispersion relations turn out to be 

2J l( i(N - m) 
2J 1 d{N - m) 



c^fcV-'- - tanh 9k (5.53) 

Cfc- 



The solution of Eq. (5.52) for d = 3 can be obtained numerically by assuming that the 
function 9 k depends only on |fe|. We plot the solution in Fig. 5.4 for iV — m = 1, 2, 3, 4, 5 
(for details see Appendix G). Recall that the baryon density is given by B = m — N/2. 
Thus for Nf < 3 we cover all baryon density short of saturation. 

The form of e k is shown in Fig. 5.5. It is easy to check that the bosons are massless, since 
r)(0) = 11,2(9,0) = implies that a k and (3k vanish at k = 0. This is a direct result of 
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0. 



CD 

jz 0. 

CO 



0. 




|k| 

Figure 5.4: Solution of Eq. (5.52) for N-m = 1,2,3,4,5. The value of tanh 9 for |fc| ^0 
drops with increasing (N — m). 

the Ward identities concerning the global U (N) symmetry. Moreover, near k = both a 
and (3 have a quadratic dependence on \k\. This means that at low momenta e k ~ \k\ 2 , 
characteristic of ferromagnetic magnons. A possible exception is the case N — m = 1. 
There we see that at low momenta, tanh# fc — > 1 and thus it is possible that e k ~ |fc| 2+p 
with p a positive integer 3 . 

Finally we recall that the relation between the number of \ fields and the number of 
physical excitations depends on the dispersion relation. Since these soft modes obey 
e fc ~ |fc| 2 j the 4(iV — m)(2m — N) fields describe only half that many ferromagnetic 
magnons. The case of iV — m = 1 might be different depending on whether p is odd or 
even. In order to determine p it is necessary to improve the resolution of the momentum 



3 £fc cannot be nonanalytic at k = since the integrals in Eq. (5.52) are regular there. 
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mesh. (This may explain the different behavior of tanh6> in Fig. (5.4) and of the resulting 
efc for this case) 



0.08 




2 4 6 

|k| 

Figure 5.5: Rescaled energy of the x bosons for for iV — m = 1,2,3,4,5. The value of 
(N — m) x efc at the boundary |fe| = of the Brillouin zone increases with N — m. 

The fields become two kinds of Goldstone boson. The ir fields are antiferromagnetic 
spin waves with a linear dispersion relation. The x fields are ferromagnetic magnons with 
a quadratic dispersion relation. This is consistent with the loss of Lorentz invariance due 
to finite density. This set of excitations falls into the pattern described by Chadha and 
Nielsen [80] and by Leutwyler [81] in their studies of nonrelativistic field theories. 

In general, a nonrelativistic system that undergoes spontaneous symmetry breakdown can 
posses two types of Goldstone boson. The energy of type I bosons is an odd power of their 
momentum. Their effective field theory has a second time derivative which means that 
each excitation is described by one real field as in the case of a relativistic scalar field. Type 
II bosons have a dispersion relation that contains an even power of the momentum. They 
can appear only in nonrelativistic theories having properties similar to a Schrodinger field 
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theory, namely theories whose action has only a first time derivative. As in the Schrodinger 
case, each excitation is described by a complex field, or two real fields. The number of the 
<p Goldstone fields is of course equal to the number no of generators broken by Eq. (5.51). 
The counting of these massless excitations is summarized by the Chadha-Nielsen counting 
rule, 

rij + 2n n > no, (5.54) 

where rijinjj) are the number of type I and II Goldstone bosons. A well-known example 
is the SU{2) spin system with a collinear ground state, where no = 2. The antiferro- 
magnet has two spin waves with a linear dispersion relation (ni = 2, riu = 0) while the 
ferromagnet has one magnon with a quadratic dispersion relation (nj = 0, nu = 1). As 
mentioned above in our case both nj and njj are non-zero. 

In addition to these massless fields there exist another set of massless fields. These are the 
Goldstone bosons that correspond to the breakdown of U (2m — N) to 0{2m — N), and of 
U(N-m)x U(N-m) to U(N-m), by the 0(l/N c ) condensate (5.50). In this case the 
Goldstone bosons are not represented by any of the fundamental x or n fields. Similar 
to the BCS case, they are collective excitations that correspond to poles in two-particle 
propagators of the generic form (x* e X*°X e X°)- These poles describe zero mass particles 
which are bound states of two x particles. We do not pursue their properties (symmetry 
structure, dispersion relations, etc.) in this work since as mentioned they are weakly 
coupled to the excitations we treat. 

The energy scale of the type II Goldstone bosons x is smaller by a factor of 1/N C compared 
to that of the type I bosons n. This points to a possible hierarchy of phase transitions at 
finite temperature which can be described by a classical model that has the Hamiltonian 

H = J i Tr [Wn+ji] ~ -h Tr J2 [°™o™+2a] • ( 5 - 55 ) 

nfj. np. 

For J 2 = 0, the ground state will be highly degenerate as above. A small positive value of 
J 2 removes this degeneracy and picks out a ferromagnetic alignment of the next-nearest- 
neighbor a fields. Thus for T <ti Ji, J 2 the ground state is invariant only under 4 

U(N — m) x U(N — m) x U(2m — N). (5.56) 

This symmetry pattern will persists until some finite temperature T C FM ~ J 2 . For T > T™ 
the ferromagnetic magnons melt the magnetization and restore some of the symmetry. 

4 We note here that order from disorder is not necessary here, and therefore there is no condensate 
analogous to (5.50). This means that the physics described by this classical model is different than the 
real model. 
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This situation will persists until a second temperature T c ~ J 1 where the symmetry 
will be restored completely. As long as J\ ^> J2 (which corresponds to N c 1 in our 
system), the hierarchy of phase transitions is well defined. 

Finally we mention that recent work [82] on effective field theories for dense QCD also 
predicts the existence of type II Goldstone bosons. It is tempting to identify these with our 
ferromagnetic magnons. First, however, one must reduce the artificial U{ANf) symmetry 
of the action to the physical chiral symmetry. This important issue is discussed in the 
next chapter. 



Chapter 6 



Next-nearest-neighbor interactions 
in the nonlinear sigma model 



The previous chapters are all mainly restricted to nearest-neighbor (NN) interactions only, 
which were symmetric under the global group of U (4iV/). In this chapter we add the sigma 
model next-nearest-neighbor (NNN) interactions that reduce this symmetry to the chiral 
symmetry U(Nf) x U(Nf), thus completing our picture of the lattice theory that has the 
symmetry of continuum QCD. Taking the NNN coupling to be much smaller than the 
NN one, we perturb the NN ground state and its excitations, found in Chapters 4 and 5 
with these NNN terms. These interactions reduce some of the global degeneracy of the 
NN vacuum, and leave only a part of chiral symmetry intact. We find that the global 
symmetry of the ground state varies with Nf and with the background baryon density. 
In all cases the condensate breaks the discrete rotational symmetry of the lattice as well 
as part of the chiral symmetry group. We also find two types of massless excitations: 
type I Goldstone bosons with linear dispersion relations and type II Goldstone bosons 
with quadratic dispersion relations. Some of the type I bosons exist also at zero density. 
The others were originally type II bosons that received energy contributions, linear in 
momentum, from the NNN terms, and became type I bosons as well. Both kinds of 
bosons can have anisotropic dispersion relations. 
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6.1 Action and symmetries of next-nearest-neighbor 
interactions 

The NN effective Hamiltonian and sigma model discussed up to now in Chapters 3-5 cor- 
respond to an underlying kernel of naive fermions. Although these fermions suffer from 
species doubling in weak coupling, this doubling is not seen directly in strong coupling. 
The reason is that their we do not have fermion excitations, but rather only gauge in- 
variant degrees of freedom, i.e., the mesonic Q\ operators and the baryonic operators 
(see Chapter 2). Doubling of species can therefore reflect itself only in the number of 
excitations that correspond to these color singlet field, and is therefore determined by the 
global symmetry of the effective Hamiltonian. 

The global symmetry group of the a model depends on the underlying fermion kernel of 
the Hamiltonian. For Nf flavors of naive fermions we get an interaction between nearest- 
neighbor (NN) sites that is invariant under U (N) with 

N = AN f . (6.1) 

This symmetry is too large and is indeed indicative of species doubling. This symmetry 
is also present in the original Hamiltonian, as can be seen by spin-diagonalizing the 
interaction. In order to extract physical consequences from our model, it is imperative 
that it has the same global symmetry as QCD. We achieve this goal by taking a different 
underlying fermion kernel in Eq. (2.35) with the reduced, correct symmetry. We add 
next-nearest-neighbor (NNN) interactions to the kernel and reduce the symmetry to 

U(N f ) L x U(N f ) R , (6.2) 

which is almost the symmetry of the continuum theory. The additional U(1)a is inevitable 
if one starts with a local, chirally symmetric, underlying theory [23]. It can be easily 
removed by hand in the a model by adding U(l) ^-breaking terms [32], and we ignore it. 

The action of the a model with nearest-neighbor next-nearest neighbor interactions is 
[See Chapter 4] 

N r J' _ J' 

s = YJ dT ~ £ Tr Au &u n + f E Tr K^+a) + f E Tr Ka^+^og 

(6.3) 

Here J[ and J' 2 are the couplings of the classical Hamiltonian Eq. (4.29), multiplied by 
2/N c . They scale as l/g 2 N c [recall the discussion in at the end of Section 4.2]. For brevity 
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we drop the primes henceforth. If one derives the fermion Hamiltonian by truncating the 
SLAC Hamiltonian [30], then the couplings J\ and J 2 obey J 2 = Ji/8. If we argue, 
however, that the strong-coupling Hamiltonian is derived by block-spin transformations 
applied to a short-distance Hamiltonian, then we cannot say much about the couplings 
that appear in it. We will assume that couplings in the effective Hamiltonian fall off 
strongly with distance, that is < J 2 <C J\. 

The NN term is invariant under the global U(N) transformation U — > VU, or a — > VaV^. 
The NNN term is only invariant if V^a^V = o M for all /i. This restricts V to the form 

V = exp[t(6 a v + l5 e a A )X a ], (6.4) 

where A a are flavor generators. This is a chiral transformation in U(Nf) x U(Nf). The 
NNN term couples (discrete) rotational symmetry to the internal symmetry, viz. 

a n — > R^a^R, n = Tin. (6.5) 

Here 1Z is a 90° lattice rotation and R represents it according to 

®l Nf - (6.6) 



R = exp 




cr,, 



6.2 Ground state of the next-nearest-neighbor the- 
ory 

The overall N c /2 factor in Eq. (6.3) allows a systematic treatment in orders of 1/N C . 
In leading order, the ground state is found by minimizing the action, which gives field 
configurations that are r independent and that minimize the interaction terms. For the 
NN theory, minimizing the single link interaction 

E=^Tr<r 1 (T2 (6.7) 

allows us to construct the vacuum by placing <j\ and a 2 on the even and odd sites. 

We impose a uniform baryon density, B n = B > 0, on the effective Hamiltonian by setting 
a fixed value of m > 2Nf on every site. To minimize the single-link energy (6.7) we first 
choose a basis where 

" = A =(V-iL)- (6 ' 8) 
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The analysis in chapter 5 shows that at the classical level the ground state is locally 
degenerate, but taking into consideration 0(1/N C ) fluctuations lifts this degeneracy. As 
a result we can choose the the value of a-i with out loss of generality to be 



0"2 



/ WiV \ 

-l N -rn 

V l N _ m J 



(6.9) 



This is a Neel structure, with two sublattices. The vacuum expectation values (vevs) of 
the even sigma fields break U(N) to U(m) x U(N — m) and then the vevs of the odd 
sigma fields break the symmetry further to U (2m — N) x U(N — m) x U(N — m). 

Now we add the NNN interactions to the effective action. At the classical level, they do 
not by themselves remove the classical degeneracy of the NN theory (for nonzero density). 
We have to introduce the 0(1/N C ) fluctuations first in order to stabilize the Neel ground 
state of the NN Hamiltonian. More precisely we assume that 

< J 2 < J ± /N c . (6.10) 

This means that we take as our starting point the (globally degenerate) vacuum deter- 
mined in Chapter 5 for the NN theory with 0(1/N C ) corrections. The NNN interaction 
is a perturbation on this vacuum and its excitations. 

We begin, then, by assuming a Neel ansatz that minimizes the NN term in the action (6.3). 
The NNN term acts within each of the two sublattices. Writing <r e o for the sublattice 
fields, the NNN contribution to the energy per 2x2x2 lattice cell is 

E NNN = y Zl Tr baa^a^] • (6.11) 
" a=e,o pL 

Here a a is a global unitary rotation of the solution to the NN theory given by Eqs. (6.8)- 
(6.9). 

We can find a lower bound for E NNN . Writing in each term E° = a a and = a ajli cr Q;^ , 
we have 

E NNN ='^Y, T ^ a i^T- (6.12) 

Note that S 12 are unitary rotations of the reference matrix A. Each term in Eq. (6.12) 
may be bounded from below by allowing these unitary rotations to vary independently 
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over the entire U(N) group. This is just the minimization problem posed in Eq. (6.7) 
above. The solution is given by Eqs. (6.8)-(6.9), whence the bound 

E NNN > 3J 2 (4m - 3N). (6.13) 



We minimize E NNN by writing an ansatz for a a that saturates the lower bound. At this 
point we choose to work in a basis where 75 is diagonal, 

75= (Y (6-14) 



N/2 



Our ansatz is 




N-m J 










— IjV-m 








ljV-TO 



(6.15) 



(Note that A G is a rotation of A e .) This is a global rotation (via U) of a particular 
configuration that minimizes the NN energy, as seen above. We further suppose that U 
takes the form 

with u e U(N/2). 

The minimization of the energy via the ansatz is a problem of vacuum alignment [83]. 
We begin with the unperturbed NN problem, where the symmetry group G = U(N) is 
broken to H = U (2m — N) x U (N — m) x U (N — m) . Our reference vacuum is given by 
a e>0 = A e o , giving a specific alignment of H as the invariance group of this vacuum. We 
determine the rotation matrix U [within the ansatz (6.16)] that minimizes the energy of 
the perturbation E NNN , which we write as 

E NNN = ^ Tr Ka^A a al, (6.17) 

where = U^a^U = a^. 

We denote the generators of H collectively as T and the remaining generators of G as 
X. The T matrices commute with both A e and A G , while the X matrices do not. The 
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symmetry-breaking term in the energy is E^nn, given by Eq. (6.17) in terms of the 
rotated Hermitian matrices a M . We project each a M onto the T and X subspaces, giving 
the decomposition 

a. 



— T — X 



Using the invariance of A 6i0 under T, 



A a ,al = 0, 



and the orthogonality of the T and X subspaces, 



Tr 



a T a x 



0, 



we have 



E 



J 2 



NNN 



£Tr \(al) 2 + A ajA a aJ 



i,a 



(6.18) 



(6.19) 



(6.20) 



(6.21) 



Following the block form of A ej0 , we divide the broken generators X into 3 sets, denoting 
them as X a with a = 1, 2, 3. Their structures are respectively 



X x 
Xj 




X 2 


xt 



and 




X 3 



(6.22) 



X 3 T 

Writing a* = Xi a a^ a , the following relations can be proved easily: 



[Ai,a^] =0 {A 2 ,c^} = 
{Ai,«?} = [A 2 ,ojf»] =0 

{A 1 ,a^} = ° {A 2 ,< 3 } = °- 

Using these together with the further orthogonality conditions, 



(6.23) 



Tr 



"jit "/i 



0, 



we bring Eq. (6.21) to the form 

E NNN = J 2 £T: [(^ 2 _(-x 3 y 



(6.24) 



(6.25) 
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The rotation U given in Eq. (6.16) saturates the lower bound for the energy. We will 
proceed to prove this for the case to > 3N/4. In view of Eq. (6.16) we can write a M in 
the form 

° a ' L ) , (6.26) 




where a M = u^a^u. It is straightforward to check that for to > 3N/4 we have 



? = 0, (6.27) 

(6.28) 




with a'^ composed of the first (2m — 3N/2) columns of <r M , 



= f {a lx ) pq for g = l,...,2m-3iV72, 
pq 1 else. 



(6.29) 



The energy is 



2m-3N/2 

E NNN = 2J 2 J2 \(^) Pg \ 2 = 2J 2 E E (^) = 3J 2(4m - 3iV), (6.30) 

i,pq i q=l 



which is exactly the lower bound. Note that according to Eq. (6.30), any value u within 
U(N/2) saturates the lower bound. Nevertheless, some values of u give the same vacuum 
as does u — 1 . More precisely, all the elements u EU (2m — 3N/ 2) xU(N — m) xU(N — to) 
act trivially in Eq. (6.15). Vacua that are associated with different nontrivial choices of u 
are in general inequivalent, and give different realizations of the U(Nf) x U(Nf) symmetry 
of the theory. Since these vacua are not related by symmetry transformations, there is 
nothing to prevent lifting of the degeneracy in higher orders in l/N c . In the sequel, we 
set u = In/2- This gives the vacuum with the largest symmetry accessible via the ansatz 
(6.16). 1 

The simplicity of this calculation depends on the assumption to > 3N/ 4. For to < 3N / 4, 
both a* 3 and are non-zero. Moreover the index structure of the projections is more 
complex. 2 Therefore in these cases we resort to numerical minimization of Eq. (6.17) 
over u G U(N/2). In each case we find that the ansatz (6.15)-(6.16) yields a minimum 



1 We have found, however, one exception where another ansatz, with different symmetries, gives larger 
symmetry than that of Eq. (6.15) for the case of B — 2 and Nf — 3. 

2 An exception is the m = N/2 case (B = 0), which was solved in Chapter 5. 
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that saturates the lower bound (6.13). In view of what happens for m > 3N/4 this may 
be only one point in a degenerate manifold of ground states. We emphasize that the 
degeneracy of these vacua is not related to the global U(Nf) x U(Nf) chiral symmetry. 
It is an accidental global degeneracy of the ground state. 

Upon calculating the a fields using Eqs. (6.15)-(6.16), it is straightforward to ascertain 
the symmetry of the vacuum. The U{Nf) x U (Nf) generators that commute with both a e 
and a Q form the unbroken subgroup of the NNN theory. The rest are broken generators 
that correspond to Goldstone bosons. 

The symmetry of these ground states is summarized in Table 6.1. In general, both chi- 
ral symmetry and discrete lattice rotations are broken; in some cases a symmetry under 
rotations around the z axis survives. Note that if we remove the unphysical axial U(l) 
symmetry from the a model, all its realizations will also drop from Table 6.1, namely, 
there will be no axial U(l) unbroken symmetries (third column) and no Goldstone bosons 
corresponding to a broken axial U(l) (fourth column). Anticipating these simple modifi- 
cations we see no reason to proceed in this direction. 

Table 6.1: Breaking of SU(Nf) L x SU(Nf) R x U(1)a for all baryon densities (per site) 
accessible for Nf < 3. Results for B = are from Chapter 5. 



Nf \B\ Unbroken symmetry Broken charges 

= 1 

11 - 1 

2 UilU 

SU(2) V 4 

1 U(l) l3 6 

2 2 SU(2) V 4 

3 U(l) l3 6 

4 SU(2) L x SU(2) R x U(1) A 

SU(3) V 9 

1 U(1) Y x SU(2) V 13 

2 U(1) Y 16 

3 3 SU(3) V 9 

4 17(1) /3 x U(1) Y 15 

5 U(l) h x U(1) Y x U(1) A , 14 

6 SU(3) L x SU(3) R x U(1) A 
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In addition we note the following: 

• Since the baryon background is fixed, we cannot tell whether the U (1) corresponding 
to baryon number is broken. [The U(1)b group acts trivially on Q\ and on a n .] 

• For each value of Nf, the case B = 2Nj is a completely saturated lattice. Each 
site is in a singlet under the chiral group, and there is no spontaneous symmetry 



• The axial U(l) is not a symmetry of the continuum, and must be broken by hand. 
Wherever it appears in Table 6.1 it should be neglected, whether as an unbroken 
symmetry or as a broken charge. 

• The U(1)a' appearing in the table for Nf — 3 is not the original U(1)a group but 
rather is generated by 75 <8> A' with 



This is the only case where an axial symmetry survives spontaneous symmetry 
breaking. If U(1)a is broken by hand, so is U(1)a'- 

For all non-zero densities short of saturation the vacuum breaks rotational invariance. 
This is easily checked for m > 3N/4 by noting that the ansatz for o e o fails to commute 
with the rotation operator (6.6). For the other cases this is easily checked numerically. 
(In some cases a discrete symmetry around the z axis remains unbroken.) Since this is not 
a continuous symmetry it will not give rise to additional Goldstone bosons. The broken 
rotational invariance will of course affect the excitation spectrum. In particular, whereas 
the NN theory possesses excitations with linear and quadratic dispersion relations [60], 
the NNN theory will produce anisotropic dispersion relations, like those seen in [84] and 
Section 6.3 below. 

In comparing our results to those of continuum CSC calculations, one must keep in mind 
that we study systems with large, fixed, and discrete values of B, rather than with large, 
continuous /i. Moreover, we use large-iV c approximations which necessarily ignore the 
discrete properties of the SU(N C ) group that are essential to baryons. Quantum effects at 
finite N c , treated correctly, should yield effects that are not accessible through the 1/N C 
expansion. 



breaking. 



I 1 

A' = 

V 




(6.31) 
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The values of (N c , Nf) that are of interest for CSC are N c = 3 and Nf = 2 or 3. In the 
two- flavor case the favored qq condensate is a flavor singlet and a color triplet, so that 
while color is partially broken, chiral symmetry is unbroken. We do not see this for any 
density at Nf = 2. Plainly our results are due to a qq condensate; whether there is a qq 
condensate as well cannot be ascertained. 

For Nf > 2 the situation is similar. Schafer [34] has considered the color-flavor structure 
of the condensates that arise for N c = 3 and Nf > 3, and he has found that both color 
and flavor are partially broken, with a condensate that locks one or more subgroups of 
the flavor group to the color group. Since we work at large N c , we should stand the 
argument on its head. A plausible qq condensate would lock successive subgroups of the 
color group to the flavor group and hence to each other, leaving unbroken the diagonal 
SU (Nf) L+R+c and some leftover color symmetry. Judging by the global symmetry of the 
vacuum, perhaps we see this for (Nf — 3, B — 3). The other cases could conceivably arise 
from a combination of qq and qq condensates, but whether the latter actually occur is an 
open question. 

Finally we note that according to diagrammatic power-counting arguments [85, 86], CSC 
is suppressed in the 't Hooft limit, where N c — > oo and g 2 N c is kept fixed. Deryagin, 
Grigoriev, and Rubakov showed that in the perturbative regime of small g 2 N c , a Fermi 
sea is unstable towards formation of a chiral density wave with momentum \p\ = 2pp, 
where pp is the fermi momentum [85]. Our ansatz does not permit this possibility. 

6.3 Excitations in the next-nearest-neighbor theory 

The Goldstone bosons of the NN theory were discussed in above. There are 2(N — m) 2 
bosons of type I with uj ~ Ji|fc| at low momenta; these are generalized antiferromagnetic 
spin waves (and are the only excitations at zero density). There are also 2(2m—N)(N—m) 
bosons of type II, that derive their energy from quantum fluctuations in 0(1/N C ). These 
are generalized ferromagnetic magnons with uj ~ (Ji/N c )\k\ 2 . 

Now we calculate the effects of the NNN interactions on the spectrum. In view of 
Eq. (6.10), the NN contributions to the propagators, found in above, remain unchanged. 
In particular we can take over the self-consistent determination of the self-energy of the 
type II bosons. We need consider the NNN contributions to the propagators in tree level 
only. We proceed to calculate these for m > 3N/4. In these cases, the calculations sim- 
plify and we perform them analytically. 3 We believe that the spectra of the other cases 

3 An exception is (N — 12, m = 10), where we have no analytic solution. See below. 
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have similar features. 

In the NN theory the a fields represent fluctuations around the vacuum (6.15) with U — 1. 
We parametrize them as in Eqs. (5.24)-(5.25). If <fi = 0, we have <7 eve n,odd — A evcnj0dd , 
which is the ground state of the NN theory. We adapt Eqs. (5.24)-(5.25) to the NNN 
theory by rotating them, 

a -> UaU\ (6.32) 

with U as given in Eq. (6.16). Now = corresponds to the ground state of the NNN 
theory. 

We substitute Eqs. (5.24)-(5.25) into the action (6.3). The rotation U disappears from 
the kinetic term and from the NN interaction — they are both U (N) invariant. This means 
that the bare spectra found in Chapter 5, when U was absent, remain intact. 

We write the NNN energy as 

-Emm = — J"^ Tl a fJ. Cr a,Na il Cr a ,N+2p,, (6.33) 

where a = (even, odd) and N denotes a site on the corresponding fee sublattice. We 
rescale <f> — > (j)/y/N~ c and expand Eq. (6.33) to second order, 

+0 ( 1) . (6.34) 



We have defined a M = Wa^U, and A^ 1,2 ) correspond to the linear and quadratic deviations 
of the a fields from their ground state values. The latter are given by 

/ -2 X \ ( -2 X X ] \ 

A™ = -2?r , = -27rx f -27T7rt t (6.35) 

V -2x f -2vrt / V +20V J 

on the even sites, and 

-2 X \ ( -2 X X j -2xvr t \ 

A£>=V\ -2vr V\ A(, 2) = 1/ -27rx t -27T7rt 
-2x t -27T+ / V + 2 0V / 

(6.36) 
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on the odd sites. Here 



V 



1 
1 
0-10 



(6.37) 



is the matrix that rotates A even to A odd . It is easy to show that the terms linear in 
vanish. 

In view of the block structure of A eveno( jd and of U, as given in Eqs. (6.15)— (6.16), and of 
a^, it is convenient to decompose x f° r m > 3Af/4 as 



X 



Xi 

X2 



(6.38) 



Here Xi nas N/2 rows, and xi has 2m — 3N/2 rows. Both have N — m columns. Sub- 
stituting into Eq. (6.34) and omitting the ground state energy we find that the 0(1) 
contribution of the NNN energy depends only on xi- The tt and X2 fields do not enter 
the NNN energy at this order and remain of type I and of type II, respectively. 



We now define a new N x N matrix x that contains only xi, 

( xi \ 



X = 





\x\ o J 



(6.39) 



and use it to write 



jvefcc 



E 4Tr 



M A e a M A e {x e N f + a^Ka^K (Vx° N V^ J. (6.40) 



-4Tr 



Next we expand x = X v ^ r ' ', where x n are rea l an d ^ v are the hermitian generators of 
U(N), normalized to 



tr 



(6.41) 



The form (6.39) of x implies that x n 7^ f° r those generators whose elements (r* 7 )^ are 
non-zero for a G [1, iV/2] and j3 G [m + 1 , N] or vice versa. Thus 



E min — ^ ^ 

aN 77,77' 



x N (N a r r x* + £ XN w x & 



2/i 



(6.42) 
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where 



(N r' 



With a Fourier transform, 



= J 2 Tr [Ha^r^aJ , 

= j 2 Tr [vrw j <a l yr r >'v j <a IJi 



x% 



N s (3 k 



Xk 



art ik-N+iwr 



(6.43) 

(6.44) 

(6.45) 
(6.46) 

(6.47) 



we write the energy in momentum space as 

£nn„ = E \N e + N? + M e(&) + M*(k) 



keBZ 
ui>0 



+X°\ N + Nj + M (-k) + Ml(-k) 



X- k - 



(6.48) 



HereM a (fc) = E M M a ^. 

The NN action, including the time derivative and 0(1/N C ) self-energy, was written down 
in Chapter 5 in terms of the Fourier transform \ °f the (2m — N) x (N — m) matrix field 
X [Eq. (6.38)]: 



E Tr " Ei,*) + (-^ - Ei, fe ) xixi* 

u;,fc 



-E 2 , fe (xtX- k + cc.) 



(6.49) 



The self-energies E a are of order J\/N c and depend on iV and m. We set X2 = and 
repeat the steps leading to Eq. (6.48) to write S nn in terms of xl, 

(6.50) 



S nn = E XVKeeXl + X°\KooX- k + Xf KeoX-k, 

k 



with the matrices K ee , K OQ , and K eo given by 



(K ee ) 
(K eo ) 



-Tr 



- -Ei fc 5 w , 
2 ' 



-^Tr 
2 

-E 2 , fc Tr 



AeHr" 



-El fc 5 W , 

2 ' 



(6.51) 

(6.52) 
(6.53) 
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Equation (6.50) is to be added to Eq. (6.48) to give the quadratic action of the type II 
Goldstone bosons xi- 

Diagonalizing the quadratic form is straightforward but tedious. As noted above, only a 
subset of the generators T v of U(N) appear in the expansion of Eq. (6.39). Since Xi nas 
dimensions N/2 x (N — m), there are N(N — m)/2 pairs of generators in the sum, which 
we write (similar to the Pauli matrices a x , a y ) as T^, T^, with rj = 1, . . . , N(N — m)/2. 
Their coefficients are similarly written as Xx,Xy- Thus for each k,r] we have 



v 
Xk 



( Xxk \ 

eq 

Xyk 
or\* 

Xx, — k 
oh* i 
V Xy'-k J 



(6.54) 



and the action is 



S =J2T,XkG k 1 Xk, 



The inverse propagator is 



( -£i, fe + n e + m e 

UJ 



lu>0 k 



-UJ 



V 



o 



-T, 1>k + n e - m e 

-S ljfc + n + m 



(6.55) 



\ 

S 2 ,fc 
— UJ 

-Si )fc + n Q - m J 

(6.56) 



where 



in. 



- -2J 2 ]Ttr [{Tlfa,K a a,K] 5™', 



(m a ) W = 2J 2 £tr [Tla.Tia, 



cos k^. 



(6.57) 
(6.58) 



Si )2 and uj contain a factor of 5 v ^i. If we write the 4x4 matrix Gf, 1 in terms of 2 x 2 
blocks, 

' A B ^ (6.59) 



Gk 1 



C D 



then its determinant is easily calculated via 

= \C\ \B- AC- l D\ 



(6.60) 
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that gives the following determinant. 



oo 2 + (n e + m e - Si)(n + m - S x ) - 
u){n + m — n e + m e ) 



u){n e + m e — n Q + m ) 
uj 2 + (n e — m e — Si)(n G — m — Si) — Y, 2 , 

(6.61) 



For iV = 4Nf < 12 (and m > 3N/A), the matrices n e ,n ,m e ,m all commute, except for 
the case (N = 12, m = 10). Dropping this last from consideration, we are left with the 
values of (Nf,B) listed in Table 6.2. For each case, the simultaneous diagonalization of 
n efi and m Ci0 gives the eigenvalues shown. 



Table 6.2: Simultaneous eigenvalues of n e>0 and m ej0 (in units of U^) for all values of 
Nf and B considered. The resulting spectra fall into four classes. Here v z = cosk z and 




Nf B n e n m e 



m degeneracy class 



3 3 12 v z 



1 1 1 2 v z 

2 1 v± 

2 2 12 v z 



3 5 



2 3 1 



1 

1 v z 
2 

2 v± 



1 v z 
2 

2 v± 



2 1 v± 

1 2 -v z 

2 1 -v± 



2 1 v± 

1 2 

2 1 -Ul 



W_L 4 

w z 4 

u± x3 4 

^ x3 4 

-f_L 4 

-u a 4 

^ 2 

2 

v± 3 

3 

f_L x6 4 

w z x6 4 

-v± x3 4 
x3 4 

x2 1 

w z 2 

2 

fj_ 3 

3 
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The zeros of the determinant (6.61) determine the spectrum u(k), giving (after u — > iu) 



Si 2 - S 2 2 + ^ (n 2 - m 2 + n 2 - m 2 ) - £i(n e + n ) 



±< 



1 I 2 
SiK - n e ) + -(n 2 - m 2 - n 2 + m 2 ) 



+S 2 



^m e + m ) 2 - (n e - n Q f 



1/2 



(6.62) 



Because of the symmetries of Eq. (6.62), the spectra of the various cases shown in Table 
6.2 fall into four classes. We examine each class in turn. 



Class 1: Here there is no contribution at all from the NNN interaction. As shown in 
[60], Si i2 are proportional to J\/N c and for small momenta they are quadratic in |fe|; 

the same holds for the NN energy y^Ei 2 — E 2 2 . These fields remain isotropic type II 
Goldstone bosons as in the NN theory. 



Class 2: Fields that correspond to the minus sign in Eq. (6.62) remain type II, but with 
anisotropic dispersion laws of the form 

u 2 = c 2 k 4 ^l + a5^y (6.63) 

The plus sign in Eq. (6.62) gives a linear dispersion law, again anisotropic, 

uj 2 = 4c l J 2 [k 2 x + k 2 y + {l + b5)k 2 ] (6.64) 

The anisotropy in both cases is proportional to the ratio 5 = J 2 /(12J 1 /N C ) of NNN to 
NN couplings. The coefficients c and c\ are defined as 

Cl - -(§^0. (6,6) 
They are of order J\/N c . The coefficients 

a= ^£i^f! and b= Ik! (6 . 67 ) 

6Ji 2 Cl c 2 6J1C1 V ; 

are of order 10 2 for most cases at hand. These coefficients as well as c and C\ are given in 
Appendix K. 



6.3. Excitations in the next-nearest-neighbor theory 



123 



Classes 3 and 4: Taking k = in Eq. (6.62) we find that the fields that correspond 
to the plus sign get a mass equal to 2J 2 . This is a result of the explicit breaking of the 
U(4Nf) symmetry by the NNN interaction terms; these particles are no longer Goldstone 
bosons. The massless bosons in Class 3, corresponding to the minus sign, are type II 
bosons described by 

u 2 = c\k\ (6.68) 

(This is different from Class 1 where u 2 = c 2 k 4 .) The massless bosons in Class 4 are again 
anisotropic, obeying Eq. (6.64), and are of type I. 

These dispersion relations are correct to 0(5) for momenta of 0(5) or smaller. In all cases 
the dispersion relation to 0(1) for k 2 3> 5 is quadratic and isotropic, unchanged from the 
NN result presented in Chapter 5. Since 5 is a small parameter, this means that in most 
of the Brillouin zone the propagator maintains its NN form. This is the reason why the 
self consistent calculations in Chapter 5 that yield Si j2 do not change when we add the 
NNN interactions. 



Chapter 7 

Summary and discussion 



In this work we have studied the strong coupling limit of lattice QCD with non-zero baryon 
density. We used the Hamiltonian formalism with nearest-neighbor (NN) fermions and 
next-nearest-neighbor (NNN) fermions. Also we do not take the continuum limit, but 
rather we regard this Hamiltonian as an effective Hamiltonian that describes QCD at 
large distances. 

Expanding in powers of the inverse coupling one obtains an effective Hamiltonian for color 
singlet objects. In this work we restrict ourselves to second order in the expansion, where 
the effective Hamiltonian is a generalized antiferromagnet for U (47V/) spins. It has NN and 
NNN interactions, which are invariant under U(4Nf) and U(Nf) x U(Nf) respectively. 
Also at this order the local baryon number B n is conserved and one diagonalizes the 
antiferromagnet for a fixed background of baryon number, which fixes the representations 
of the U (47V/) spins. 

When we write the partition function of the effective Hamiltonian with boson coherent 
states, it becomes tractable in the limit of large N c and Nf, with the ratio r = N c /Nf 
kept fixed. In this limit we extend Arovas and Auerbach's result for zero density and 
dimension d — 2 to three dimensions. We find that for r > r c , the antiferromagnet is 
in a Neel ordered phase, which corresponds to a phase with spontaneous broken chiral 
symmetry. The value of r c in three dimensions is found to be 0.31. When we analyze 
this limit for a periodic configuration of baryon number, where a concentration c = | 
and c = \ of the lattice sites is saturated with baryons, we find that r c is 0.35 and 0.47 
respectively. This means that the Neel phase is shrinking with increasing density. 

For a random distribution of sites saturated with baryons, we average over the sites' posi- 
tions in a quenched approximation. This is valid only if the concentration c of these sites 
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obeys c C 1, and we calculate the free energy to first order in the concentration. Com- 
paring this to the unquenched calculations at c = |, |, we find that the two calculations 
deviate near the phase transition point r = r c . In fact, in the quenched calculation, r c 
decreases with increasing density. 

Despite this deviation, we note that the two approaches are in good agreement far from 
r c . This opens the way to the systematic investigation of the system's properties with the 
quenched approach. For example, relying on the condensed matter analog of magnetic 
impurities, low lying excitations should develop non-zero decay widths and decreased 
velocities for non-zero density. Note that in Section 3.4.2 we were concerned with the 
disordered phase, where the excitations are massive. In certain cases, however, there 
exists a Holstein-Primakoff transformation for the generalized spins in terms of bosons 
(see Appendix H). Using this transformation and taking the limit of large N c and fixed 
Nf, one can perform the quenched calculation in the ordered phase, on lines similar to 
Section 3.4.2. It will be interesting therefore to see how the Goldstone bosons decay and 
slow down in the presence of baryons. 

Next we wrote the partition function using generalized spin coherent states to obtain a 
nonlinear sigma model. In Section 4.3 we analyzed the classical counterpart of the sigma 
model by dropping the Berry phase term. This analysis presumably describes the high 
temperature region of the system, where quantum fluctuations are not very important. 
It also serves as a tool to see how the symmetry realization changes as we add baryons 
to the lattice. Although the classical sigma model can be treated numerically (its action 
is real), we postponed this to future research and used mean field theory instead for a 
toy model that has N = 3. In that case we extracted the classical phase diagram in 
temperature-density plane. We found that U(3) can be broken either to U{2) x U(l) or 
to U(l) x 17(1) x U(l), depending on the density. The order of the phase transition also 
depends on the density, as does the critical temperature, which decreases with increasing 
density. 

We studied the full quantum sigma model in its semiclassical limit of large N c and finite 
Nf. In the vacuum sector, we rederived Smit's result for the lowest-energy configuration 
of alternating B n — ±(m — N/2) sites. The classical degeneracy in the number m is 
removed by quantum fluctuations of 0(1/N C ), and the ground state indeed has B n = 0; 
the U {4:Nf) symmetry of the nearest-neighbor theory is spontaneously broken to U (2Nf) x 
U(2N f ). We extended this result to the NNN theory and found that its U(N f ) x U(N f ) 
chiral symmetry is broken to the vector U(Nf) flavor subgroup. We also found that the 
classical degeneracy of the nearest-neighbor interactions can be removed by the NNN 
terms themselves, with no quantum corrections. These two procedures for lifting the 
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degeneracy led to the same ground state. 

Adding net baryon number to the system, we examined the case of uniform baryon density, 
B n = B = m — 2Nf > 0. When we limit ourselves to the nearest-neighbor theory, we find 
a Neel-like ground state that breaks U(N) to U(2m - N) x U(N - m) x U(N - to). The 
number of broken generators n G thus depends on the baryon number B as 

n G = 2(2N f -B) (2N f + 3B). (7.1) 

When we calculate the low energy excitations of the NN theory, we find two types of 
Goldstone bosons [80, 81]. The type I bosons are ni = 2 (27V/ — B) 2 antiferromagnetic 
spin waves that have a linear dispersion relation. We also have n u = 4B(2Nf — B) type 
II bosons 1 . The numbers rij, rin, and n G obey the Nielsen-Chadha counting rule 

m + 2n H > n G . (7.2) 

We call the type II bosons ferromagnetic magnons since they emerge due to a ferro- 
magnetic alignment of NNNspins. As typical for magnons, their energy is quadratic in 
momentum 2 and of 0(Ji/N c ), where J\ is the coupling of the nearest-neighbor interac- 
tions. 

The mechanism that removed the classical degeneracy at zero and non-zero density, and 
gave rise to type II Goldstone bosons in the latter case, is known in the condensed matter 
literature as order from disorder. It takes place in systems that possess a classically 
degenerate ground state. For example we mention the double exchange model [54] and 
the Kagome antiferromagnet [55]. There, the classical ground state energy is invariant 
under a rotation of local groups of spins, making it exponentially degenerate as in our 
problem. 

Finally when we perturb the nearest-neighbor ground state with the NNN interactions, we 
find that some of the global degeneracy of the nearest-neighbor ground state is removed, 
and that the vacuum expectation values of the sigma fields break chiral symmetry in a way 
that depends on Nf and the baryon number. In all cases with non-zero B, discrete lattice 
rotations are spontaneously broken as well. We summarize the patterns of spontaneous 
symmetry breaking in Table 6.1. 

1 Recall that the local baryon number B ranges from to 2Nf with unit increments. 
2 A possible exception is for N — m = 1, where the energy may depend on a different power of the 
momentum. Correspondingly, the number of excitations will be different. 
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In comparing Table 6.1 to results of continuum CSC calculations, one must keep in mind 
that we study systems with large, fixed, and discrete values of B, rather than with large, 
continuous /i. Moreover, we use large-iV c approximations which necessarily ignore the 
discrete properties of the SU (N c ) group that are essential to baryons. 

The values of (JV C , Nf) that are of interest for CSC are N c = 3 and Nf = 2 or 3. In 
the two-flavor the ground state does not break chiral symmetry. We do not see this for 
any density at Nf = 2. Our results are due to a qq condensate; whether there is a qq 
condensate as well cannot be ascertained. 

For Nf > 2 the situation is similar. The CSC ground state for N c = 3 and Nf > 3 
partially breaks color and flavor, with a condensate that locks one or more subgroups of 
the flavor group to the color group [34]. Since we work at large N c , we should stand the 
argument on its head. A plausible qq condensate would lock successive subgroups of the 
color group to the flavor group and hence to each other, leaving unbroken the diagonal 
SU (Nf) L+R+c and some leftover color symmetry. Judging by the global symmetry of the 
vacuum, perhaps we see this for (Nf — 3, B — 3). The other cases could conceivably arise 
from a combination of qq and qq condensates, but whether the latter actually occur is an 
open question. 

Finally we note that according to diagrammatic power-counting arguments [85, 86], CSC 
is suppressed in the 't Hooft limit, where N c — > oo and g 2 N c is kept fixed. In the 
perturbative regime of small g 2 N c , a Fermi sea prefers pairing of particle and holes that 
gives a chiral density wave with momentum |p| ~ 2/x, where \i is the chemical potential. 
Our ansatz for the ground state does not permit this possibility. 

When we determined how the NNN terms change the dispersion relations for the Gold- 
stone bosons we found that the physics of the NN theory is mostly undisturbed by the 
NNN interaction. At leading order, the properties of the type I bosons (ir) and of some of 
the type II bosons (x 2 ) do not change. The type II bosons grouped in the xi field suffer 
a variety of fates, falling into four classes that appear for different values of Nf and m. 
Class 1 bosons are unaffected by the NNN perturbation. Class 2 bosons split into type 
I and type II, and all gain anisotropic contributions of 0(J 2 / (Ji/N c )) to their energies. 
Here J 2 is the coupling of the NNN interactions in the antiferromagnet. Some of the Class 
3 bosons become massive while other remain unaffected. In Class 4, some become massive 
while the others become anisotropic type I bosons. 

The symmetry of the theory, in all cases, is severely broken by the NNN terms — from 
SU(4N f ) to SU(Nf) x SU(Nf) x U(l) A . Not surprisingly, a simple count shows that 
the total number of massless real fields is far greater than the number of spontaneously 
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broken generators of SU(Nf) x SU(Nf) x U(1)a, as shown in Table 6.1. The particular 
NNN interaction we use is simply unable to generate masses in lowest order for many 
of the particles unprotected by Goldstone's Theorem. This is partly reflected in the 
accidental degeneracy of the ground state, which we mentioned below Eq. (6.16). Just as 
this degeneracy should be lifted in higher orders in 1/N C [beginning with 0(J 2 /N c )], the 
corresponding massless excitations should develop masses. The only particles protected 
from mass generation are the minimal number needed to satisfy Goldstone's theorem (or 
the Nielsen-Chadha variant). 

Another effect that is missing is the mixing of type I and type II Goldstone bosons, 
which is certainly permitted when the NNN interaction is turned on. In the NN case we 
proved that such a mixing is forbidden, since the two types of boson belong to different 
representations of the unbroken subgroup. We present the classification in the NNN 
theory in Appendix L, where it is easy to see that it is less restrictive, and permits mixing 
of the bosons. Whether mixing occurs is a dynamical issue that can be settled only by 
calculating to higher order in l/N c . 

To conclude we note that other recent work [82, 84] on the high density regime of QCD — in 
the continuum — also predicts type II Goldstone bosons and anisotropic dispersion. There, 
the starting point is an effective field theory that describes the low energy dynamics of 
QCD with non-zero chemical potential \i. For /i ^ 0, Lorentz invariance is broken, and the 
field equations become nonrelativistic. This leads to the emergence of type II Goldstone 
bosons. In addition, the ground state in [84] can support a non-zero expectation value of 
vector fields. This breaks rotational symmetry and makes some of the dispersion relations 
anisotropic. 

Possible directions of future research begin with gaining more freedom in fixing the baryon 
configurations. In this work we restricted ourselves to either having sites saturated with 
baryons (in the boson coherent states approach) or with having the same baryon number 
on all sites (in the nonlinear sigma model approach). It is interesting to see what is the 
ground state of the antiferromagnet for other baryon configurations, that have 2Nf > 
B n > only on a sparse sublattice, along the lines of Section 4.3. This can be easily 
investigated for the classical counterpart of the sigma model, whose action is real, by 
Monte Carlo simulations. A general baryon distribution has a classical degeneracy which 
may be removed by order from disorder. This can be seen by calculating 1/N C corrections 
that correspond to quantum order from disorder, as in Chapter 5. Another way to lift the 
degeneracy is with thermal fluctuations. These can be studied in the classical counterpart 
of the sigma model as well. 

Another interesting direction is to fix the chemical potential \i instead of the baryon 
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Chapter 7. Summary and discussion 



number. This requires calculating the grand canonical partition function Zgc(h), which 
is related to the canonical partition function Z({B n }), treated in this work, by 

Z G c= E e^n^ Z ({B n }). (7.3) 

{B n } 

This will make the comparison to results of continuum and other strong coupling calcu- 
lations more straightforward. Again a first step can be to calculate Zqc for the classical 
sigma model. Also it will be interesting to see how order from disorder manifests itself in 
the grand canonical ensemble. 

An ultimate goal is the incorporation of the third-order term in the effective Hamiltonian 
in order to have a theory with dynamical baryons. This Hamiltonian looks like the t—J 
model of condensed matter physics, but is in fact more complicated since the baryon 
operators do not obey canonical commutation relations. An instructive half-measure 
would be to study the antiferromagnet in the presence of a disordered baryon background. 

One can also replace the baryon operators with simpler operators that mimic only some of 
the baryons' features. Reintroducing all the properties of the real operators can be done 
using methods of Hilbert space mapping (see for example [69]), which were introduced 
in the context of atomic and nuclear physics to treat quantum fields that do not obey 
canonical commutation relations. 



Appendix A 



Anticommutation relations of baryon 
operators 

In this appendix we present the calculation of the anticommutation relations of the baryon 
composite operators introduced in Chapter 2. For simplicity we slightly deviate from the 
notations of Chapter 2 and denote color indices by a,b,c,... = 1,2,3 and Dirac-flavor 
indices by i, j, k, . . . = 1, . . . , 4iVy. The baryon operators are given by 

b l3k = e abc r^l (A.l) 
where ^? are fermion operators obeying 

{r i ^f} = 5 ab 5 l3 . (A.2) 

We wish to calculate the following anticommutation relations of the baryon operators. 
Note that we are only interested in operators residing on the same site, and only in the 
anticommutators between b and ¥ . All other relations are trivial. 

We wish to calculate 

{b\ jk , b mnl } = e ahc e a , Vd (4V]Vl>^^Vf + rMi'i>WM a ) ■ (A.3) 

We now concentrate on the second term, which we want to bring to a form similar to the 
first term, with a minus sign. The two will then cancel, and leave the expression we are 
after. We do that by moving the ip^s in the second term of Eq. (A.3) to the left, and 
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Appendix A. Anticommutation relations of baryon operators 



have 



+ba>ch m bil&ac'^^ h n + *P m \^ n ' {hl^cc'^f - kj^c'b^k 



(A.4) 



Moving ip\ a in the first term of Eq. (A.4) leftward, we cancel the first term in Eq. (A. 3) 
and obtain 

+ijfij b nS a > c 5 km 5u5 ac > + ij^fyfi}? a9a 5 k i5 cc i 

-rMM a9a Wbc> + r^f4 a Sn k 6b>c 

-rM^i'^nM} ■ (A.5) 

Using the 5 functions and the two e's in front of the brackets in Eq. (A.5), we have the 
following result for the anti-commutation relations, 

{^Ijfc? ^ranij = 63(ijk,lnm) ~ 2 {&kl&imQ jn + &kl&jmQin + ^jl^kmQin) 

3 {5km$ilQjn $km$inQ jl $km$inQjl ~t~ $kl$jnQim 
~^~3kl$inQjm $jl$knQim $nk$ilQ jrn $nk$imQ jl $jl$imQkn) 

4 (5i m 5j n Q kl ~\~ $in$jlQ km ~l~ $in$jmQ kl ~\~ $jm$ilQ kn 
J r$il$jnQkm ~\~ $nk$jmQil $nj$kmQil) 

~^~$im {QklQjn ~\~ QknQjl) ~l~ ^in {Q klQ jm QkmQjl) 
~^~$il {QkmQjn ~l~ QknQjm) ~l~ $jm {QklQin ~l~ QknQil) 
+8jn {QkmQil + QklQirn) + &jl {QknQim + QinQkm) 
~^~^kn {Q jmQil QjlQim) &km {QjnQil ~t~ QjlQin) 

~\~$kl {QjnQim + QinQjm) ■ (A-6) 
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Here we define 

3(ijk,mnl) $kl$jm$in ~\~ 0~klO~jnO~im ~\~ 0~ jlO~kmO~in ~\~ $jl$kn$im ~\~ $kn$jm$il ~t~ 0~njO~kmO~il ■> 

Qij = (A.7) 

£1=1 

We note in passing that in order to get Eq. (A. 6) we used the following relations 

= -ipj'W + SSij, (A.8) 
[QijM = -5 ik ^, (A.9) 

^l c Qij^i = QkiQij-SuQkj, (A.fO) 
^tQij^i = -QijQik + 35i k Qij + S5 ik 5ij - 5 ik Qij (A. 11) 

As mentioned in Chapter 2, the baryon operators are not canonical anticommuting oper- 
ators, and their anti-commutation relations obey = A n / ^ 8ip. 



Appendix B 

Discretization of Euclidean time 



In this section we stress the importance of the slicing prescription one uses for Euclidean 
time. In particular we work with a nearest-neighbor action defined on a bipartite lattice. 
We show the difference between using the exact prescription that results from the coherent 
state formalism of Section 3.1, and a naive prescription. 

The coherent state formalism results in an action of the form 



$ = X/ ] X/ Z r,n ( Z T,n ~ Z r+l,n) + 6 ^ Z T,n Z T+l,r 
r=l I n 



"I" Z T,m {^ Z T,m Z T+l,m) e ^ Z - 
m 

~ e X Z T,nQnmZ T m + Z T Tn Q 



o* o 

TjITl^T+ljITl 



* e* 
mn z r,n 



f*G- l f. 



fB.l) 



Here the index n(m) runs from 1 to N e (N ) and denotes the sites on the even(odd) 
sublattice. / has two indices; t — 1, . . . , N T — 1, and a site index n that takes values from 
1 to N e or N , 

( <i \ 



/t,t 



Z T,Ne 



\ z °: No i 



(B.2) 
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Appendix B. Discretization of Euclidean time 



In this basis the inverse propagator G 1 is 

= + J T-) • (B.3) 

where the diagonal blocks have dimensions N e N T , and N Q N T . The matrices <9, and D are 

<9r,r' = &r,r'-l — ^t.t'j ^t.t' = &t,t'-1- (B-4) 

times the unity in the lattice indices. Q is a unit matrix in Euclidean time indices multi- 
plied by an N e x N a matrix that describes the interactions between the two sublattices. 
Note that we do not take N e = N Q , to account for systems with vacancies. In order to 
calculate the determinant of the propagator, it is convenient to work in Matsubara space. 
There we have 



d„ x = <5 nX (-l + e-^), 
and the determinant of Eq. (B.3) becomes 

c n l Ne eQ 



D i — e~ 



(B.5) 



det G~ l 



n det V eQ* <1 I 

Ildet (c* n l No ) det (c n l Ne - e 2 g«)- 1 lQ t ) 

n 

Il«) 7Vo " 7Ve det(|c n | 2 l^- e W). 



(B.6) 



Here c n = (1 — e ^ n£ ) + e itJne eA. Next we diagonalize the N e x iV e matrix QQ^ and write 
the determinant as 

We 

\N„-N e TT 



|2 2 2 
" £ 9a 



(B.7) 



detG" 1 = ni c n| 

n 

Here g 2 are the N e eigenvalues of QQ^ . 

We now calculate the determinant for the continuum naive prescription, where the time 
derivative d — > — iu n e5 ntn ', and D — > In this case we have c n = e{ioj n + A), and the 
products in Eq. (B.7) become 



n 



oc exp 



n(i< 



1 2 2 2 

| ~tq a ) oc exp 



Vlog(c 2 + A 2 ) 

. n 



(B.8) 
(B.9) 
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Using the following relation [87] 

oo 

53 log (ul + £ 2 ) = -2 log 2 sinh (3E/2 + constant, 



(B.10) 



where the constant depends only on (3, we find 



det G continuum y 2s mh(p\/2) 



N -N e 



n 



a V 2 sinh (/3w a /2), 



;b.ii) 



Returning to the exact prescription we see that the products in Eq. (B.7) become 



Y[ \c n \ oc exp 



Vlog(4(l- e A)sin 2 (c n 6/2) + e 2 A 2 ) 



n(l c n| 2 " e2 ^) oc exp ^log(4(l-eA)sin 2 Ke/2) + e 2 cu 2 ) 



(B.12) 
(B.13) 



We separate the sums of the form given above to 

53 log (4(1 - eA) sin 2 (cu n e/2) + e 2 X 2 ) = -(3\ + ]T log (4 sin 2 (u; n e/2) + e 2 F 2 ) . (B.14) 

n n 

Here Y 2 = X 2 /(l — eA). We have verified numerically that the sum on the right hand 
side of Eq. (B.14) is given by the sum in Eq. (B.10), by replacing E with X. One can 
understand this point as follows. Divide the sum into two, 



53log(4sin 2 Ke/2) + e 2 r 2 )= £ + £ , 

n \n\<m \n\>m 



(B.15) 



where m N T . In the first sum we replace 4 sin 2 (ci; n e/2) with u 2 , and Y with X. For 
N T — > oo, m becomes very large, and since this sum is convergent, one can replace m with 
oo. 

The second sum the sin 2 (w„e/2) term is no longer of 0(e 2 ), and Y 2 can be dropped. 
Collecting only the X dependent term we are left with a sum of the form of Eq. (B.10). 
This gives 



det G- 1 = 



1 - e-^ 



N -N e 



n 



o/3A/2 



2 sinh (puj 2) J ' 



(B.16) 



with u a = \J\ 2 - q 2 a . 
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Appendix B. Discretization of Euclidean time 



As a result we see that the effective action received after integrating the bosons S c g ~ 
log det G~ l obeys 

SX ct = SX tinuum - (N e + N )\/2. (B.17) 

For practical reasons it is more convenient to perform calculations with the naive pre- 
scription of d — > — icu, and in order to restore the correct result we use the rule of thumb 
given in Eq. (B.17). 



Appendix C 



Minimizing the 
next-nearest-neighbor term 



The sfi signs are defined only when the a v are written in the basis M v = T A <g> A a = 
p a <&aft®\ a . We choose a chiral basis for the gamma matrices, so that 75 = p 3 , «j = p 3 cr l , 
and (3 = p 1 . The energy (4.76) is a sum of squares, 



(C.l) 



with the constraint J2 T ){ a ' r> ) 2 = N/2. The coefficients A v = Z) M take on the values 
{—3, —1, 1, 3}. The minimum of e' occurs when all <j v are zero except those corresponding 
to A v = -3, namely, those for M v = (3 ® A a = p 1 ® X a and M v = ^75 ® A a = p 2 ® A a ; the 
energy is independent of these 0^, and will be equal to 



e' = -3]TK) 2 = -3JV/2. 



(C.2) 



Thus the set of solutions can be written in the form 



cr = 



/ [/ \ 

u 

W 

V c/t y 



p 1 + ip 2 



u + 



p 1 — ip 2 



c/t. 



(C.3) 
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Appendix C. Minimizing the next-nearest-neighbor term 



Recalling that <7q = 1, we have UW — 1, so U e U(Nf). A chiral rotation cr — > VVoV, 
with 

/ C/ \ 

! =-(l+p 3 )®f/ + -(l+p 3 )®l, (C.4) 

V o o 01/ 



turns do into p 1 = 7o, which is invariant only under vector transformations generated by 
1 <g> A a . 



Appendix D 

The function $ for mean field theory 
of the U(3) toy model at non-zero 
density 



In this section we give <£> for all configurations that appear in Table 4.1 with non-zero 
baryon number. The lower indices of the magnetizations correspond to the numbering in 
Table 4.1. 

As for zero density, we reduce the number of independent variables in <3> by substituting 
for some of them with their corresponding MF equations. Also we assume that all magne- 
tizations commute, and are therefore diagonal in some basis. This restriction is justified 
for all configurations except no. 2 where we have to assume this as an ansatz. 



Configuration no. 2 



In this case the form of the MF equations does not allow us to take the magnetizations 
as diagonal, and we make the ansatz that in this case 

/ 1 \ 



rrii 



2 = 1,2,3,4, 



(D.l) 



1 
V -2 J 

which means that we allow fif ^ only. This gives the following form for $ as a function 
of ill 

$/N s = \/3/4K(ni - 3^3 - iii + /x 4 ) - K/A(3fj, 1 fj, 2 + 3/i 3 /i 4 + 6/xi/i 3 ) - Lp(6K\/3iii), 
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<p(Xi) 



A*2 
A*3 



, fe- Xl -l l\ , /e^-l 1 



X? 



+3 log 
+3 log 



' e -(§X3 + §S 2 ) _ ]_ 

(§x 3 + |x 2 ) 2 ' §x 3 + \x 2 



+ 



e (|a;i + 3X4) 



+ 



7X4 



-/(-2K v / 3(2/i 1 + /i 4 )), 



where we have dropped the "8" index from all /i's, and 

, , _2_ fS e x {l-2/x + 2/x 2 )-2/x 2 
f[X) ~VZ e*(l - 1/x) + 1/x ' 



(D.2) 



(D.3) 
(D.4) 



(D.5) 



We find that there is a first order phase transition, similar to the zero density case, at 
K r = 0.94. 



Configurations no. 3, 5, 6 

Here taking only diagonal magnetizations into account is justified and we obtain the 
following $'s. (From here on we denote (/i 3 ,/x 8 ) by \x, for all magnetizations) 

1. Configuration no. 3 

There are three independent magnetizations j2f, i — 1,2,3 and we have 

$/iV s = V3/2K(3 f ,l-2 f xl- P Ll)-3K/A(2fl 1 -il 2 + il 2 -pL 3 )-<p(6Kil i ), (D.6) 



1, 



lo § ( — — S ~ „ 3 + 31og 2^ „ 3 2„ „ 3 



+3 log 



§(x +j3 + 2x+,i)(x| + 2x\) + 2a:_,i)(x| + 2x1) 

(D.7) 

^ = f 3 (AKfil,AKfi 8 2 ), fi^MAKulAKfil), (D.8) 
^ = f 3 (6K^ 2 ,6K^ 2 ), »l = f 8 (6K»l6K»l), (D.9) 

where / 3;8 are given in Eqs. (4.64)-(4.65). Also x± ti = V3xf ± x\ , and we use this 
definition in further cases as well. 
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2. Configuration no. 5 

There are three independent magnetizations fa; i = 1, 2, 3 and we have 

$/iV s = y/3/2K(3f4 + 2/4 + /if) - 3K/4 (2fa ■ fa + fa • fa) - <p(6Kfa). (D.10) 
Here ip, and fa >3 are given by Eqs. (D.8)-(D.9). 

3. Configuration no. 6 

Here we have two magnetizations fa; % = 1,2 and we have 

$/iV s = VS/2K(fi 8 2 + pLf)-3K/Afa-fa-<p(6Kfa), (D.ll) 



/ e ^+,2 _ x e x-,2 _i\ I e \x +A _ x e \x- A _ x 

l0 § 3 3~ + 3 l0 S — 3 2 3" 

y X +t 2%2 X -,2 X 2 / \ 3 X +,1 X 1 s x -,l x l 



(D.12) 

//? = f^K^Kul), n\ = f^Kul^K^). (D.13) 

In these three configurations we find that $ has three equivalent nontrivial minima on 
the lines Eqs. (4.32)-(4.34) in the /i 3 -/i 8 plane. These become lower in energy than the 
symmetric solution for K > K c . The phase transition is first order and K c is 1.03, 1.14, 
and 1.72 for configurations 3, 5, and 6. 



Configuration no. 4, 7, 8 

Here we also find that one can take the two magnetizations to be diagonal. In all these 
configurations we have two magnetizations fa; i — 1,2 

1. Configuration no. 4 
$ is given by 

®(fa)/N s = -3Kfa-fa-2KV3(^ 1 + ^ 2 )- l p(6Kfa), (D.14) 



' f e x+,i-i e -*-,i_i\ f e x+,2_i e x +. 2 - T 

l0 § 3 3~ )~ l0 S 3 3~ 

\ x +,l x l x -,l x l / \ X +,2 X 2 X +,2 X 2 , 



H\ = f 3 (6Kfii6Kf4), f4 = M6Kfii6Kf4). (D.15) 

In this case we find that the phase transition is of second order. Also, there are six 
nontrivial minima. All of them have the same <£>, and become preferable over the 
symmetric solution for K > 1. Four of these minima are invariant under U(2) xU(l), 
and two are invariant only under U(l) x U(l) x U{\). 
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2. Configurations no. 7,8 

Since these configurations are the two and one dimensional counterparts of config- 
uration 4 [see Table 4.1], the MF equations and <£> are given by replacing 6K in 
the MF equations and $ of configuration no. 4 by AK and 2K, respectively. This 
replacement simply increases the value of K c by factors 3/2 and 3. 



Appendix E 



The matrices C and D in the 0(1/ N c ) 
calculations 



To analyze the matrices C and D n , given by Eqs. (5.7) and (1.18), we begin with the U(N) 
generators M v that lie outside the subalgebra H = U (m) x U(N — m) that commutes 
with 

A - ( V -iL ) • (E l) 

We choose for them a basis M pql and M P92 with p = 1, . . . , m and g = m + 1, . . . , N, 
given by 

(M^) a/3 = l(6 pa 6 qf) + 5 pf3 5 qa ) (E.2) 

(M^ 2 ) a/3 = ^(^V-VW- (E.3) 

Since the coset space U(N)/H is a symmetric space, the commutator [M pga , M p ' q ' a '] lies in 
if; in order for Tr (A[M pqa , M p q a }) to be nonzero, the commutator must have a nonzero 
component in the Cartan subalgebra of H. This is only possible if a ^ a' and (p, q) = 
(p', q'). Thus in this basis C takes the form 

C = i { , ^(N-m) \ (E 4) 

\ — Im(AT-m) U y 

Diagonalizing C gives eigenvalues ±1. The generators corresponding to the basis that 
diagonalizes C are 

(M pq+ U = (M pql +iM pq2 U = 5 pP 5 qa (E.5) 



145 



146 



Appendix E. The matrices C and D in the 0(1/N C ) calculations 



(M^-)ap = (M pql -tM pq2 ) a p = 5 pa 5 q p. 



(E.6) 



As for D: The only anticommutators among the M pq± that do not vanish [note the bounds 
on (p,q)} are between M pq+ and M vq ~ , viz., 



{M pq+ ,M pq ~U = 5 qq 
Noting that M pq ~ = (M pq+ )\ we find that D n takes the block-diagonal form 



jjpq±,p'q'± 



' E Tr 

m(n) 



(m) 







{{M pq± )\M p ' q,± } a ™ 

\ U — ±N-m 



m(n) 



We summarize this by writing 



D= En ° 
E„ 



1 



N-mi 



where E n is an m x m matrix given by 



E n = E (in."*) 

T7a(Tl) 



(E.7) 

(E.8) 
(E.9) 



(E.10) 



(Ell) 



It is easy to prove that the eigenvalues of E range from to 4d. In particular, D is 
positive. 



Appendix F 

Feynman rules to 0(1/ N c ) for the 
effective action of the nonlinear 
sigma model 



The following are the Feynman rules extracted from Eqs. (5.28)-(5.31). A = (e,o) stands 
for even and odd respectively. Latin indices take values in the range [1, N — m] while 
Greek indices take values in the range [1, 2m — N}. The propagators are given in Figs. F.l 
and F.2. The vertices extracted from the cubic and quartic interactions [Eqs. (5.29)- 
(5.31)] are given in Figs. F.3-F.4. 




Figure F.l: Propagator of the it fields. 




Figure F.2: Propagator of the x fields. 
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(6) 




= 2Jd N s N c p S ^SfiS jk Si g [S Ae S C eSBoSDo + (e <-> o)] 7 p _ q , 




(d) 



7T. 



7T 



c 

/.9 




X = 2Jd N JV p SgiSjedpsShf [SAeSceSBoSDo + (e <-> o)] 7 p _q, 



Xe/3 



( e ) 4f-!r>*~ xg = x!f =#• 





7T 



= -Jdj^^SasSjkSleSif [5 Ae S B eSco5De + (e <-> o)] 7 p . 



149 




*03 Xjf} 



= -^ J d^^p8 a p5jf5 g i \5a 5bJco - (e <-> o)] 7 fe 
Figure F.4: Vertices corresponding to the cubic interaction. 



Appendix G 



Integral equations for the self energy 
of the zero modes 



In this appendix we derive the integral equations represented by Fig. 5.3 and solve them 
approximately in order to obtain the dispersion relations of the \ bosons to 0(1/N C ). In 
analogy with the propagators of the antiferromagnetic spin waves we assume the following 
structure for and £ 2 ,fc (we denote the frequency by Q), 

Ei >Jfc = A 1 (tt 2 ;k)+ittB 1 (tt 2 ;k), (G.l) 
S 2 , fc = A 2 (n 2 -k)+i{lB 2 (n 2 ;k), (G.2) 

where Ai j2 ,Bi :2 are real and A\ > A\ for stability. In order to perform the frequency 
integrals we make the following rational ansatz. 

A - -^o^m <0? (G - 3) 

2J 1 d P 2 (Q 2 ;fc) 

M ~ I^Q^ky (G - 4) 

Where Pi )2 (Q 2 ; k) and Qi^(Q 2 ', k) are polynomials of the same order in uj 2 . Thus we have 

(rX] 1 (iQ(l + B 1 )-A 1 A 2 + iQB 2 \ 

\^ k )AB Q2[ ( 1 + Bi)2 _ B 2 ]+A 2_ A 2 y A 2 -iQB 2 -iQ(l + B 1 )-A 1 ) AB 

(G.5) 

Since B li2 ~ 0(1/N C ) we drop them in the denominator. The self-consistent equations for 
A and B then decouple, while only Ai :2 are needed to find the poles of the x propagators 
to leading order. 
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Appendix G. Integral equations for the self energy of the zero modes 



G.l Evaluation of diagrams 



We begin with the diagrams that contribute to S x . In this case we have contributions 
from all three diagrams in Fig. 5.3. Using the Feynman rules we evaluate each diagram 
and present the results for the self-energy matrix's upper left block with x e external legs 
in Fig. G.l. (The result for external legs of x° is the complex conjugate of that) 1 



...Q.. 
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Figure G.l: Schwinger-Dyson equation for Ei. 



Here 

E q = y/T=i*. (G.6) 

The factor of (N — m) is the multiplicity of the diagram, and represents the number 
of particles that flow in the loop and have the same contribution. The diagrams that 
contribute to E2 are only the second and third in Fig. 5.3. The result is given in Fig. G.2. 
Note that in this case there is no multiplicity to the diagrams. This is a result of the 
structure of the vertices. For example the first term in Eq. (5.31) that contributes to the 

iHere we define J d 4 q ee J bz JT„ 



G.l. Evaluation of diagrams 
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Figure G.2: Schwinger- Dyson equation for E 2 . 



second and first diagram in Fig. G.l, and Fig. G.2 has the following form. 



TrrfxV 



E (x e U(x* e Mx hAx*°: 



aj ■ 



(G.7) 



a/3g[l,2m-JV] 
m] 



Its contribution to Ei (of the even sites) comes from the contracting x° and x*° ■ Since this 
contraction is proportional to Sapdjj, a sum over j gives a factor of N—m. The contribution 
to E 2 comes from contracting x*° to x* e - This gives 5^5^, and has multiplicity 1. The 
same reasoning applies to all the diagrams in the figures we present. 

Using the parametrization (G.3)-(G.4) we get the following integral equations 2 . 



Pi(V 2 ;k) r (dq\ d /■« duo 

Qi(n 2 ;k) ~ ( ~ m 'jBzU^J i-oo2T 



+ 



^Pi(a; 2 ;g)/Qi(a; 2 ;flf) 



(uj-n) 2 + E 2 q ' uj 2 + E 2 (cu 2 ;q) 
(" - n) (t? - 72) + i^ 2 ; <?)/Qi(u; 2 ; g) (2 7fc7q7(J _ fc - 7 2 - 7 2 ) 



+ #V ;9 )] [( W - n)2 + 



(G 



2 We have defined q = 2J\du and ko — 2J\dQ. 
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where the three terms come from the three diagrams in Fig. G.l and 



P 2 (^ 2 ;fc) 



du j rc P2{uJ 2 ;q)/Q2{u> 2 ]q) 
2tt uj 2 + E 2 {uJ 2 ;q) 

27fc7q ~ lq-k {ll + 7fc) _ 

- ^) 2 + K_ k 

which comes from the second and third diagrams in Fig. G.2. Here 



x 



E(u 2 ;q) = — 



n 


Pi 


1 p 


Ql 


Ql 


~ N C Q 



(G.9) 



(G.10) 



G.2 Evaluation of frequency integrals 

It remains to evaluate the 0(1/N C ) contribution of the following integrals 



1 2 — 



r % . , (G.i2) 

N c J-oo27T (uj2 + e 2 )[(u-Q) 2 + E 2 } 

h = r Pi — A {U ~ Q) — • (G.i3) 

J-oc 27T ^2 + E 2^ [( w _ ty2 + £2] 

All these integrals are convergent at large ui and can be evaluated with residue calculus. 
Beginning with II, 

ji _ 1 r &l p iQiQU (cu) 

1 N c J-oo27fu 2 Q 2 + ^P 2 ' 1 ' ; 
The poles of the integrands are given by 

c 2 Q V) + -^P V) = 0. (G.15) 

c 

To leading order in 1/N C the roots are determined by either 

u 2 = -^P 2 (0)/Q 2 (0) ee ul (G.16) 



G.2. Evaluation of frequency integrals 
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or by 



QV) = o. 



(G.17) 



Since the polynomials P { and Q { do not depend on N c and Qi(0) can be chosen to be 1, 
the solutions of Eq. (G.17) are all 0(1). Moreover since Q 2 > 0, all these roots appear 
in complex conjugate pairs (00,00*). For defmiteness we choose Imw > and close the 
contours from above. 



The result for P l is 3 
1 2ni 
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+ 0(1/N C ) 



1 P(o) 



2g,(o) 



Mr - fMV 

Wi(o)J VQ 2 (o)y 



+ 0(1/7V C ). 



(G.18) 



The leading contribution to Eq. (G.ll) comes from the poles given by Eq. (G.15). For the 
fourth row we used the fact that all other poles are of 0(1), and in the fact that Q 2 (0) 
can be written in terms of a product of its zeros. For the fifth row we used the definitions 
of Q(0) and P(0). 

The calculation of I2 and I3 is similar, but one has two more poles to consider, at 00' = 
fl ± iE. In the case of I2, the leading contribution is again due to the pole at uj only. 



1 2 
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W r 27 



Pi(Vn)Qi(Un)Q&il 



00 r, 



1 



2i\m ijj n fl 
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mJ \ r. 



CO 



*2' 



X 



(O0 r , 



tt) 2 + E 2 



3 Note that the denominator Q 2 oj 2 + 1/N 2 P 2 > is a polynomial in J 1 and can be written as Jlm^ 2 — 
w^)(w 2 — u>^), where ±iv m are the (complex) zeros of Eq. (G.15). We also have Q 2 — Yl' m ((^ 2 — w^)(w 2 — 
), where the product includes only the zeros of Eq. (G.17). 
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Pi{<J)Q&du') 



1 Pi(0) 
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2Qi(0) 
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(G.19) 



The integral J3 is written as 
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(G.20) 



29? + E 2 

Here we used the the fact that u n ^ Q are the zeros of Q 2 , which means that the residue of 
these poles is zero. 

Using our freedom to normalize the polynomials <5i,2(0) > 0, we choose Qi,2(0) = 1 and 
write the integrals to leading order as 

1 Pi(0) 



2 V /P 2 (0)-P 2 2 (0)' 
1 P(0) 



1 P 



(G.21) 
(G.22) 
(G.23) 



2Q 2 + P 2 ' 

Thus we find that to 0(1/N C ), all integrals depend only on the values of the polynomials 
at cu = 0. Finally Eqs. (G.8)-(G.9) simplify to 



Qi(9 2 ,k) 



- (N - m) LM. 



cosh 9„ ( 1 + 



27fc7g7g-fe ~ 7q ~ 7fe N 
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- E n 



P 2 (ft 2 ,fc) 

Q 2 (ft 2 ,fc) Jbz \ i- 
Where we have defined 



sinh Or, 



tanh 9 r 



EU + ft 2 



7q-fc 



A(0,g) 



(G.24) 
(G.25) 

(G.26) 



9 Pi (0,9)' 

Here we also point out that Si j2 have zeros at /c = 0. For example the function multiplying 
the the cosh in Eq. (G.24) is exactly zero at ft = 0, and q = 0. This is also the case for 
the rest of Eq. (G.24), and for the function multiplying the sinh in Eq. (G.25). This is 
a direct result of the Ward identities regarding the U (N) generators that correspond to 
fluctuations in the x directions. 

Taking ft = 0, and dividing Eq. (G.25) by Eq. (G.24), we have 

d 



(N - ra)tanh0 fe = 



( !*1 

bz \ An 



I 2 (q, k) sinh 9 q 



Here we have further defined 

V(k) 
h(q,k) 

h(q,k) 



( 

bz \ An 



r Idq 
Jbz \ An 



h(q, k) cosh^q — rj(k) 



(G.27) 



lq-k ~ Ik 



+ E a 



= 1 + 



^Iqlklq-k ~ l\ ~ lV 



E l- k 



^Iqlk - lq-k (l\ + ll 



F 2 

^q-k 



lq~k 



(G.28) 
(G.29) 

(G.30) 



A solution of this equation must obey 



/bz(S) 7 i(^fc)cosh0 q -77(fc)>O, 
so that the stability condition (5.49) will be satisfied. 



(G.31) 
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After solving the integral equation (5.52) one can extract the spectrum of the x bosons by 
calculating the poles of the propagator (5.48). These poles are given by equation (G.15). 
As seen there are two kinds of solutions. The 0(1/N C ) poles of the form (G.16) give rise 
to a low energy band [see Eq. (G.10)] 



2 T d 

±iu = k ) ~ p 2 (0, k) (G.32) 



' c 



Using Eq. (G.26) and Eq. (5.52) we have 

2J 1 d(N-m) I —r- 2J 1 d(N-m) s 

±iu = - J - 1 ^ J -a h VI - tanh 2 9 k = - J ~jj >-e k . (G.33) 

The 0(1) poles that solve Eq. (G.17) represent uninteresting massive excitations. 



G.3 Solution of integral equations 

In order to solve the integral equations we assumed that the function 6 q is a function of 
the absolute value q = \q\. Taking f2 = 0, and dividing Eq. (G.25) by Eq. (G.24), we have 
(for d = 3) 



dqq d cos(0) dip h(q, k) sinh 8 q 

(N - m) tanh 9 k = JBZ . (G.34) 

/ dqq 2 dcos(0) dipl\(q, k) cosh# 9 — rj(k) 
Jbz 

and q = q (sin(0) cos(</?), sin(0) sin(<^), cos(0)), where < ip < 2n, < < 7r. The 
integration is over the first Brillouin zone of the fee lattice. The BZ is presented in 
Fig. G.3. Since it is not spherically symmetric, the maximum value that \q\ and \k\ can 
take depend on the direction of the vectors. In direction (1, 0, 0) it is 2tt, whereas in 
direction (1,1,1) it is V3n. The maximum allowed value is y/En, and is attained for 
example in direction (0,1,2). In order to perform the momentum integrations we have 
generated a spherical mesh of the 6 coordinates k and q that has a radius of VoV while 
restricting it to the BZ by multiplying all integrands by a step function that is 1 if a point 
in this six dimensional space obeys the following restrictions, and zero otherwise. We first 
restrict |/c x ,j/,z| < 27T and |? x ,j/,«| < 27T, and further 

\q x + % + Qz\ < 3tt, I - q x + q y + q z \ < 3tt, 

1 5s - Qy + Qz\ < 3tt, \q x + q y - q z \ < 3tt, (G.35) 



G.3. Solution of integral equations 
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Figure G.3: The first Brillouin zone of the fee lattice. The side of the outer cube is 
An. The BZ surface is composed out of six small and eight larger faces. The dis- 
tance between the origin and the of the faces are \/3tt and 2ir. (Reproduced from 
http://140.114.18.41/micro/chap2/fig/fig0206.jpg) 



and the same for k. This forms the shape shown in Fig. G.3. Using the spherical as- 
sumption we reduce the integral equation to be one dimensional (with regard to the 
integration). We define 

ai(q,k)=q 2 / dcosU) diphiq, fe), a 2 (q,k)=q 2 / dcosU) d(pl 2 (q, fe), 

JBZ(q,k) JBZ(q,k) 

(G.36) 

where the boundaries of the integration depend on the value of q and k, hence the notation 
BZ(<7, k). Finally we have 

/ dq a 2 (q, k)sinh 6 q 
(iV-m)tanh# fc = -^ff . (G.37) 

J dq a\(q, k) cosh 9 q — i](k) 
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Eq. (G.37) is only an approximation to Eq. (5.52). On the right hand side, there is an 
explicit dependence on the direction of k (via and 77), while on the left hand side the 
dependence should vanish by assumption. We first examine the solution to Eq. (G.37) 
for k in the direction (0, 1,2), so that \k\ can take all values up to y/bn like \q\ on the 
left hand side. This solution is presented in Fig. 5.4. The next step is to see how this 
solution behaves in other directions. We do not expect the equation to be satisfied for 
all directions. More exactly we expect that the behavior away from the BZ boundary, 
near k = 0, will indeed is spherical, so that the equation is approximately correct for this 
region. We give the results of this check in Fig. G.4. Here the solid line is the original 



1 




Figure G.4: Inspection of the integral equation with the spherical approximation in dif- 
ferent directions. Here we show the case of c = 5. 

solution of the integral equation, in direction (0, 1, 2), with (9, 0) = (tan _1 (0.5), 0). There 
\k\ has a maximum value of \^Eir. The stars, boxes, and circles are for directions of 
(6,<j>) = (tan- 1 (l/2),7r/10),(tan- 1 (2/5),0),(7r/10,0). The dashed line is for (0,0) = 0. 



G.3. Solution of integral equations 
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In all these cases, the solution deviates from the original one for high enough momenta. 
As we get further away from the original direction, the region in \k\ space in which the 
solution is correct becomes smaller. Common to all directions is that near k = they 
agree, which was expected, and points to a spherical behavior there. 



Appendix H 



Hamiltonian approach to order from 
disorder for N — m = 1 



In this appendix we derive the self-consistent equations in a Hamiltonian formulation 
using a generalized Holstein-Primakoff transformation for N — m — 1. The Hamiltonian 
is 

# = -A £ Q a p{n)Q Pa {n + jx). (H.l) 

np,af3 

where the U(N) generators Q a /3(n) obey 

[Qapin), Qys(m)} = (Qas(n)5f S ~ f - Q yf s(n)5 a s) 5 nm . (H.2) 

For N — m = 1 there is a simple representation of these operators in terms of N — 1 boson 
operators boson operators [49] , 



Q(n) a p = — 



( ljv_i-2#t 



1 2,/i - 0t^t _i + 2 0U y 

\ V ' o 

Here is a (N — l)-component boson field defined on site n that obeys 



= Q\n)p a . (H.3) 



{n)i,<j>\m)j 



-^SijS^, Vi,je[l,N-l]. (H.4) 



The relation (H.4) has an important consequence. It means that at the large N c limit the 
field operators 4>% are classical c-numbers, and the problem is not quantum. In the path 
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integral formalism, Eq. (H.4) is equivalent to the fact that the Berry phase of the action 
is proportional to Sb ~ N C J dr (f>\d T (f>i. Indeed at large N c , stationary phase arguments 
hold and the kinetic term is zero, making the system completely classical. 

The Hilbert space is created by repeated applications of $ on the no-quantum state |0), 
<f>\0) = 0. It is restricted to be a U (N) representation with iV c columns and N—l rows, by 
restricting the eigenvalues of the number operator <p^(f> to be smaller than or equal to one. 
This representation is the generalized Holstein-Primakoff transformation. For N — m > 1 
the generalization is more complicated (see the Appendix in [21]). 

We first verify that the relations Eq. (H.2) are satisfied by Eq. (H.3). We have four kinds 
of operator given by 

Qij = ^(^-204), (H.5) 
TV 

Qnn = -f{-l + 2h), (H.6) 

QiN = NcfaVl - n, (H.7) 
Q Ni = N c y/T=h<j>\ = Q\ N . (H.8) 

Here we define the number operator n — J2i <l>\<j>i- We need to verify the following com- 
mutation relations. (All operators are on the same site.) 



[Qij-, Qki] 


= Qirfjk — 


QkjSil, 


(H.9) 


[Qij, Qnn] 


= QiN<)jN - 


- QNj$iN = 0, 


(H.10) 


[Qij, Qun] 


= QiN^jk - 


" Qkj^iN = QiN$jk, 


(H.11) 


[Qnn, Qin] 


= Qnn^n 


— QiN^NN = —QiN, 


(H.12) 


[QiN, QjN\ 


= QiN^jN 


— QjN^iN = 0, 


(H.13) 


[QiN, QNj] 


= Qij$NN 


- QNN^ij = N c (l - n) - N c <j>i4>). 


(H.14) 



All other commutators are automatically satisfied once the above list is satisfied. We 
first note that n (and therefore Qnn) commutes with all operators that do not change the 
number of bosons. As a result the second commutator is zero, as expected from Eq. (H.10) 
(i cannot be equal to N). We now evaluate the rest of the commutators using Eq. (H.4). 

[Qij, Qki] = (-2^) 2 [04, <j> h <j>\] = N^t, <j> k ]<t>\ + N 2 C M&, <t>\\<t>) 

= QilSjk — QkjSil, (H.15) 

[Qij, QkN\ = -iV c 2 [0i0], (j) k ]Vl - n = N c (pi5jkVl - n = QiN^jk - Qjk^iN, (H.16) 
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[QnN, Qin] 
[QiN, QjN\ 
[QiNi QNj\ 



N*[n, <pi]\/l-n = -N c (f)iVl - n = Qnn^n - QiN^nn, (H.17) 
NZ(l-n)[<f> i ,<j> j ]=0, (H.18) 
N^Vl-n, Vl - h</)}} = N 2 C (0,(1 - n)<i>) - Vl - n^fay/l - n) 



= N 2 C (0,(1 - h)4>) - 0, t l v / l^) = iV c 2 [0,(l - n), 0]] 



NJl-n)-N c 



;h.i9) 



Next we write the operator as 







X 

7T 



(H.20) 



The (N— 2)-component field x w iU represent the zero modes and the fields n will constitute 
the antiferromagnetic spin waves. We write the generators on the even sites by inserting 
Eq. (H.20) into Eq. (H.3). As in Section 5.2.1 we replace Q by 

Q'(0) = VQMVi (H.21) 

on the odd sites. V is given in Eq. (5.21). It is easy to see that the algebra (H.2) on the 
odd sites is obeyed by Q' as well. 

Now, we scale — > 4>/yfN~ c , substitute Q and Q' into Eq. (H.l), and expand up to 
0(1/N C ). We are left with the following normal ordered form of the Hamiltonian, 

H cS = N 2 E + N C H 2 + ^N C H 3 + H A , (H.22) 

where 



H 2 
H 3 



nfj. 



^ 1 X] (xli+^iXni^n+p, ll nXn+ (j.,iXni + h.C.^j , 



(H.23) 

(H.24) 



n.fj.,1 



nfj. I i] 



2 {XniXni^n+fj^n + X TJ ,+ / x j jX'n.+/x,j 7r n 7r n+/l + /l-C.^ 

i 

- (xliXn^l+^n+ix + (n <-> n + £)) 



(H.25) 



+- 



(7rJ l 7r Tl 7r Tl 7r Tl+/i - 7rJ l 7r^ +/i 7r Tl 7r TH . A + /i.e.) + (n <->• n + ft) 
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H.l The ground state at 0(1) 



At lowest order in 1/N C , the Hamiltonian is given by H 2 , which does not depend on Xn- 
Next we introduce the spatial Fourier transform 

where N is a site in the fee lattice, and k belongs to its first Brillouin zone. In momentum 
space one diagonalizes H 2 with a Bogoliubov transformation to find 

N C H 2 = 2.hdN c ]T ^l-ll (at a, + b\b q ) . (H.27) 
q 

7q is given by Eq. (5.31) and 



n q | _ ( cosh(/9 9 sinh(/9q 
7rl° / V sinh^q cosh<£>„ J \ b[_ 



(H.28) 



-q J \ ^ox^q y ^ y 

where the even function (p q = <p_ q is given by 

tanh2^ q = 7 9 , (H.29) 

and the fields a q and b q obey the usual commutation relations of ladder operators. The 
0(1) ground state |0) is thus given by 

a q |0) = &,|0) = 0. (H.30) 

To 0(1) the excitations of |0) are antiferromagnetic spin waves with a linear dispersion 
relation. The ground state |0) has a local degeneracy, corresponding to creating an ar- 
bitrary number of x bosons on any site (creation of a x boson costs zero energy at this 
order). It is therefore more appropriate to write 

|0> = |6 W ) ® \ X ), (H.31) 

which emphasizes that the ground state projection to the x Fock space is arbitrary. This 
is exactly the result of Section 5.2.1. 



H.2. Self energy and 0(1/N C ) effective Hamiltonian 
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H.2 Self energy and 0(1/N C ) effective Hamiltonian 

We now calculate an effective Hamiltonian that splits the spectrum of the degenerate 
subspace discussed above and removes the local degeneracy. We use Rayleigh-Schrodinger 
perturbation theory [67] where H± and if 3 are perturbations to Hi- In fact contributes 
at first order and if 3 at second order, and both give a contribution of 0(1/N C ). The 
effective Hamiltonian is therefore 

H cS = PH,P + PH 3 ^Lh 3 P = + H®. (H.32) 

Here P is a projection operator to the degenerate sector The projection operator is 

P=\b n ){O v \®l x , (H.33) 

and D is the energy denominator. H c g is an effective Hamiltonian for the x fields alone. 
Any bilinear of the 7r fields in is replaced by the bilinear's vacuum expectations 

value (vev), which we calculate using (H.28) and (H.30), 

(b\n^\b) N = ^E^ fc) (o|4Xlo) 

= J- E z lN[q ~ k) sinh ^ sinh ^(0|6_ fe 6L q |0) 

= ^E sinh V 9 , (H.34) 

1 ^ s q 
s k,q 

= J" E e^-*/ 2 ^- fe ) sinh y?_ fc sinh ^_ q <0| a_ fc aL g |0) 

s fc,q 



= E Sillh2 W 

- iv s qr 



(H.35) 



k^W6> = ^E^ (,J+fc)+J9( ^ i/2) (o|^|o) 



J_ ^ e «v(,+fe)+wW-*/2) cogh ^ ginh V9 _ g (o|a fc a t _ g |6) 

s fc,q 
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-^-^sinh^cosh^e^-^, (H.36) 

* »s q 



k,q 

= Y E e-™^-^-*™ sinh y? fc cosh V _ q {Q\b_ k b\ |6> 

s fc,q 

= sinh cosh ^,e- <9( ^ /2) . (H.37) 

Here we have used the following notation. N is a site on the fee lattice of even sites. On 
each fee site one has two degrees of freedom, 7r e , and it°. The neighboring sites of N are 
N + p, with 

p = 0,z,z±x/2,z±y/2. (H.38) 

Inserting Eqs. (H.34)-(H.37) into PH 4 P we have 

PH 4 P = |0 ff )(0 w |fZ- 4 |0 ff )(0 w | (H.39) 

Np I ij 



E (xSwXak + ,i) T7- E sinh2 



1 



2 

(xtei + xR+^XAT+fti) tt E sinh cosh ip^-m + ^. c . 



^ 1 jy E | 2 E XpjXkiXq+k-p,iXqj ' 7fc-p 



-2d (xfciXfci + XkiXli) x (sinh 2 ^ - 7<7 sinh y9 q cosh y9 q ) } . (H.40) 

Here we have omitted an irrelevant constant that represents terms quartic in tt. Finally 
we use Eq. (H.29) to write as 



+2 E kSx^U-pX- ' 7fc ^ \ (H - 41) 



fcpq 
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We see that has harmonic and anharmonic terms. In order to evaluate we write 
(1 - P) as 

1 - P = l x <g> \k a )(k a \ + \h)(k b \, (H.42) 

k 

where, for example \k a ) = a' k \0 n ) ® \x)- Here \x) is |fc a )'s component in the x Fock space. 
In general (H.42) is incorrect since it takes into account excitations of one spin wave only. 
Here it suffices since if 3 connects \0 n ) with excitations of that sort only. Next we replace 
the energy denominator with 

D = - 72, (H.43) 

and obtain the following expression for if eff ' , 

He's = ^ 2dJ N ^ T\ 2 ^ ^ (xn+p'iXNiKN+p ~ ^nX^n+p^Xni + h.c)j X 



|fca)(fca| + |^6)(fcfe| x ^2 (^N'+p,i'X e N'i' 7r N'+p' ^N'X^N'+p' ,i'XN'i' + h.C.^j P. 
N'p'i' 

(H.44) 

In order to calculate Eq. (H.44) we use Eq. (H.28) and note that 

(6 w |7rfe|fc a ) ~ <6 W | (a* or b) a[\0 n ) = 0. (H.45) 

On similar lines one can easily show that (O^lvr^lfca) = (0^|7r^|fc 6 ) = {O^n^lkb) = 0. 
The only non-zero elements are 

(6*1^1*0 = ^coshy9 fc e* fcJV , (H.46) 

(0,\^\k a ) = ^Asinh^e-^-^ 2 ), (H.47) 

(0„\rt\h) = ^smhip k e ikN , (H.48) 

(0M\k b ) = ^coshy^^/ 2 ). (H.49) 
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Using these relations in Eq. (H.44) we see that is a quartic interaction in the x fields. 
For example the contribution of \k b )(k b \ is 



iV s k yi-ii nn> 

v pp' ii' 

(cosh ^ k e^ N '-^\% t a% W - sinh Vk e~ ikN ' w) 



N s k 



(cosh (p k e tk ^ i/2) - sinh y? fc ) (cosh ip k e lk ^' 5/2) - sinh </? fc ) 

JV « fc NN> 

pp 1 a 1 

(coshy^e*^/ 2 ) - sinhy9 fc ) (cosh ip^-^-'l^ - sinhy? fc ) 



-</i^E cosh V?fc E XN+p,iX°N+p',i x (coshy^e^ £/2) -sinhy9 fc ) 

s k Npp'i 

(cosh ^ fe e- ife ^'- 5 / 2 ^ - sinh Vk ) . (H.50) 
We write the last (quadratic) term as 

-^E cosh2 ^Ex!X X ( cosh2 ^e i{k+q) ^ ) + sinh 2 ^ k+ ^~^ 

s qki pp' 

- COsh if k Sinh ^(k-lKp-zm+ilip'-zm + e i( q -fe)(p'-5/2)-^-5/2)^ 

= -4d 2 Ji-^- 2 xJiXqi cosh 2v? fc (7 2 +fc (cosh 2 ip k + sinh 2 y9 fc ) - 2-f q -f q+k cosh y9 fc sinh y? fc ) 

s qfcj 

= -4rf 2 Ji^^XqiXqi cosh 2v? fc (7 q +fc cosh 2v? fc -7 9 sinh 2v? fc )7 9 +fc (H.51) 

s qki 

In order to treat the quartic interactions in H e g we make an ansatz for the ground state 
of the x Fock space, and assume the following vacuum expectation values of the x fields. 

(6|xSx e fe 'j6> = (0|xgx^|0) = lA lk 5 kk Aj, (H.52) 



2 

(0\xliX O - k ' d \0) = (Olxl-X^lO)* = -l^kS^j, (H.53) 
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or in real space 

(0\x%X%>M = (0| X gx^|0) = £ A lk e ik ^'-^ • (H.54) 



(H.55) 



Here A 12 are assumed to be real and even in the momentum k. Note that a non-zero 
value of A 2 breaks the stationary group of the condensates of Q e and Q Q , U (2m — N) x 
U(N - m) x U(N - m), to 0(2m — N) x U(N - m). 

As in the treatment of order from disorder in condensed matter (see for example [52]) 
we decouple all the quartic terms in H e g to be quadratic terms times a condensate in all 
possible ways. The quartic term in becomes the following four quadratic terms. (We 
go back to real space.) 

^ X (XjXi) N (XiXj) N a — Ji E (x)Xi)n ixlxj) N d + (xjXi) N (XiXj}N+p 

Np 1 Np 1 F 



13 IJ 



~ Jr XNiX°N+p,j(XNjXN+p,i) + (XNiX°N+p,j)XNjXN+ 



p,l 



~ ^ ]\f 51 [^Q {x*N+p,iX°N+p,i + X^nXni) A 2q e lq ^ Z l' 2) (xNiX°N+p,i + X^NjX^N+p,i) 
s qNpi 

= 2dJ i^ X [ A i, + - A 2g79 - fc (xl t X°- k , t + • (H.56) 

s qrfc,i 

The next step is to use the ansatz (H.52)-(H.53) in Eq. (H.56). This gives 



frW _ j . 2 v ( fc v / A lfl + 1 - ^/l - Tf -A 2g7q _ fc \ / ^ 

^eflf - J 1«TT 2^ I Xfci X-fci J a Ail R 2 v to 

^ 9fcli V 7 V -A 2(? 7q-fc A lq + 1 - Jl - 7 | / V x_ fcj 



The same procedure is performed to the quartic terms of H^g ■ Here the procedure in 
much more cumbersome, and the four quartic terms of Eq. (H.50) alone give 4 x 4 = 16 
different quadratic terms. The |fe )(fe a | terms contribute 16 more. For example, in order 
to calculate the contribution of Eq. (H.50) we first decouple the quartic term, 

X^N+p,iX N'+p',i'X^N'i'XNi ~ ¥ (x\v+p,iX < N'+p',i')X^N'i'XNi + x\v+p,iX < N'+p',i'{X^N'i'XNi) 

Jr (XN+p,iXN'i')X°N'+p / ,i'Xm + X^N+p,iX^NU'(X°N'+p' ,i'XNi) ■ 

(H.58) 
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Next we substitute (H.54)-(H.55) into (H.58) and into its \k a )(k a \ counterpart, and find 
(after a lot of algebra) 



with 



Jiu = / BZ (J) (Ai, + l)7i( 9 ,fe)-i7(fc), (H.60) 

J 2fc = / Bz (g^A 2q / 2 ( Q ,fc). (H.61) 

I 1;2 (q,k) and r/(fc) are as given in Eqs. (G.29)-(G.30). The Hamiltonian (H.59) can be 
diagonalized by a Bogoliubov transformation 



X% \ _ ( coah\e q sinh±0 q W c q , 

v x% ) ~ [ sinh \e q cosh \e q ) [ dl q ) ' {ti - b2) 

with 



tanh(# fc ) = -J 2h /J lk . (H.63) 

The result is 

H cS = Jxd £ ^Jl k -Jl 2 (clc fc + dld fc ) , (H.64) 

which means that the energies of the x bosons are smaller by a factor of N c than the 7r 
fields' energies. 

To obtain self-consistent equations we calculate the vevs (H.52)-(H.53) using (H.59) and 
as in the calculation of the 7r vevs (H.34)-(H.35) one has 

(0\x[ e X q \0) = (0|x t 9 V 9 |0) = sinh 2 ^ <J = i(cosh^-l), (H.65) 

(6|x q X°- q |6> = (<6|x?xt°,|6>) = sinh U q cosh \o q = \ sinh 6 q , (H.66) 

which leads us to identify A^, = cosh#fc — 1, and A 2 k = — sinh 0%.. Substituting this into 
Eqs. (H.60)-(H.61) we write Eq. (H.63) as 

/ (%L) I 2 (q,k) smh0 q 

tanh£ fc = ' BZ ; ^ > , (H.67) 

L S h( q ,k)cosh0 q - v (k) 
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which is exactly the integral equation Eq. (5.52) for iV — m = 1. 

This Hamiltonian calculation sheds light on the meaning of the function 9 k that appears 
in Chapter 5 in terms of vevs of the x operators. A stable solution that obeys (5.49) 
has non-zero vevs. Mathematically, the integral equation (5.52) has also a trivial solution 
with 6k = that leads to A 12 = 0. However, this solution does not obey (5.49) and is 
thus unphysical, giving J7i < and making H c s not positive definite. 

Finally, we note that the Holstein-Primakoff transformation for iV — m > 1 is much more 
complicated [21], making the Hamiltonian approach unattractive for this case. The path 
integral approach is more general. 



Appendix I 



Effective ferromagnetic interaction 
at 0(1 /N c ) between odd sites in the 
nonlinear sigma model with non-zero 
density 



In this section we derive an effective interaction in 0(1/N C ) between the classical parts of 
the sigma fields on the odd sites, and show that this is a ferromagnetic interaction which 
prefers that these fields align to the same direction within the degenerate submanifold 
discussed in Chapter 4. The effective interaction comes from integrating over over fluc- 
tuations around the classical N c = oo ground state and derive an effective action for the 
classical part (4.95) of the odd spins. We let the a n on the even sites fluctuate around A; 
we let the a n on odd sites roll freely around the U(m)/[U(2m — N) x U(N — to)] man- 
ifold covered by Eq. (4.95), and also execute small oscillations off the manifold into the 
U(N)/[U(m) x U(N — to)] coset space. The counterpart of the action (5.1) for our prob- 
lem has an antiferromagnetic spin-spin interaction, with no (— l) n factors. We separate 
it into odd, even, and coupled terms, 



S' 



'even 



•odd 



s 




neven 
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iAF 



/*EW (i.i) 



Here a° n = Y, m (n) a m, where m(n) are the nearest neighbors of the even site n. We 
expand the field on the even sites around a n = A in the manner of Eq. (5.4), 

1 

a n = A + L n - -L 2 A (n even), (1.2) 
while for the odd sites we write (see Appendix J) 

cr n = U n ( l2m ~ N _ A , _ iL , 2A , J Ul (n odd), (1.3) 



with 

— In-tu 



A'= \~ m ° (1.4) 



and 

/ 77( m ) \ 

u " = ("o J-J- < L5 > 

L„ describes the fluctuations of the even spins around their classical value A. U n ro- 
tates the odd spins within the manifold of their classical values, while L' n describes their 
fluctuations outside that manifold. We further define 

•S = U. ( V —A' ) * = ( 1" lL ) ■ CL6) 
the classical field on the odd sites. 

We leave S odd alone and expand S cvcn and S AF around the classical values of the fields, 
S cvcn = -\fdrJ2 Tr AL n d T L n , (1.7) 

rieven 

S AF = So + y frfr^ (Tr L„oVj, - ^Tr L 2 n Aa n - Tr L n L n ) 

rieven ^ ' 

+rf Jx [drJ2 ("Tr AL„. + ^Tr iA . (1.8) 

nodd 



Here 



- = f» ( n r° ) ^ (1-9) 



L' n 
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is the rotated fluctuation field on the odd sites, and the Hermitian matrices a n and L r 
are sums over the odd neighbors of the even site n, 



m(n) 



cl 

mi 



m(n) 



(1.10) 

(1.11) 



Since both a n and A are block diagonal, the first trace in each integral in Eq. (1.8) is zero. 
Now we organize the partition function as follows: 

Nr. f^AA „ dJi 

2 

Vnodd / L \ " nodd 



Z = J ( II d *») exp [-^ (s odd + S + ^ / rfr £ Tr L 



(1.12) 



with 



JJ dL n j exp 
1 



V TIC VCIl 



4/* 



neven 



: Tr AL n d T L n 



^-Ty L 2 n Aa n - ^-Ty L n Z 7 



(1.13) 



Z even is a product of decoupled single-site integrals. Again we expand in the group algebra, 

L n = llM\ (1.14) 

where the sum is over the 2m{N — m) generators of U(N) that are not in U(m) x [/ (N—m). 
L n can be expanded similarly and we obtain the following form for the integral over the 
even fields: 



where 



Z cvcn = J Dl n exp -^Jd^ (%Mgl£ + 

The matrix C is the same as in Eq. (5.7), 

C-W = Tr (A [AT 7 , M r '']), 



(1.15) 
(1.16) 

(1.17) 



while the new matrix D varies with the site n according to the average a n of its neigh- 
boring spins, 

D W = -TY({M\M^'}ka n ). (1.18) 
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We study the two matrices in Appendix E. Diagonalizing C as before, we arrive at 

Zcvcn = I] {/ Dl eX P - Y l ^ + ^ + } , (1-19) 



where M. n is the m x m matrix 



M n = - (d r + .hE n ) . 



(1.20) 



The quantities l q (and l q ) for each q = 1, . . . , N — m are complex m-component vectors; 
they are rotations of the 2m(N — m) real components l v (and l v ) into the basis that 
diagonalizes C. The matrix E n is given by Eq. (6.11); it carries the dependence on a n . 
Performing the gaussian integration we have 



n 

nq 



1 



Det M, 



cxp 



NcfJi 
2 V4 



(1.21) 



Finally we separate the integral (1.12) over the odd spins into an integral over the classical 
field <7^ and an integral over the fluctuations around it. We obtain 



= J Dat exp (s + (N - m) £ Tr log M^j 



x J Dl m exp { S' 



odd 



^1 |/<? |2 _ 
r> rml 

^ TTK| 



(1.22) 



Here m stands for an odd site, n for an even one. 

Equation (1.22) gives an effective action for the classical odd spins a^. These enter the 
exponents through M. [via Eqs. (6.11) and (1.20)] and through [via Eqs. (1.9) and (1.11)]. 
The action in the first exponent is minimized when each matrix E^af^ has the largest 
number of zero eigenvalues, each of which makes Tr \ogM. n approach — oo. It is easy to 
check that E n has 2m — iV zero eigenvalues (the maximal number) when the on all 
the odd sites m(n) align, i.e., 



= a E U(m)/[U(2m - N) x U(N - mj\. 



(1.23) 



Moreover, when Eq. (1.23) holds, all the Z^'s align parallel to each other and is maxi- 
mized; also the eigenvalues of are maximized (to +oo). Thus the action in the second 
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exponent also has a minimum at this point in configuration space. These effects add up 
to an effective ferromagnetic interaction among the 2d nearest neighbors m of each even 
site n. This effective interaction will align the classical spins on the odd sublattice to the 
same direction in their submanifold, U(m)/[U(2m — N) x U(N — m)]. 

The divergences in the effective action have their origin in the fact that the semiclassical 
corrections are calculated as gaussian integrals in the even fluctuation fields L n , and the 
coefficient matrix M. n acquires zero eigenvalues. The correct range of integration over 
L n is of course not infinite, but rather the volume of the U(N)/[U(m) x U(N — m)] 
manifold. This will regulate the divergences, but leave the effective action for the odd 
spins attractive. 



Appendix J 



Fluctuations on the odd sites at 
non-zero density 



In this section we show how to separate the fluctuation outside the degenerate submanifold 
of the odd spin. 

In the classical analysis, the fields on the odd sites take values in the sub-manifold 
U{m)/[U(2m - N) x U(N - m)\ of the manifold U(N)/[U(m) x U(N - mj\. We de- 
note these values <r cl , 

'■' ' 1 (J.l) 



Here 

with e U(m), and 



V IjV-m / 
a (m) _ jj(m) _^(m)^y(m)t_ 



1-2m-N 




-N-m 



(J.2) 
(J.3) 



(jM contains 2(N —m) (2m — JV) independent degrees of freedom. Any a e U (N)/[U (m) x 
U(N — m)] can be written as 



cos 



f2v / aat N ) 

V > Va'a 

sin (2v / ata) 



a} a 



cos 



sin (2 Vat a) \ 

V at a 

(2v / at^) 



(J.4) 
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[cf. Eq. (4.19)], which coincides with Eq. (J.l) if 



a = U {m) 





(tt/2)1 



N-m 



(J.5) 



Recall that a is an to x (N — to) matrix, so the zero block has dimensions (2m — N) x 
(N-m). 



We allow motion out of the sub-manifold by allowing a to vary further, 



a = £/ (m) 







(J.6) 



The 2(N-m) 2 de grees of freedom in a complement the 2(N — m)(2m- N) degrees of free- 
dom inherent in to give 2m(N — to), the dimensionality of the entire U(N)/[U(m) x 
U(N — to)] coset space. 



Writing a with the generalized a we have 



a = U 








cos 
sin 



(2 v / ^t) 



—a 



sin (2V a^aj 



with 



We can also write this as 



U 



a = U 



> Vfl'fl 

in(2vata) , / _ 

— y _ ' qt -cos (2^ 

V at a v 



' f/ (m) \ 

\ IjV-m / ' 
hm-JV \ j- 

aP^— )l(a,at) r • 



(J.7) 



(J.8) 



(J.9) 



a [2(N m)} j g & ma ^ r i x j n the man if ld U(2(N — m))/[U(N — m) x U(N — to)]. Indeed for 
a = (7r/2)l A r_ m , we have 



r [2(JV-m)] 



"ljV— m 
ljV-m 



-A. 



(J.10) 



Since U(2(N — m))/[U(N — to) x U(N — to)] is a self-conjugate manifold, its structure 
near o"[ 2 ( Ar_m )] = —A is the same as its structure near o"' 2 ^ - " 1 ^ = A, which corresponds 
to a = 0. Expanding o"[ 2 ( JV-m )] around —A gives Eq. (1.3). 



Appendix K 

Coefficients in the dispersion 
relations of Goldstone bosons 



Here we present the various coefficients that appear in the dispersion relations given by 
Eq. (6.62). These are c, ci, a, and b found in Eqs. (6.63)-(6.68). These coefficients are 
extracted from the nearest-neighbor theory, and depend on the number N — m only. 

The constants c and c\ are the zero momentum curvatures of the following functions, 



a and b are defined as 



o = i i , (k.i) 

fc=0 




a = - -rrrk ■ (K.2) 



N c? - c 2 

(K.3) 



6Ji 2c lC 2 ' 

(K.4) 

6 Ji ci V ; 



We present the constants in Table K.I. 
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Table K.l: The constants c, ci, a, and b, for all values of N—m. Recall that for N—m = 1, 
the solution of the integral equation is different from the other cases, and it is possible 
that the nearest-neighbor spectrum follows a power of momentum higher than 2. This is 
the source of the sizeable differences between this case and the other cases. 



N — m 


c 


Cl 


a 


b 


1 


8.3 x 1(T 5 


0.12 


8.4 x 10 6 


17.1 


2 


0.002 


0.01 


1.2 x 10 3 


187 


3 


0.003 


0.009 


490 


220 


4 


0.003 


0.008 


350 


240 


5 


0.003 


0.008 


280 


250 



Here we note that the case of iV — m = 1 has an exceptionally large value for a, with a 
very small curvature c. This is a result of the different behavior of the solution to the 
integral equation for this case [see Fig. 5.4]. This behavior should disappear if we take a 
better momentum mesh, especially near k — 0. 



Appendix L 



Classification of fields according to 
unbroken subgroups 



In this section we consider the a fields that represent Goldstone bosons of the effective 
theory. The a multiplet contains fields called n that are type I Goldstone bosons in the 
nearest-neighbor theory, and fields called x that are type II bosons there. Moreover, the 
antiferromagnetic nature of the theory leads us to consider separately the fields on the 
even and odd sites of the lattice. 

The next-nearest-neighbor theory breaks its chiral SU {Nf) x SU(Nf) symmetry sponta- 
neously as shown in Table 6.1. The following two tables show the classification of the 
various fields according to their representations and charges under the unbroken groups. 
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Table L.l: Classification of a fields on the even sites under the unbroken symmetry for 
3iV/4 < m < N. As indicated in Eq. (6.38), not all x fields are present for a given 
combination of Nf and B. We neglect U(1)a- All fields are complex. 

Symmetry rep 

Nf B Symmetry Field Matrix structure or U(l) charge (s) 



1 


1 







1 X 1 






2 


2 


SU{2) 


7T 


2x2 




1 ©3 








Xl 


4x2 


1® 


1 ©3 


©3 


2 


3 


u(i) h 


7T 


1 x 1 




+ 1 














f +1 ^ 












4 x 1 















Xle 




+ 1 
















I o J 












2 x 1 














X2e 




[") 
















3 


3 


SU(3) 


7T 


3x3 




108 








Xl 


9x2 


1© 


1©8 


©8 



3 4 1/(1) /3 x U(1) Y vr 



Xu 



X2e 



2x2 



6x2 



2x2 



(+1/2,-3) 
(+1,0) 
(+1,0) 
(0,0) 
(+1/2,-3) 
(+1,0) 
(0,0) 
V (+V2,-3) 



(0,0) 
(+1/2, +3) 
(+1/2, +3) \ 
(-1/2, +3) 

(0,0) 
(+1/2, +3) 
(-1/2, +3) 

(0,0) 



(+1,0) (+1/2, +3) 
(0,0) (-1/2, +3) 



5 U(l) h xU(l) 



Y 



7T 



Xle 



lxl 



6 x 1 



-1/2, +3) 
+1/2, +3) \ 
-1/2, +3) 

(0,0) 
+1/2, +3) 
-1/2, +3) 

(0,0) 
+1/2, +3) 
-1/2, +3) 
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Table L.2: Classification of x fields on the odd sites. The classification of the odd tt fields 
is the same as of the even tt fields for all values of Nf and B. For Nf — B — 1, 2, 3, the 
classification of the odd x fields is the same as of the even x fields. Here we present the 
classification of the odd x's for the other cases. All fields are complex. 



Nf B Symmetry Field Matrix structure 



U(l) charge (s) 



2 3 



U(l) 



h 



Xlo 



4 x 1 



X2o 



2 x 1 



( ^ 

-1 


v- 1 / 



-1 



3 4 U(l) l3 x U(1) Y xio 



6x2 



X2o 
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